Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4940TG52.8% 29.8 / 51.1 36intergenic (‑512/+91)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (3/14);  new base G (19/0);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.79e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-05

CAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987351‑987644
                                                                                                                                               |                                                                                                                                                      
cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaaa                                                                                                                                                  >  1:832637/1‑150 (MQ=255)
 aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                                                                                                <  2:939455/151‑1 (MQ=255)
 aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTaaa                                                                                                                                                  >  2:863287/1‑149 (MQ=255)
              gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACt                                                                                                                                    <  2:904224/150‑1 (MQ=255)
               ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAAttacgt                                                                                                                                        >  2:455056/1‑141 (MQ=255)
               ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTTCTTTACtt                                                                                                                                   >  2:952387/1‑150 (MQ=255)
                  aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTCTTGACtt                                                                                                                                   >  1:949689/1‑147 (MQ=255)
                    aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaa                                                                                                                             <  1:622133/151‑1 (MQ=255)
                      aTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAtttgttttctttaaaa                                                                                                                              >  1:468809/1‑134 (MQ=255)
                                      gCTATTAGCTAAGAAAAAACTAATCCGTAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaa                                                                                                           <  2:255402/151‑1 (MQ=255)
                                       cTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCGAAAAGATTGCTTGTCTTTAAAAAAAAATaaaatgaaag                                                                                                              >  1:941178/1‑139 (MQ=255)
                                         aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGa                                                                                                        <  1:900820/151‑1 (MQ=255)
                                                aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                                                                                                   <  1:26837/149‑1 (MQ=255)
                                                   aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGAATTTTTGGCTTTACAAACAAAAAAATTTGCAAAAAAAATaaaaaaa                                                                                                >  1:902615/1‑149 (MQ=255)
                                                        aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCCTGACTTAACAAACAGATAACTTGCCAGAAAAGATAAAAAAAGAGCCAg                                                                                         <  1:597859/151‑1 (MQ=255)
                                                          tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAAATTATTGACTTAAAAAAAAGATAAAATTACAGAAAAGAaaaaaaaaaaaac                                                                                           >  1:689012/1‑142 (MQ=255)
                                                              ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                    <  2:606706/150‑1 (MQ=255)
                                                                gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGCTTGACTTAAAAAAAAGATAAATTGACAGAAAAAATaaaaaaaaa                                                                                              >  2:88590/1‑136 (MQ=255)
                                                                           gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                       <  1:565196/150‑1 (MQ=255)
                                                                           gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAAAAAATTGCTTTACTTAACAAAAAAATAAATTTACAGAAAAAAAAAAAAAAAAACCAGCGGGGGCTGGCTCtt                                                                          >  2:515951/1‑147 (MQ=255)
                                                                           gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTTTTTGGATTAACAAAAAAATAACTTTACAGAAAAGATAAAAAAAAAAGCAGCGGGGGCTGGCTCTTTaa                                                                       >  2:287149/1‑150 (MQ=255)
                                                                            cGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                     <  1:785414/151‑1 (MQ=255)
                                                                               aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGc                                                                  <  2:222367/151‑1 (MQ=255)
                                                                                 tAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                                 <  1:431850/150‑1 (MQ=255)
                                                                                           gcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGCGAATAAATTTGTTTTCTTTACAAAAAAAAAAATTTAAAGAAAAAATAAAAAAAAAACCAGCGGGGGGCGGGTCTTTAAAAAGCtaaaataaaaa                                                       >  1:555313/1‑149 (MQ=255)
                                                                                             gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGCTTGACTTTACAAACAAATAACTTGACGGGAAAAATAAAAAAAAAAGAAACGGGCGCGGGGCCCTGAAAAAGCCAAATTATaaaatg                                                    >  2:225884/1‑147 (MQ=255)
                                                                                              tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTCTTGAATTAAAAAACAGATAAATTGACAGAAAAAAAAAAAAAAAAAACAGCGGGGGGCGGGTCTTGAAGACGCTAAAATAAAAACTgg                                                   >  2:622571/1‑151 (MQ=255)
                                                                                                   gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTGCTTGACTTTACAAAAAAATAACTTTACAGAAAAGAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCggt                                               >  2:440256/1‑148 (MQ=255)
                                                                                                      gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTCTTTTCTTTAAAAAAAGATAAATTGACAGAAAAAAAAAAAAAAAAGACAGCGGGCGCTGGCTCTTTAAAACGCTAAAATAAAAACTGGCGTTAcg                                            >  1:371501/1‑150 (MQ=255)
                                                                                                          tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCGAAAAAATTGCTTTACTTAAAAAAAAAAAAACTTTACAGAAAAAGAAAAAAAAAAGCAAACGGGGGGCGGGTCTTGGAGAAGACAAAATAAAAACTGGAGGTACGCgg                                         >  1:6842/1‑149 (MQ=255)
                                                                                                                atTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTTTTGACTTTACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAACGGGCGCGGGCTCTTTTAGAAGCCAAAATAAAAAATGGCGGTTCGCGGAGGAc                                    >  1:892715/1‑148 (MQ=255)
                                                                                                                     gTAACAGGATGATCGTTATGGCAGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAATACGCTAAATTAGAAGCCGGCGTTACGCGGAGTACGTGGGaa                             <  2:624865/150‑1 (MQ=255)
                                                                                                                      tAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTTGTTTACTTAAAAAAAAAATAAATTGACAGAAAAGATAAAAAAAAAAACCAAGGGGGCTGGGTCTTTAAGACGCTAAAATAAAAACAGGCGGTACGCGGAGGACGTgggga                             >  1:21126/1‑147 (MQ=255)
                                                                                                                                       tGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacg          >  2:921430/1‑151 (MQ=255)
                                                                                                                                          cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt       >  2:348190/1‑151 (MQ=255)
                                                                                                                                              gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa   <  2:168586/151‑1 (MQ=255)
                                                                                                                                               ttGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc  >  2:903500/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987351‑987644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑