Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3580TG41.7% 55.1 / 10.0 36T66P (ACC→CCC) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (4/17);  new base G (15/0);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 7.21e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266214‑266506
                                                                                                                                                |                                                                                                                                                     
cctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGAcc                                                                                                                                                 <  1:4016937/151‑1 (MQ=255)
         ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAagta                                                                                                                                       <  2:2560978/151‑1 (MQ=255)
          cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGAAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGGAAAAACCGCACGGTGGGACCggggaaata                                                                                                                                        >  1:1331097/1‑143 (MQ=255)
            gTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCCGCAAAGGGGGACCGGG‑GAagtag                                                                                                                                      >  1:1468415/1‑149 (MQ=255)
                    ttATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAGGGGGGGCCGGG‑TAAGGAGGCaa                                                                                                                                  >  1:2022935/1‑145 (MQ=255)
                     tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCtgtg                                                                                                                            <  2:3210874/150‑1 (MQ=255)
                      aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGAc                                                                                                                          <  2:1413459/151‑1 (MQ=255)
                         ggCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGG‑AAAGGAGTCAAACTGTGACCa                                                                                                                        >  2:1077942/1‑150 (MQ=255)
                             aTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCAttctt                                                                                                                   <  1:1911143/151‑1 (MQ=255)
                               acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCAttcttc                                                                                                                  <  2:2013558/150‑1 (MQ=255)
                                     cgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCCAAGTGGGGCCGGA‑GAAGGAGTAAAACTGTGAACAATCTTTGCCGAc                                                                                                            >  1:1283069/1‑150 (MQ=255)
                                            atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCCAGGGGGGACCGGA‑TAAAAAGACAAACCGGGACCATTCTTTGCCGAca                                                                                                           >  1:2792481/1‑144 (MQ=255)
                                                  ggCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGCGGAAAACCCGCCAGGGGGGGCCGGAGAAAGAAGTCAACC‑GGGAGCCTTCTTTCCTGACCTCAACCtcgggg                                                                                               >  1:2271750/1‑145 (MQ=255)
                                                      ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTg                                                                                            <  1:1310480/149‑1 (MQ=255)
                                                      ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAGGGGGGGCCGGG‑TAAGGAGTCAAACTGGGACCATTCTTCTCTGACATAATCCtgg                                                                                                  >  2:1943340/1‑143 (MQ=255)
                                                        aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGTAAAAACTGCAAGGGGGGGCCGGA‑TAAGAAGTCAACCTGTGACCCTTCTTCGCTGACATCATCCTGGGGGCTGAcc                                                                                         >  2:2493800/1‑150 (MQ=18)
                                                             tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACTGCAAGGGGGGGGCGGG‑TAAAAAGTCAACCCGTGGCCATTCTTCTCTGACATCATCCCGGTGGGTGACcac                                                                                       >  1:3646766/1‑147 (MQ=255)
                                                              ttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCa                                                                                   <  1:2310683/150‑1 (MQ=255)
                                                                  ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAAGGGGGGCCGGG‑TAAATAGTCAAACTGTGGCCAATCTTCGCTGACATCATCCTGGGGGCTGGCCACCTCACTGa                                                                               >  2:2872224/1‑150 (MQ=255)
                                                                     gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAc                                                                            <  1:2170195/150‑1 (MQ=255)
                                                                      cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGACCCGG‑TAAGGAGTCAACCTGTGACCCTTCTTCCCTGACATCATCCTGGTGGCTGACCACATCACTGaaa                                                                             >  2:919341/1‑148 (MQ=255)
                                                                         atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGGGGGAACGGG‑GAAAGAGTCAAACTGTGACCATTCTTCGCTGACATCATCCTGGTGGGTGACCACATCACTGAAAAGCgg                                                                        >  2:349875/1‑150 (MQ=21)
                                                                                aGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACTGCAAGGGGGGGCCGGG‑TAAGTAGTCAAACTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGCAACGCGATgggggg                                                                 >  2:706672/1‑148 (MQ=21)
                                                                                       agtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCt                                                          <  2:862801/150‑1 (MQ=255)
                                                                                          aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGGGGGGCCGGC‑TAAATAGGCAAACCGTGACCCTTCTTCGCTGACATCATCCTGGTGGCTGGCCCCATTACTGAAACGCGGGGGGGGGAGGGCCtt                                                         >  2:511824/1‑148 (MQ=255)
                                                                                             tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGACCGGG‑TAAAAAGTCAACCTGTGACCATTCTTTGCTGACATCATCCTGGTGGCTGGCCCCATCACTGAAACGCGGGGGGGGGGGGGCCCTAAAttt                                                   >  2:3385181/1‑151 (MQ=255)
                                                                                               ccATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc                                                 <  2:1623514/151‑1 (MQ=255)
                                                                                                       gCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGcc                                         <  2:1383739/151‑1 (MQ=255)
                                                                                                               cGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtc                                 <  2:2639095/151‑1 (MQ=255)
                                                                                                                gggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAg                                    <  2:853537/147‑1 (MQ=255)
                                                                                                                gggCGACTGGCGGGTACGTAAAACCCGCAAGGGGGGACCGGA‑TAAGAAGGCAACCTGTGACCAATCTTCCCTGACATCATCCCGGGGGCTGACCACATCACTGAAAAGCCGTGGGGGGAGGGCCTTAAATTTACTGATGGGCTTCAgaga                                 >  1:853226/1‑147 (MQ=255)
                                                                                                                                 gTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc                >  1:1681428/1‑150 (MQ=255)
                                                                                                                                   aaaaCCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt              <  1:3329902/150‑1 (MQ=255)
                                                                                                                                    aaaCCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGc            <  1:4091633/151‑1 (MQ=255)
                                                                                                                                           cAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgaga     <  2:926219/151‑1 (MQ=255)
                                                                                                                                              ggTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc  >  2:88352/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGA‑TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266214‑266506

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