Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7610CT41.7% 46.7 / 18.5 36G56E (GGA→GAA) ralRRac prophage; restriction alleviation protein
Reads supporting (aligned to +/- strand):  ref base C (4/17);  new base T (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

CCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAA  >  NC_000913/1413616‑1413911
                                                                                                                                                 |                                                                                                                                                       
ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTg                                                                                                                                                    <  2:432288/151‑1 (MQ=255)
   gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACCCCTCCATTATTATTTCCAGTGGtt                                                                                                                                                 <  1:113580/151‑1 (MQ=255)
      tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACCCCCCTTTTTTTTTTCAGTgg                                                                                                                                                   >  1:544837/1‑146 (MQ=255)
      tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTc                                                                                                                                                <  2:380114/149‑1 (MQ=255)
             tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTATAAATTAAAACCCCCCTTATTATTTTCAGTGGTTCGtttttt                                                                                                                                         >  1:618588/1‑149 (MQ=255)
               ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAACAACCCCCTTATTTTTTTCAGTGGTTCTtttttt                                                                                                                                         >  2:472115/1‑147 (MQ=255)
               ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt                                                                                                                                     <  2:274862/151‑1 (MQ=255)
                     agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                                                                                                                               <  2:301764/151‑1 (MQ=255)
                          tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAAACTCCCTTTTTATTTTCAGTGGGTGGTTTTTTTCTTCTTTTGAGTg                                                                                                                           >  1:469362/1‑150 (MQ=255)
                          tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATTAACACCCCCATTATTATTTTCAGTGGGTTGTTTTTTTCATTTTTCgcg                                                                                                                             >  1:603074/1‑146 (MQ=255)
                           ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATAGAAACCTCCCTTTTTTTTTTTCGTGGGTCGGTTTTTTCCTCTTTTGAg                                                                                                                             >  1:624266/1‑147 (MQ=255)
                            tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACAACCCCATTATTATTTTCAGTGGGTTGTTTTTTTCCTCTTTTTAGTGgtt                                                                                                                        >  1:599734/1‑151 (MQ=255)
                                                 gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAAATGAACACCCCCATTATTATTTTCAGTGGTTTGTTTATTTCATCTTTCGAGGGGTTCTTTTTTACTTCCCCCa‑aaa                                                                                                   >  1:315388/1‑150 (MQ=255)
                                                   ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAAcc                                                                                                 <  1:583782/150‑1 (MQ=255)
                                                               tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAAACCCCCTTTTTTTTTTCAGTGGGTCGTTTTTTTCCTCTTTCGAGTGCTTTTTTTTTTCTTTCCCCA‑TAAACCGGTCggggt                                                                                       >  1:82806/1‑146 (MQ=255)
                                                                    tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACCCCCCTTATTTTTTTCCGTGGTTTGTTTTTTTCATCTTTTTAGTGCTTCTTTTTTTCTTTCACCC‑TAACCGGGTCGGGGTTTGCATc                                                                                >  1:217622/1‑150 (MQ=255)
                                                                     gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACACGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATcc                                                                               <  2:146978/150‑1 (MQ=255)
                                                                       ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACCCCCCCTTTTTTTTTCCGTGGGTTTTTTTTTTTTTCTTTTGAGGGCTTCTTTTTTTCTTTCCCCA‑TAACCCGGTCGGGGTTCGCCTTCGt                                                                             >  2:380352/1‑150 (MQ=255)
                                                                       ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCCCCATTATTTTTTTCAGTGGGTTGTTTATTTCATCTTTTGCGTGCTTCTTTTTTTCTTCCTCCA‑TAACCGGTTCTGGATTCGCCTCCGt                                                                             >  2:151887/1‑150 (MQ=255)
                                                                       ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGtt                                                                            <  1:128878/151‑1 (MQ=255)
                                                                            tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGAAATTTTTTTTAATAACAACCCCCCTTTTTTTTTTAGTGGTTTTTTTTTTTCCTTTTTTGAGTGCTTCTTTTTTTCTTCCTCCATTAACCGGTT‑GGGGTTCGCATCCCtt                                                                            >  2:514910/1‑146 (MQ=255)
                                                                              gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTACATAAAAACCTCCCTTATTATTTTCAGGGGGTCGGTTTTTTTTTCTTTTCTGTGCTTCTTTTTTTCTTTCACCA‑TAAACGGGTCTGGGGTCGCCTCCCtt                                                                            >  2:381031/1‑144 (MQ=255)
                                                                                          ccttccttTTCAATAGTGGCGGGAATTTTTTATATGAACACCTCCCTTTTTTTTTTCCGTGGTTCGTTTTTTTCCTCTTTTGCGTGGTTCTTTTTTTCTTCCACCC‑TAAACGGGTCGGGGTTTGCATCCGTTACACTTCGCTccg                                                               >  1:402823/1‑143 (MQ=255)
                                                                                                 tttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaa                                                   <  1:465331/150‑1 (MQ=255)
                                                                                                  ttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaaa                                                  <  2:507967/150‑1 (MQ=255)
                                                                                                  ttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaaa                                                  <  2:266617/150‑1 (MQ=255)
                                                                                                      aTAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCt                                             <  1:503580/151‑1 (MQ=255)
                                                                                                               gTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt                                    >  1:502261/1‑151 (MQ=255)
                                                                                                                     ttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt                              >  1:607236/1‑151 (MQ=255)
                                                                                                                        taCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg                            >  1:516408/1‑150 (MQ=255)
                                                                                                                             gAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAAcc                      <  1:398701/151‑1 (MQ=255)
                                                                                                                             gAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAAcc                      <  1:138969/151‑1 (MQ=255)
                                                                                                                              aaCACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACca                     <  1:533407/151‑1 (MQ=255)
                                                                                                                                cacCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACcaca                   <  1:347686/151‑1 (MQ=255)
                                                                                                                                            tattTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATg        <  1:429123/150‑1 (MQ=255)
                                                                                                                                             attTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATgg       >  2:181044/1‑150 (MQ=255)
                                                                                                                                                 ccAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTaa  <  2:48670/151‑1 (MQ=255)
                                                                                                                                                 ccAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTa   >  2:400345/1‑150 (MQ=255)
                                                                                                                                                 ccAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTa   >  2:581101/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAA  >  NC_000913/1413616‑1413911

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