Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,965,1010TG30.0% 65.6 / 10.5 40N19H (AAC→CAC) thyAthymidylate synthase
Reads supporting (aligned to +/- strand):  ref base T (9/19);  new base G (12/0);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 6.16e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGTCGCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTG  >  NC_000913/2964956‑2965248
                                                                                                                                                 |                                                                                                                                                   
gtgtCGCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGttttt                                                                                                                                                  <  1:1389179/149‑1 (MQ=255)
     gCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCtgtg                                                                                                                                             <  2:1130341/149‑1 (MQ=255)
     gCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCtgtg                                                                                                                                             <  2:702188/149‑1 (MQ=255)
             gAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGGCttag                                                                                                                                     >  2:1730746/1‑147 (MQ=255)
                    cagcagTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCa                                                                                                                              <  1:1085152/149‑1 (MQ=255)
                    cagcagTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCa                                                                                                                              <  1:1086160/149‑1 (MQ=255)
                        agTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCACAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCAATCGTATAGCaaatt                                                                                                                          >  2:1107966/1‑145 (MQ=255)
                           tCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTg                                                                                                                       <  2:293494/149‑1 (MQ=255)
                             aTGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCa                                                                                                                     <  1:744715/149‑1 (MQ=255)
                                 atgatgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAg                                                                                                                 >  2:692350/1‑149 (MQ=255)
                                   gatgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGtt                                                                                                               <  1:1032142/149‑1 (MQ=255)
                                    atgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGACAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTc                                                                                                              >  2:1685952/1‑149 (MQ=255)
                                       gAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGGCGTCGTGGGGGACTTTTTTGATCAGTTGTaa                                                                                                           >  1:892219/1‑149 (MQ=255)
                                         aCGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTTTTCCTTCGTCGTGCACTTTTTTCATAAGTTGTAAAt                                                                                                         >  1:2100482/1‑149 (MQ=255)
                                                   cAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTTTGTCTTGGTCGAGCACATTTTTAATAAATTATAAATACTGTTTCaa                                                                                               >  1:1276057/1‑148 (MQ=255)
                                                   cAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGGCATAGTCGAGCACGTTTTTCATAAGGTGTAAATAATTTTTAAt                                                                                               >  2:700557/1‑149 (MQ=255)
                                                   cAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTATGTTTTCCTTCGTGGAGCACGTTTTGAATAAGTTATAAATAATTTTTAAt                                                                                               >  2:699991/1‑149 (MQ=255)
                                                     aCGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATTTGTGGCATACTCGACCAAATTTTGGATCAGGTATAAATACTGGTTAAttg                                                                                             >  2:1017186/1‑147 (MQ=255)
                                                              tGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAgg                                                                                    <  1:1426583/149‑1 (MQ=255)
                                                                 cACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATTGTTCCTCAGGaca                                                                                 >  2:2061962/1‑147 (MQ=255)
                                                                 cACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATGTGTGCCATCGGCGAGCACTTTTTGCATCAGTTCTAAATAAAGTTTAAAGGTTACTCAGGaaa                                                                                 >  2:1006192/1‑149 (MQ=255)
                                                                     aGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACgtg                                                                             <  1:537667/149‑1 (MQ=255)
                                                                        gggAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTTGCTTAGTCGACCACGTTTTTGAAAAGTTATAAATAAATTTTAATGGTTGCTCAAGAAACATGTTg                                                                          >  2:240335/1‑149 (MQ=255)
                                                                             tCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCtgtg                                                                     <  1:1274478/149‑1 (MQ=255)
                                                                                       aGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTTTGCCTTCGTGGAGCCCGTTTTTCATAAGTTCTAAATACTGTTTGAAGTTTACTCAGGTAACGTGTTTCGGTGGGGTggggcg                                                           >  1:1067811/1‑144 (MQ=255)
                                                                                        ggTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTTTGTGTGGCGTCGTCGGGCACATTTTTGATAAGTTGTAAATAATTTTTAAATGTTGCTAAGGGAAACTGTTTGTGTGGGGGGcggcga                                                          >  1:763935/1‑149 (MQ=255)
                                                                                        ggTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTGATGGTTGCTAAGGAAACGTGTTGCGGTGGGCTGcgacga                                                          >  2:523996/1‑149 (MQ=255)
                                                                                          ttAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATTTGTGCCTTAGTCGAGCACTTTTTTGATCAGTTCTAAATAATGTTTGATTGTTACTAAGGAAACGTGTTGCTGTGGGGTGCGACGata                                                        >  2:96299/1‑149 (MQ=255)
                                                                                          ttAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGCCGTAGTGGTGCACAGTTTTGAACAGGTCTTAAAAATGTTTGATTGTTGCTGAGGAAACGTGTTTGTGTGGGGGGCGGCGata                                                        >  2:785203/1‑149 (MQ=255)
                                                                                           tAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATAAATTTTGAAGGTTCCTCAGGAAACGTGTTGCGGTGGGCGGCGACGatat                                                       >  2:779426/1‑149 (MQ=255)
                                                                                               aCGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGccc                                                   <  2:993512/149‑1 (MQ=255)
                                                                                                       gatgaCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACcatc                                           <  2:538300/149‑1 (MQ=255)
                                                                                                         tgaCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACcatcat                                         <  1:1170552/149‑1 (MQ=255)
                                                                                                            ccAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGAt                                      <  1:699991/149‑1 (MQ=255)
                                                                                                            ccAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGAt                                      <  1:700557/149‑1 (MQ=255)
                                                                                                             cAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATc                                     <  1:565139/149‑1 (MQ=255)
                                                                                                                aaaTGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATcaca                                  <  1:303801/149‑1 (MQ=255)
                                                                                                                                ggttccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCg                  <  2:2100482/149‑1 (MQ=255)
                                                                                                                                   tcaggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGa               <  1:268779/146‑1 (MQ=255)
                                                                                                                                              cGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAAtt    >  1:1947129/1‑149 (MQ=255)
                                                                                                                                              cGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAAtt    >  2:1852933/1‑149 (MQ=255)
                                                                                                                                                tttttCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTg  <  2:1837559/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTGTCGCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTG  >  NC_000913/2964956‑2965248

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