Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913589,9340TG51.3% 42.2 / 64.6 39T341P (ACC→CCC) nfrAbacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (2/17);  new base G (20/0);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 3.35e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-06

ATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  NC_000913/589786‑590072
                                                                                                                                                    |                                                                                                                                            
aTCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTgg                                                                                                                                            <  1:153367/151‑1 (MQ=255)
         cAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggggggggcgg                                                                                                                                    >  1:967465/1‑142 (MQ=255)
                        aGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTGCGCt                                                                                                                           >  2:761478/1‑144 (MQ=255)
                          gCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGTGTGGGGGGGGacgaacctttgccc                                                                                                                            >  1:696412/1‑128 (MQ=255)
                             ggTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGa                                                                                                               <  2:323226/151‑1 (MQ=255)
                               tAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCCGGTTTTTCGCCGCGGCGTAccgtcc                                                                                                             >  1:564313/1‑146 (MQ=255)
                                aaGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACct                                                                                                            <  2:484316/151‑1 (MQ=255)
                                   aTTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACcttct                                                                                                         >  2:396486/1‑151 (MQ=255)
                                      tGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCttt                                                                                                       <  1:196919/150‑1 (MQ=255)
                                       gCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCCGTTTTTTCGCTGCGTCGTACTGACCTTCttttt                                                                                                     >  2:304603/1‑150 (MQ=255)
                                          ggTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCCGTTTTTTCGCGGCGGCGGTCTGAACTTCtttttaa                                                                                                   >  1:118665/1‑147 (MQ=255)
                                            ttCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGTGAGCAGGTTTTTGGCTGCGTCGTACTGGCCTTCtttt                                                                                                      >  2:838272/1‑145 (MQ=255)
                                              cTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCCAGCATTTCATTGGCGGGGAGGGGGGGGAGGAGGTTTTTCGCTGCGGCGTACTGGCCTTCTTTTaa                                                                                                    >  1:685634/1‑145 (MQ=255)
                                                    gATACAGGAATCGTGCGAGACGCAGAGCTGCAGGCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCGCGAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTGCTGACCTTCTTTGAACAGGGCCGGTAg                                                                                        <  1:27491/151‑1 (MQ=255)
                                                    gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTa                                                                                         <  2:628915/150‑1 (MQ=255)
                                                       aCAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGt                                                                                     <  1:107975/151‑1 (MQ=255)
                                                        cAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGGTACGGGTCGCCCCGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGt                                                                                     <  2:987746/150‑1 (MQ=255)
                                                           cAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTcgcg                                                                                 <  2:138810/151‑1 (MQ=255)
                                                           cAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGTGAGCAGTTTTTGCGCCGCGGCGTACTGGACTTTTTTTTAAAGCACCCGTAGCGAcgc                                                                                  >  2:435155/1‑150 (MQ=255)
                                                            aaTCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGCCGTTTTTTGGCTGCGTCGTACTGACCTTCTTTTAAAAGCACCGGGAGCGTcgcg                                                                                 >  2:68514/1‑150 (MQ=255)
                                                                  gCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTGCGCTGCGCCGTACTGACCTTTTTTTTACAAAACCCGGAGCGTCccgc                                                                                >  1:194087/1‑145 (MQ=255)
                                                                     aGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACAt                                                                       <  1:514635/151‑1 (MQ=255)
                                                                        cGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGCAGGTTTTGCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGGCGCGCCAACAACATa                                                                      >  2:926631/1‑149 (MQ=255)
                                                                         gCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAAAAGCACCGGGAGCGTCGCGCcaacaa                                                                          >  2:612562/1‑144 (MQ=255)
                                                                             agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑gg                                                              <  2:92913/151‑1 (MQ=255)
                                                                              gCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGTTTTTTCGCTGCGTCGGTCTGACCTTCTTTTTAAAGCACCGGGAGCGGCGCGCCAACAACAAACTGGC‑ggg                                                             >  1:650186/1‑150 (MQ=255)
                                                                                  cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑t                                                         <  2:364719/150‑1 (MQ=255)
                                                                                    gCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGAGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTAAAAGCACCGGTAGCGTCGCGCCCACAACATACTGGC‑GGgtg‑t                                                         >  2:758870/1‑148 (MQ=255)
                                                                                       ttGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGGGCAGTTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGGACCGGGAGCGGCGCGCCAAAAAAACACTGGC‑ggggtgtcggg                                                     >  1:503943/1‑141 (MQ=255)
                                                                                               ggTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGCGAGAGCCGGTTTTTCGCTCCCTCGTACTTTCCTTCTTTTTTCAGCACTCGGCGCGTCGGCCCAACCACCAACtggt‑ggggggccggcgac                                                  >  2:902514/1‑133 (MQ=255)
                                                                                                gTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGAGCAGGTTTTGCGCTGCGTCGTAATGACCTTCTTTTTACAGCACCGGGAGCGGCGCGCCAACAACAAACTGGG‑GGGTG‑TCggggcaaagtgac                                           >  2:53344/1‑139 (MQ=255)
                                                                                                   gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGtat                                       <  1:621058/151‑1 (MQ=255)
                                                                                                                  tAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGAGCAGGTTTTGCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGGCGCGCCAACAACATACTGGC‑GGGTGTTTGGGAAACTTTACCGGAAAATTTCCCAAcccc                         >  2:575096/1‑151 (MQ=255)
                                                                                                                      gCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAAc                    <  1:866390/151‑1 (MQ=255)
                                                                                                                         cctcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACggg                 <  1:701175/151‑1 (MQ=255)
                                                                                                                            cAAGCATTTCATTGGCGGGGGGGGGGGCGAGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCCACAACATACTGGCGGGTTT‑TGGGGAAACATTACCGCTAAATTTCCCAACCCCGGACCgggg                >  2:130713/1‑150 (MQ=255)
                                                                                                                             aaGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTa             <  1:733897/151‑1 (MQ=255)
                                                                                                                                     cATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTa     <  2:538782/151‑1 (MQ=255)
                                                                                                                                        tGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGGTT‑TTGGCAAAAAGTAACCGATAAATTGCCCAAGCCCGAAAGGGGGTTGGGCTgggttt  >  2:383994/1‑146 (MQ=255)
                                                                                                                                           cGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGGAAACTGTACCGTATAATTTGCCAACGCCCGAACGGGGGTAGCGCTGGAttt  >  2:356232/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
ATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  NC_000913/589786‑590072

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