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

CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAA  >  NC_000913/266214‑266505
                                                                                                                                                |                                                                                                                                                    
cctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGAc                                                                                                                                                 <  1:686750/150‑1 (MQ=255)
   ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCgg                                                                                                                                              <  1:584520/150‑1 (MQ=255)
     tATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGAgggggcgg                                                                                                                                              >  1:532282/1‑143 (MQ=255)
         ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAagta                                                                                                                                       <  1:569643/151‑1 (MQ=255)
                 aTGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACTGCAAGGGGGGGCCGGATAagtagggca                                                                                                                                  >  2:954264/1‑144 (MQ=255)
                   gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCtgt                                                                                                                             <  1:136843/151‑1 (MQ=255)
                   gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCtg                                                                                                                              <  2:471554/150‑1 (MQ=255)
                         ggCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCAt                                                                                                                      <  1:914818/151‑1 (MQ=255)
                           cAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGGGGGGCCGGATAAGAAGTCAACCTCTG‑aac                                                                                                                        >  1:378241/1‑145 (MQ=255)
                           cAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGAAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGCGGTAAAACTGCCAGGGGGGGCCGGAGAAACAGGCAACCTGGG‑GCCAtt                                                                                                                     >  2:314148/1‑150 (MQ=255)
                              tACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCCAGGGGGGGCCGGGTAAGAAGTCAACCCCTG‑ACCAtt                                                                                                                     >  2:329448/1‑147 (MQ=255)
                               acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAAGGGGGGCCGGGTAAATAGTCAACCTGTG‑ACCcttctt                                                                                                                  >  1:112539/1‑149 (MQ=255)
                                    gcgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGAc                                                                                                           <  1:443484/151‑1 (MQ=255)
                                      gcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGAcat                                                                                                         <  2:767291/151‑1 (MQ=255)
                                      gcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGGCCGGATAAATAGTCAACCTGTG‑ACCCATTTTTGCTGAc                                                                                                           >  1:618002/1‑149 (MQ=255)
                                       cATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGAcatc                                                                                                        <  2:133367/151‑1 (MQ=255)
                                         ttCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGAcatca                                                                                                       <  1:174475/150‑1 (MQ=255)
                                           catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACatcatc                                                                                                     <  1:557778/150‑1 (MQ=255)
                                              caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACCGAAAACCTTCAAGGGGGGGCCGGCAAAGTAGGCAACCCGTGAACCATTCTTC‑CTGAcaac                                                                                                        >  1:321669/1‑142 (MQ=255)
                                                   gCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGGTAAGTAGTCAAACTGTG‑ACCATTCTTCTCTGACATCATCCtggtgg                                                                                              >  1:83936/1‑149 (MQ=255)
                                                              ttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACCGCCAGGGGGGGCCGGGTAAATAGTCAACCTGTG‑ACCATTCTTTGCCGACATCATCCTGGTGGCTGAcc                                                                                        >  2:920551/1‑144 (MQ=255)
                                                                ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACAGCCGGGTGGGGCGGGGTTCGGAGGCATACTGTG‑GCCATTCTTCGCTGACATCCTCCGCGGCGCTGCCCACATc                                                                                   >  1:453728/1‑147 (MQ=255)
                                                                  ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGGGGGGCCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGGTGACCACATCAc                                                                                 >  1:828769/1‑147 (MQ=255)
                                                                     gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAc                                                                           <  2:593351/150‑1 (MQ=255)
                                                                       atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGAGTAAAAACCGCAAGGGGGGGCCGGGTAAAGAATCAAACCTTG‑ACCATTTTTTGCTGACATAATCCTGGGGGGTGACCCCAACAATGAAAagcg                                                                        >  1:355867/1‑151 (MQ=255)
                                                                       atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGGCCGGATAAATAGGCAACCTGTG‑ACCAATCTTCGCTGACAACATCCCGGGGGCTGACCACAACACTTaaa                                                                            >  1:300585/1‑147 (MQ=255)
                                                                          tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGACCGGAGAAATAGACCACCCGTG‑TTCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATAACCGAAACGCggcg                                                                     >  1:254824/1‑149 (MQ=21)
                                                                             aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTgggg                                                                  <  1:29313/151‑1 (MQ=255)
                                                                               gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTg                                                                <  1:794962/151‑1 (MQ=255)
                                                                                aGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                               <  2:459412/151‑1 (MQ=255)
                                                                                    agtggtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGt                                                            <  1:502220/150‑1 (MQ=255)
                                                                                             taccATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGACCGGGTAAAAAGTCAACCTGTG‑ACCATTTTTTGCTGACATCATCCTGGGGGCGGACCACATCACTAGAAAGCGGGGGGGGGGGGGCCTTAAAttt                                                  >  2:174689/3‑151 (MQ=255)
                                                                                             tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACCTAAAAACCGCCAAGTGGGGCCCGATAAATAAGCAACCCGTG‑ACCATTCTTTGCTGACATCCTCCCGGTGGCTGACCACATCACTGAAAAGCGGGGGGGGGAGGGCCTTTAAttt                                                  >  2:98929/1‑151 (MQ=255)
                                                                                              gCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGGTAAATAGTCAACCTGTG‑ACCCTTTTTCGCTGACATCATCCTGGGGGCTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCTTAAAatt                                                  >  1:444926/1‑150 (MQ=255)
                                                                                                 aTCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAACCTGCAAAGGGGGGCCGGGTAAAGAGTGAACCTGTG‑ACCCTTCTTTGCTTACATCATCCTTGTGGCTGGCCACATCACTGAAAAGCGGGGGGGGGGGGGCCTTTAAttttc                                                >  1:544745/1‑147 (MQ=255)
                                                                                                        cAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCt                                       <  2:104720/151‑1 (MQ=255)
                                                                                                             aTCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAgg                                  <  1:811799/151‑1 (MQ=255)
                                                                                                                           gggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTc                    >  2:376455/1‑151 (MQ=255)
                                                                                                                           gggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTc                    <  2:671893/151‑1 (MQ=255)
                                                                                                                                             aGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaaa  <  1:605889/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTG‑ACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAA  >  NC_000913/266214‑266505

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