Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913141,1270CT23.7% 75.5 / 18.7 38W33* (TGG→TGA) gcdglucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (12/17);  new base T (4/5);  total (16/22)
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.37e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAGAGACACCATTTCAATTCATTAAT  >  NC_000913/140990‑141252
                                                                                                                                         |                                                                                                                             
ccccAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGAccagcc                                                                                                                              >  1:3967/1‑139 (MQ=255)
            gTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGTCAGCCTCCGCCaa                                                                                                                  <  2:561137/139‑1 (MQ=255)
            gTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCaa                                                                                                                  <  1:249659/139‑1 (MQ=255)
                            cTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGccc                                                                                                  >  1:571425/1‑139 (MQ=255)
                               cGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCa                                                                                               >  2:222813/1‑139 (MQ=255)
                                  ataGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAg                                                                                            <  2:101703/139‑1 (MQ=255)
                                                 cgcgTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGc                                                                             <  2:571425/139‑1 (MQ=255)
                                                             gcCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCaaaa                                                                                  <  1:702332/122‑1 (MQ=255)
                                                             gcCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCaaaa                                                                                  >  2:702332/1‑122 (MQ=255)
                                                              ccACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCaa                                                                                                                  <  1:152117/89‑1 (MQ=255)
                                                              ccACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCaa                                                                                                                  >  2:152117/1‑89 (MQ=255)
                                                              ccACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACg                                                                <  1:222813/139‑1 (MQ=255)
                                                                      tCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAAtcg                                                        >  1:229483/1‑139 (MQ=255)
                                                                         aGGCGACGCCGAGCATCACAAGGCCAGAGATAGGGTAGTACCAGTAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAAtcgtcg                                                     <  2:631504/139‑1 (MQ=255)
                                                                              aCGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGc                                                <  2:229483/139‑1 (MQ=255)
                                                                                          acaAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTaa                                    <  2:3967/139‑1 (MQ=255)
                                                                                                   gCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCaa                                                                                                                  >  1:389906/1‑52 (MQ=255)
                                                                                                   gCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCaa                                                                                                                  <  2:389906/52‑1 (MQ=255)
                                                                                                    cGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAAtcgtcg                                                     <  2:214582/112‑1 (MQ=255)
                                                                                                    cGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAAtcgtcg                                                     >  1:214582/1‑112 (MQ=255)
                                                                                                        aGGGTAGTACCAAGAGCCGCCAATCGCAACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGTGCTGTTAGCGTAACGAGAAGTCGTCGCGAGCCTGTATTGTTAATTGCCATAAAgaga                      <  1:115104/139‑1 (MQ=255)
                                                                                                              gtaCCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGAt                                                                                                         >  2:684388/1‑50 (MQ=255)
                                                                                                              gtaCCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGAt                                                                                                         <  1:684388/50‑1 (MQ=255)
                                                                                                                  cAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAGAGACACCATTTCa            <  1:194911/139‑1 (MQ=255)
                                                                                                                        ccgccAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAAtt                                  >  1:350719/1‑111 (MQ=255)
                                                                                                                        ccgccAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAAtt                                  <  2:350719/111‑1 (MQ=255)
                                                                                                                         cgccAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAgaga                      >  1:75218/1‑122 (MQ=255)
                                                                                                                         cgccAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAgaga                      <  2:75218/122‑1 (MQ=255)
                                                                                                                         cgccAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAGAGACACCATTTCAATTCAtt     >  1:645798/1‑139 (MQ=255)
                                                                                                                         cgccAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAGAGACACCATTTCAATTCAtt     >  2:263454/1‑139 (MQ=255)
                                                                                                                           ccAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGTGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGc                                                >  2:91122/1‑94 (MQ=255)
                                                                                                                           ccAATCGCGACCAGTCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGTGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGc                                                <  1:91122/94‑1 (MQ=255)
                                                                                                                           ccAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCa                              <  2:621280/112‑1 (MQ=255)
                                                                                                                           ccAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCa                              >  1:621280/1‑112 (MQ=255)
                                                                                                                           ccAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAgaga                      <  2:242602/120‑1 (MQ=255)
                                                                                                                           ccAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAgaga                      >  1:242602/1‑120 (MQ=255)
                                                                                                                            cAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAGAGACACCATTTCAATTCATTAat  <  1:263454/139‑1 (MQ=255)
                                                                                                                                    aCCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGa                                                  >  2:487907/1‑83 (MQ=255)
                                                                                                                                    aCCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGa                                                  <  1:487907/83‑1 (MQ=255)
                                                                                                                                         |                                                                                                                             
CCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCACAGCATCCAGGCGACGCCGAGCATCACAAGGCCAGCGATAGGGTAGTACCAGGAGCCGCCAATCGCGACCAGCCAGCCTCCGCCAATGAGTAGATACAGCCCGCAAAGCGCTGCAAAAAGGGCTGTTAGCGTGACGAGTAATCGTCGCGAGCCTGTATTGTTAATTGCCATAAAGAGACACCATTTCAATTCATTAAT  >  NC_000913/140990‑141252

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