Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3490CA44.4% 99.5 / 31.4 55V69F (GTT→TTT) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (4/26);  new base A (24/0);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT  >  NC_000913/266204‑266499
                                                                                                                                                   |                                                                                                                                                      
tCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACTGc                                                                                                                                                    >  1:347162/1‑150 (MQ=255)
     cACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTg                                                                                                                                              <  1:654718/151‑1 (MQ=255)
     cACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAgggg                                                                                                                                              >  1:136460/1‑149 (MQ=255)
            tGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCCAAGGGGGGCCgg                                                                                                                                       >  2:468891/1‑151 (MQ=255)
              gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGa                                                                                                                                      <  2:735653/150‑1 (MQ=255)
              gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGCGTAAAAACCGCAAGGTGGGGCCgg                                                                                                                                       >  1:368760/1‑149 (MQ=255)
                aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGATa                                                                                                                                    >  1:532513/1‑150 (MQ=255)
                  tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAAGGGGGGCCgg                                                                                                                                       >  2:777144/1‑145 (MQ=255)
                   ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAagt                                                                                                                                 <  2:150616/150‑1 (MQ=255)
                   ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAagt                                                                                                                                 <  2:945274/150‑1 (MQ=255)
                    cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCCGCAAGGGGGGGCCggggaa                                                                                                                                   >  1:449654/1‑143 (MQ=255)
                       tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGTAAAAGTTGCGGGGGGGGTCCggcgagcg                                                                                                                                 >  2:397377/1‑141 (MQ=255)
                         tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                           <  1:288156/150‑1 (MQ=255)
                         tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                           <  1:840724/150‑1 (MQ=255)
                               tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCtgtg                                                                                                                     <  1:94214/150‑1 (MQ=255)
                                  cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCa                                                                                                                 <  1:257814/151‑1 (MQ=255)
                                  cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGCCGGGTAAGGAGGCAAACTGGGaac                                                                                                                  >  2:612510/1‑148 (MQ=255)
                                     cAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGATAAAGAGTCAACCTGTGACCaat                                                                                                               >  1:495315/1‑148 (MQ=255)
                                      aaTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttct                                                                                                             <  2:669612/151‑1 (MQ=255)
                                                 cATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAcatc                                                                                                  <  1:202315/151‑1 (MQ=255)
                                                    tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACCGCCAGGGGGGACCCGGTAAAAAGTCAAACCGTGGCCCATCTTCGCCGAcatca                                                                                                 >  1:657127/1‑149 (MQ=255)
                                                      atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCt                                                                                             <  2:619993/151‑1 (MQ=255)
                                                               tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGAGGAAAAACAGCAAGGCGGGGACGGGCAAGGAGTGCACCTGTGCCCATTCTTTACTGTCATAATCATGGGGGCTg                                                                                     >  2:829255/1‑150 (MQ=255)
                                                                ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGGCCGGGTAAGTAGTCAACCCGTGACCATTCTTTGCTGACATCATCCTGGGGGGTg                                                                                     >  2:783671/1‑149 (MQ=255)
                                                                      ttttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTTCAAAGTGGGGCCGGACAAAAAGGCAACCCTTGTCCATTCTTTGCTGGTATTATCCCGGTTGCTGa                                                                                    >  1:305332/1‑144 (MQ=255)
                                                                         tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAAGAAAAACCGCAAGGTGGGGCCGGATAAGGAAACAACCTTTGACCCTTTTTTGCTGTCAATATCCCGGGGGGTGGGCACATc                                                                             >  1:141953/1‑148 (MQ=255)
                                                                          ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGGTAAAAAGGCAACCTGTGACCATTTTTCGCTGACATCATCCTTGGGGCTGAcccc                                                                                >  1:624432/1‑142 (MQ=255)
                                                                            ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGATAAAGAGGCAAACTGTGACCATTTTTTCCTGACATCATCCTGGTGGCGGGCCACCTCACtta                                                                        >  1:171347/1‑148 (MQ=255)
                                                                              tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCTGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaa                                                                       <  2:253538/149‑1 (MQ=255)
                                                                               gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaaa                                                                      <  2:935844/149‑1 (MQ=255)
                                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGGGAAGTAGTCAACCTGTGACCATTCTTCGCCTACATCATCCTCGTGGGTGACCACATCACCGaaa                                                                      >  1:562725/1‑148 (MQ=21)
                                                                                     ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTg                                                               <  2:436110/150‑1 (MQ=255)
                                                                                       aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACCGCCAAGGGGGACCGGGAAAGTAGTCCACCTGTGGCCATTCCTCCCTGACCTCATCCTGGTGGGTGACCCCATAACTCAAACGCGGCgggg                                                            >  2:191775/1‑151 (MQ=255)
                                                                                        cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                           <  2:268869/151‑1 (MQ=255)
                                                                                        cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                           <  2:851356/151‑1 (MQ=255)
                                                                                              agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACGGCAAGGGGGGGCCCGGTAAGAAGTCCAACTGTGACCATTCTTCGCTGACATCCCCCTGGCGGCAGACCACATCACTTAAAAGCGgggggggg                                                          >  1:351168/1‑144 (MQ=12)
                                                                                               gtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTcc                                                    <  2:950856/151‑1 (MQ=255)
                                                                                                  gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAGGGGGGACCGGATAAGGAGTCAACCCGTGACCATTCTTCGCTGACATCATCCTGGTGGGTGACCACATCACTGAAACGCGGGGGGGGGAGGTCCttt                                                 >  1:581822/1‑150 (MQ=18)
                                                                                                  gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGTGGGGCCGGGTAAAGAGTCAAACCGTGACCATTCTTCGCTGACATCATCCTGGGGGCGGACCACATCACTGAAACGCGGGGGGGGGGGGGCCt                                                   >  2:172290/1‑149 (MQ=18)
                                                                                                    aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGTAAAAACCGCAAGGTGGGGCCCGCTAAGAAGTCAACCTTTGACCCTTCTTCTCTTACATCATCCCCGTGGGTGGCCACCACACTTAAAAGCGGTGGGGGGGGGGCCttta                                                >  2:588874/1‑148 (MQ=255)
                                                                                                     aCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAAGGGGGACCGGATAAATAATCAACCCGTGAACACTTTTTCCTGACATCATCCTGGGGGGTGACCAAATCACTGAAACGCGGGGGGGGGAGGGCCTTaaa                                               >  1:787370/1‑150 (MQ=255)
                                                                                                       tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACTGCAAGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGGCCACATCACTGAAACGCGGGGGGGGGGGGGCCTTTAAttt                                            >  2:358009/1‑151 (MQ=18)
                                                                                                         ccATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc                                          <  1:231285/151‑1 (MQ=255)
                                                                                                         ccATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc                                          <  1:201791/151‑1 (MQ=255)
                                                                                                             ccTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                      <  2:406526/151‑1 (MQ=255)
                                                                                                             ccTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                      <  2:788675/151‑1 (MQ=255)
                                                                                                              cTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGTGGGGCCCGATAAGTAGTCAACCTGTGACCATTCTTCTCTGACATCATCCTGGGGGCTGACCACATCACTGAAAAGCGGTGGGGGGAGGGCCTTTAATTTACTGAt                                      >  2:622324/1‑150 (MQ=255)
                                                                                                                   aGCCATCGGGCGACTGGCGGGGCAGATAAAACCGCAAAGGGGGGCCGGATAAAAAGTTAACCCGTGTCCACTCTTCGCTGACATCCTTCTGGTGGGTGACCACCTCAACGAAACACGGAGGGGGGAGGGGCTTAAAATTAATGATGGCCtt                                >  2:961345/1‑151 (MQ=255)
                                                                                                                   aGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                                <  2:814675/151‑1 (MQ=255)
                                                                                                                         cGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtc                          <  2:722556/151‑1 (MQ=255)
                                                                                                                                  ggCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                  <  2:622846/150‑1 (MQ=255)
                                                                                                                                  ggCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                 <  1:603455/151‑1 (MQ=255)
                                                                                                                                    cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt               >  2:198156/1‑151 (MQ=255)
                                                                                                                                       ggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa             <  2:431756/151‑1 (MQ=255)
                                                                                                                                         tACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc           >  1:250643/1‑151 (MQ=255)
                                                                                                                                                   ccTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt  >  2:826254/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT  >  NC_000913/266204‑266499

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