Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 982,565 T→A 100% intergenic (‑127/+476) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273982,5650TA72.7% 11.4 / 3.0 11intergenic (‑127/+476)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base A (3/5);  total (3/8)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01
Rejected as polymorphism: Not supported by required number of reads on both strands.

TGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACAC  >  CP009273/982426‑982670
                                                                                                                                           |                                                                                                         
tGCGCTTCATCAGTATTTATTACCCTCTTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATACaaa                                                                                                       <  1:51726/144‑1 (MQ=255)
                     aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCATGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGc                                                                             <  2:101452/149‑1 (MQ=255)
                           aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTATCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATACAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaa                                                                       >  2:2790/1‑149 (MQ=255)
                              gTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGGCTGCAGGCATCTTTTCATACAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagat                                                                    <  2:138902/149‑1 (MQ=255)
                                                gCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACTGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa                                                  <  2:72378/149‑1 (MQ=255)
                                                       ccGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAg                                                <  1:248563/144‑1 (MQ=255)
                                                                  tgttGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATACAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGt                                     <  1:2469/141‑1 (MQ=255)
                                                                       tCAATAAAAATTTACGGAACTATTGATGAGGGTTTGATGTCTTTATGTGTCTGCAGGCATCTTTCCATACAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGt                           >  2:210237/1‑149 (MQ=255)
                                                                                    aCGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAGCTTTCCATACAAACTAACGTTTCGCTACCGTTAAAGTGCTACAAAGATAAAGATTTTGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTa                   <  1:75322/144‑1 (MQ=255)
                                                                                        aaCTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATACAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTTATAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACttt          <  2:64611/149‑1 (MQ=255)
                                                                                                atgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATACAAAATAACAATTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac  >  2:215671/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                         
TGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACAC  >  CP009273/982426‑982670

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: