Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 982,566 C→A 23.2% intergenic (‑128/+475) 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,5660CA23.2% 57.6 / 11.4 35intergenic (‑128/+475)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base A (3/5);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTT  >  CP009273/982420‑982712
                                                                                                                                                  |                                                                                                                                                  
gAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCaa                                                                                                                                                  <  1:652539/149‑1 (MQ=255)
              atcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCt                                                                                                                                    >  1:691314/1‑149 (MQ=255)
               tcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTa                                                                                                                                   >  2:363981/1‑149 (MQ=255)
                       tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAGCGTTTCGCTACCGTGaaa                                                                                                                           <  1:560532/149‑1 (MQ=255)
                       tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGaaa                                                                                                                           <  2:292169/149‑1 (MQ=255)
                           aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGc                                                                                                                       <  2:552425/149‑1 (MQ=255)
                               tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                                                                                                                   >  2:212816/1‑149 (MQ=255)
                                    gTTTTGTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGCTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagat                                                                                                              <  2:558627/149‑1 (MQ=255)
                                       ttttttATGACTCCTGCCACTGCCGTCAATACGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTATGCTGGCATCTTGCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                                                                                                           <  2:386657/149‑1 (MQ=255)
                                       ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                                                                                                           <  1:84304/149‑1 (MQ=255)
                                           ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAtt                                                                                                       >  1:542817/1‑149 (MQ=255)
                                              tGACACCTGCCACTGCCGTCAAGAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTgg                                                                                                    <  2:781639/149‑1 (MQ=255)
                                                  acCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc                                                                                                >  2:94923/1‑149 (MQ=255)
                                                    cTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGGTAGGGGTATTTATGTGTCTGCAGGCGTCTTTCCTTTCAAACTAACGTTTCGCTACCGTGAAAGGGCTACCAAGATACAGGATTAGTTTCaa                                                                                              >  2:141433/1‑149 (MQ=255)
                                                            tgcGTCAATAAGTTCTGTTGATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa                                                                                      <  1:301821/147‑1 (MQ=255)
                                                            ggcGTCAATAAGTTTTGTTAATAATAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa                                                                                      <  2:295751/147‑1 (MQ=255)
                                                             ccGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTCCAAAGATAAAGATTTGGTTTAAAAAAGaaaaa                                                                                     >  1:59253/1‑149 (MQ=255)
                                                               gTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAat                                                                                   <  2:658616/149‑1 (MQ=255)
                                                                tCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCTAAAAGAAAAAatt                                                                                  >  2:430161/1‑148 (MQ=255)
                                                                 cAATAAGTTCTGACAATAAAAAATCACGAAACTATTAATGAGAGTTTGAAGTCTTTATATGTCTGCAGGCAACTTTCCAATTAAACAAACCTTACGCTACCGTGAAAGTTCTACAAAGATAAAGATTTGGAAAAAAAAAAAAAAAatat                                                                                 <  1:281014/149‑1 (MQ=255)
                                                                     aaGTTCAGTCAAAAAAAATTGACGCAACTATTGATGAGATTTTTGTTGCTTTATGTTTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAACTGCTACAAAGATAAAAACTTGCTTTCAAAAAGAATATATATGTaa                                                                             <  1:236498/149‑1 (MQ=255)
                                                                          cTGTTAATAAAAATTTACTGAACTATTGATGGGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCaaa                                                                        <  1:283490/149‑1 (MQ=255)
                                                                                                   ttGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa                                               >  2:806593/1‑149 (MQ=255)
                                                                                                   ttGATGAGAGTTTGGGGTTTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa                                               <  1:123164/149‑1 (MQ=255)
                                                                                                       ttagagGTTGGTGTCTGTATTTGGCTGCAGGCATTTTTCCATTAAAACTAACGTTTTGCTACCGTGAAAGGGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTATAATTTAAGGAACTTTGTGATCAcc                                           <  1:166592/147‑1 (MQ=255)
                                                                                                             tttGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGCACACCGTCata                                     >  1:286112/1‑149 (MQ=255)
                                                                                                             tttGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCata                                     >  1:310353/1‑149 (MQ=255)
                                                                                                                ggTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt                                  <  1:212816/149‑1 (MQ=255)
                                                                                                                                    ggCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGatat              <  1:363981/149‑1 (MQ=255)
                                                                                                                                     gCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt             >  2:631388/1‑149 (MQ=255)
                                                                                                                                     gCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt             >  1:786832/1‑149 (MQ=255)
                                                                                                                                     gCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTAAAAAGCTGAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt             >  2:817206/1‑149 (MQ=255)
                                                                                                                                                ttCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCGAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGAAGAGATGATATTTACAGCAAtta  >  2:840594/1‑148 (MQ=255)
                                                                                                                                                ttCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAAttt  >  2:773090/1‑149 (MQ=255)
                                                                                                                                                ttCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAAttt  >  1:52951/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTT  >  CP009273/982420‑982712

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