Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC61.8% 23.6 / 43.7 34T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (13/0);  new base C (0/21);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.08e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-07

TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376358‑1376655
                                                                                                                                                       |                                                                                                                                                      
ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                     >  2:348919/1‑151 (MQ=255)
                             ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                         >  2:709272/1‑150 (MQ=255)
                             ggACGGCGACCGAAAAAGGGTGAGTGATCCCCGGCGCGTACTGAAACGGGTCGCCGCGCGTACCGGACCCGTGGTTTATCTCGACCCGTAAACCCCCTATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                          <  2:224798/149‑1 (MQ=255)
                                   cgcccaaaaaaTGGTGAGTGTTCCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCCGGGGTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATGATCACCCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                  <  1:624854/144‑1 (MQ=255)
                                          aataTGGGGAGTGATCCCCTGCGCGTCATTGAACAGGTCGCCGCACGTACGGAACCCGTGGTGTCTCTCTAACCGTTAACCCCTGACCGGGCACATATGAACACCCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                            <  2:357152/147‑1 (MQ=255)
                                              tGGTGAGTGATCCCCTGCGCGTACTGAAACAGGCCCCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCTTAAACCCGCTATCAGGATCATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                        <  1:78842/150‑1 (MQ=255)
                                               ggTGAGTGTTTCCCGGCGGGTCCGGGAAAAGGTCGCGGCGCGTCCCGGACCGGGGGTGTATCTCAACCGTTTAACCCGCTATAAGGATCAAATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  1:459706/148‑1 (MQ=255)
                                               ggTGAGTGTTCCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTTTATCTCGAACGGTTAACCCGCTATCAGAATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                      <  1:544204/151‑1 (MQ=255)
                                                          cccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                           >  1:16848/1‑151 (MQ=255)
                                                             tGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  1:293807/1‑150 (MQ=255)
                                                                    aCTGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGACCCGTTAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   <  1:124525/149‑1 (MQ=255)
                                                                        gAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              >  1:629251/1‑150 (MQ=255)
                                                                         acaagttCGCCGCGCGTACCGGGCCGGTGTTTTATCTGGAACCGTAAACCCGCTAACAGGATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                             <  2:701050/144‑1 (MQ=255)
                                                                         aaagggTCGCCGCGCGTCCCGGAACCGTGTTGTACCCGGACCCTTAAACCCGCTATGAGGATCATATGATCACCCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            <  1:36554/147‑1 (MQ=255)
                                                                         aaCAGGTCGCCGCGCGTACCGGAACTGGGGTTTATCTCGAACGGTAAACCCGCTATCAGTATCATATGACCACCCCCCCCCCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              <  2:374052/149‑1 (MQ=255)
                                                                             ggTCGCCGCGCGTCCCGGAACGGTGTTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          <  1:557145/149‑1 (MQ=255)
                                                                                    gcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTATATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  1:319051/1‑149 (MQ=255)
                                                                                         tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  1:589186/1‑147 (MQ=255)
                                                                                           ccGGAACCGGTGGGTTATCCCAACCCTTAAACCCGTTAACAGAATCAAATAACCACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  2:377021/150‑1 (MQ=255)
                                                                                           ccGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           >  1:388900/1‑150 (MQ=255)
                                                                                            cGGAACCGTGGGTTTTTCTAGACCCGTTAAACCGCTACCAGAAGAATATAATCACCCCCCCCGCGTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  2:57095/149‑1 (MQ=255)
                                                                                                       ggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  >  1:375204/1‑149 (MQ=255)
                                                                                                             cactccacccttaaaaCCGCTATGAGCACCATGGTAACAACCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  1:192887/135‑1 (MQ=255)
                                                                                                              ctcGACCCGTAAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                               <  2:320944/145‑1 (MQ=255)
                                                                                                                      cttaaaaCCCCTATCAGAATCATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                   <  1:504390/146‑1 (MQ=255)
                                                                                                                        ttAACCCGCTATAAGAATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                                 <  2:216548/150‑1 (MQ=255)
                                                                                                                        ttAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                                >  1:62241/1‑151 (MQ=255)
                                                                                                                          aaacccGCTATCAGGATCATATAACCAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                               <  1:451119/147‑1 (MQ=255)
                                                                                                                           aaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             >  2:321411/1‑151 (MQ=255)
                                                                                                                            acccGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             <  1:484031/149‑1 (MQ=255)
                                                                                                                              ccGCTAACAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                                <  2:32567/146‑1 (MQ=255)
                                                                                                                                   aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                       >  2:713319/1‑150 (MQ=255)
                                                                                                                                            aaatGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt               <  2:188835/147‑1 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  2:347999/1‑151 (MQ=255)
                                                                                                                                                       acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:520706/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376358‑1376655

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