Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913287,1720GC45.5% 32.6 / 15.0 33P128P (CCG→CCC) yagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base G (3/15);  new base C (15/0);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GCTGCTCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTG  >  NC_000913/287043‑287316
                                                                                                                                  |                                                                                                                                                
gctgctCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGACCC‑CTCCCCGTTCCCCGACGACGAt                                                                                                                             >  1:69198/1‑151 (MQ=255)
 ctgctCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCacgacgac                                                                                                                                <  1:83631/147‑1 (MQ=255)
 ctgctCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTTACCC‑CTCCCTGTTCCac                                                                                                                                      >  2:399084/1‑141 (MQ=255)
         cacCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGAACC‑CTCCCCGTTCCACGACGACGATGGcc                                                                                                                         >  2:325896/1‑146 (MQ=255)
            cGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCCGTTCCACGACGACGATGGCCGCaaa                                                                                                                    >  2:27335/1‑148 (MQ=255)
            cGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGGCCC‑CCCCCCCTTCCCCCACGCCGACGGCCGCAAATa                                                                                                                  >  1:412348/1‑150 (MQ=255)
              aCGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACta                                                                                                               <  1:462695/151‑1 (MQ=255)
                 tGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGACCCGGCCCCTGTT‑CACCACGACGATGGCCGCACATaacat                                                                                                              >  2:824633/1‑145 (MQ=255)
                  gAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGACCC‑CTCCCCGTTTCACGACGACGAAGGCCGCAAATACtata                                                                                                             >  2:584939/1‑149 (MQ=255)
                     gATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTAc                                                                                                        <  2:488182/151‑1 (MQ=255)
                      aTTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGATCC‑CTCCCTGGTTCTCGACGACGATGGGCGCCAATACTATATCTAc                                                                                                        >  2:266328/1‑150 (MQ=255)
                         gTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGc                                                                                                     <  1:379660/150‑1 (MQ=255)
                               tCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGACCC‑CTCCCCGTTCCACCACGACGATGGCCGCAAAAACCATATCTACCCCCCCTggg                                                                                              >  1:237154/1‑151 (MQ=255)
                                    gTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCcg                                                                                         <  2:327229/151‑1 (MQ=255)
                                         aaaaCGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGACCC‑CGCCCTGTTCCACCACGACGATGGCCCCAAATACTATAACTACCCCCCGCGGGGGGCGCGcc                                                                                     >  2:688459/1‑150 (MQ=255)
                                             cGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCcaccac                                                                                <  2:845989/151‑1 (MQ=255)
                                                        tCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTTACCC‑CTCCCCGTTCCCCGACGACGATGGCCCCCAATACTCTATCTACCGCCCGCGGGGGCCGCGCCCCCCCACCAAccc                                                                        >  1:528600/1‑148 (MQ=255)
                                                            cGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGcaca                                                                   <  1:725425/149‑1 (MQ=255)
                                                            cGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCacaa                                                                  <  1:198650/150‑1 (MQ=255)
                                                                 ccGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGACCC‑CTCCCCGTTCCCCGCCGCCGAAGGCCGCAAAAACTATATCCACCGCCCCTGGGGGGCGCGCCCCCCCCCCAACCCCCCCCACAcc                                                              >  1:410711/1‑149 (MQ=255)
                                                                             tCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTGACCC‑CCCCCCGTTCCACCACGAAGATGGCCCCCCAAACTATATATACCCCCCCCGGGGGGGCCGCCACCACACCAACCCCCACAACAcccacgcggttc                                                    >  1:792878/1‑137 (MQ=255)
                                                                                         ggAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCg                                      <  1:766235/149‑1 (MQ=255)
                                                                                          gagcgagcCAATCCCGATGGGCAACGGCGGGTTTTACCC‑CTCCCTGTTCCCCGACGACGAAGGCCGCCAATACTATATCTACCCCCCGGGGGGGGCGCGCCACCACAGCAAACCCCACAACACCATCGTGTTTCAGGCGGTTGGCCCGc                                     >  2:783511/1‑149 (MQ=255)
                                                                                          gagcgagcCAATCCCGATGGGCAACGGCGGGTTTGACCC‑CTCCCCGTTCCACGACGACGATGGCCGCAAATACTATATCTACCCCCCCCGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTTAccc                                       >  2:53348/1‑147 (MQ=255)
                                                                                              gagcCAATCCCGATGGGCAACGGCGGGGTTTTCCC‑CTCCCCGTTCCCCCACGACGAAGGCCGCCAAAACTATATCTACCCCCCCCGGGGGCCGCGCCCCCACCGCAACCCCCACCACACCACCGTGTTACAGGCGGTTGACCCCCCGAcc                                >  1:456734/1‑150 (MQ=255)
                                                                                                      cccGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGct                        <  2:842881/150‑1 (MQ=255)
                                                                                                       ccGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGctct                      <  2:204757/151‑1 (MQ=255)
                                                                                                        cGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCCCGACGACGATGGCCGCAAATACTATATCTACCCCCCGGGGGGGCCGCGCCACCACAGCAACCCCCACAACACCATCGTGTTTCAGGCGGTTGACCCCCCGACCGGCCCGCtctc                     >  2:177722/1‑151 (MQ=255)
                                                                                                        cGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCtctc                     <  1:622818/151‑1 (MQ=255)
                                                                                                            gggCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGcc                  <  2:801890/150‑1 (MQ=255)
                                                                                                                  cggcggGTTTGACCC‑CTCCCTGTTTCACCACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCCCCACAGCCACCCGCCCCAAACCCTCGCGTTAAAGGCGGTTTAACCCCCGACCCGCCCGCCCTCcccc                 >  1:379543/1‑145 (MQ=255)
                                                                                                                          ttGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGc    <  1:252267/150‑1 (MQ=255)
                                                                                                                            gACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTg  <  1:584901/150‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                                
GCTGCTCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGTTTGACCC‑GTCCCTGTTCCACGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTG  >  NC_000913/287043‑287316

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