Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913574,4861.C20.8% 60.3 / 19.8 24intergenic (+147/+43)quuD/nmpCDLP12 prophage; putative antitermination protein/DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base . (9/10);  new base C (2/3);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCT  >  NC_000913/574350‑574621
                                                                                                                                         |                                                                                                                                       
gTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGCGGGAGGcggt                                                                                                                               <  7:18675/148‑4 (MQ=255)
                                  aCGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCttttt                                                                                                         <  4:3752/135‑1 (MQ=255)
                                  aCGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCttttt                                                                                                         >  3:3752/1‑135 (MQ=255)
                                          aaaaaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGt                                                                                    >  3:76541/1‑148 (MQ=255)
                                                 aaGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGa                                                                            <  6:64184/149‑1 (MQ=255)
                                                  aGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTa                                                                                   <  6:136998/141‑1 (MQ=255)
                                                  aGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTa                                                                                   >  5:136998/1‑141 (MQ=255)
                                                               aCAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGAt                                                              <  4:146996/149‑1 (MQ=255)
                                                                   aaGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCttttt                                                                                                         <  2:158056/102‑1 (MQ=255)
                                                                   aaGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCttttt                                                                                                         >  1:158056/1‑102 (MQ=255)
                                                                            ataGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTg‑g                                                                                                                                        <  3:119994/62‑1 (MQ=255)
                                                                            ataGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTg‑g                                                                                                                                        >  4:119994/1‑62 (MQ=255)
                                                                            ataGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCt                                                  >  4:12975/1‑148 (MQ=255)
                                                                                ttAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTacac                                             >  1:202127/1‑149 (MQ=255)
                                                                                   gAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGa                                          <  2:202127/149‑1 (MQ=255)
                                                                                     aGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGt                                        <  3:12975/149‑1 (MQ=255)
                                                                                                  ttCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTg‑g                                                                                                                                        >  5:63709/1‑40 (MQ=255)
                                                                                                  ttCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTg‑g                                                                                                                                        <  6:63709/40‑1 (MQ=255)
                                                                                                         ggAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAAt                    <  4:76541/149‑1 (MQ=255)
                                                                                                                gCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGa                                          >  2:181689/1‑120 (MQ=255)
                                                                                                                gCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGa                                          <  1:181689/120‑1 (MQ=255)
                                                                                                                 cTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTa                                                <  5:183571/113‑1 (MQ=255)
                                                                                                                 cTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTa                                                >  6:183571/1‑113 (MQ=255)
                                                                                                                               agatatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATtcttct  >  4:122006/3‑145 (MQ=38)
                                                                                                                               agatatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATtcttct  <  3:122006/143‑1 (MQ=38)
                                                                                                                                 atatGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAAt                    >  6:26118/1‑125 (MQ=255)
                                                                                                                                  tatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgg              >  3:83498/1‑130 (MQ=38)
                                                                                                                                  tatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgg              <  4:83498/130‑1 (MQ=38)
                                                                                                                                         |                                                                                                                                       
GTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCT  >  NC_000913/574350‑574621

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