Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,0230AT72.1% 54.4 / 59.4 68G123G (GGT→GGAnohANohA
Reads supporting (aligned to +/- strand):  ref base A (12/7);  new base T (22/27);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-03
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTG  >  NC_000913/1633933‑1634115
                                                                                          |                                                                                            
ccttttCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCAcc                                                                                            <  1:1295662/90‑1 (MQ=1)
               aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCACGAATGCTGGCTACTTCTCCGGCCCTCCGCGACAg                                                                             >  1:5013349/4‑93 (MQ=0)
               aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                             >  2:7136342/4‑93 (MQ=0)
               aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                             >  2:1079211/4‑93 (MQ=0)
               aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                             >  1:3908292/4‑93 (MQ=0)
               aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                             >  1:3139736/4‑93 (MQ=0)
               aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                             >  2:4081035/4‑93 (MQ=0)
               aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                             >  1:2214147/4‑93 (MQ=0)
                               agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                             <  1:3632982/90‑1 (MQ=0)
                               agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                             >  2:4812722/4‑93 (MQ=0)
                               agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                             >  1:8094508/4‑93 (MQ=0)
                               agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                             >  2:8323188/4‑93 (MQ=0)
                               agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                             <  1:5577826/90‑1 (MQ=0)
                               agctCCGGGAAGCCCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTATTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                             <  2:2922592/90‑1 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     >  2:3838494/5‑93 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     <  2:7734907/89‑1 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     >  2:1554131/5‑93 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     <  2:5708663/89‑1 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     >  2:8946125/5‑93 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     <  1:5021187/89‑1 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     >  1:1171938/5‑93 (MQ=0)
                                       gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                     <  1:1803465/89‑1 (MQ=0)
                                        aagcgccgcTGAACCGACATAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                    <  1:7750695/90‑1 (MQ=0)
                                        aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                    <  2:4714319/90‑1 (MQ=0)
                                        aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                    <  2:463259/90‑1 (MQ=0)
                                        aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                    >  2:4317570/4‑93 (MQ=0)
                                        aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                    <  2:2876776/90‑1 (MQ=0)
                                                 tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa                                           <  1:945715/90‑1 (MQ=0)
                                                                     atcaagaATTCTGGCTACTACTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAttttt                       <  2:2636927/89‑1 (MQ=0)
                                                                     atcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAttttt                       <  2:7320715/89‑1 (MQ=0)
                                                                     atcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAttttt                       <  2:2058666/89‑1 (MQ=0)
                                                                      tctagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCTGAATGCGGTGTCCACCACCTCAGCGTATTCTTTTGCAtttttt                      >  2:8965858/4‑92 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGCCAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGTATCTTTTGCAtttttt                      >  1:4014866/4‑92 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                      >  2:7517864/4‑92 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                      <  2:5583767/90‑2 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                      >  2:4751983/4‑92 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                      <  1:2659324/90‑2 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                      <  1:2659325/90‑2 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                      >  1:8561269/4‑92 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                      <  2:9878503/90‑2 (MQ=0)
                                                                      tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCCCAGCGGGATCTTTTGCAtttttt                      >  2:7225435/4‑92 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                >  1:8722093/4‑93 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  2:2177975/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  1:2804173/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  2:592915/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  1:9522514/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  1:3960162/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  2:97096/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  1:647544/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  1:1375327/90‑1 (MQ=0)
                                                                            atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc                <  2:813397/90‑1 (MQ=0)
                                                                                  gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg          <  1:3825496/93‑1 (MQ=31)
                                                                                   cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt         >  2:745874/1‑93 (MQ=32)
                                                                                    caTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc        <  1:4799219/92‑1 (MQ=18)
                                                                                    aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc        >  2:9711953/1‑93 (MQ=33)
                                                                                    aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc        <  2:2158454/93‑1 (MQ=33)
                                                                                     aTTTCACGTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGTGATTTTTCAGCTCCTGTGCGTCg       >  2:6168764/1‑93 (MQ=255)
                                                                                     aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg       >  2:1813895/1‑93 (MQ=33)
                                                                                     aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg       <  1:8091176/93‑1 (MQ=33)
                                                                                     aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg       >  1:4485965/1‑93 (MQ=33)
                                                                                     aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg       >  1:3838734/1‑93 (MQ=33)
                                                                                      tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGAGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg      >  1:1019359/1‑93 (MQ=18)
                                                                                       ttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGc     >  2:4238577/1‑93 (MQ=33)
                                                                                       ttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGc     >  2:1283783/1‑93 (MQ=33)
                                                                                        tCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGcc    >  2:5103996/1‑93 (MQ=33)
                                                                                        tCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGcc    <  1:4619439/93‑1 (MQ=34)
                                                                                         cACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCt   >  2:5840472/1‑93 (MQ=33)
                                                                                         cACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCt   >  1:4409395/1‑93 (MQ=33)
                                                                                         cACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCt   <  2:9896842/93‑1 (MQ=33)
                                                                                          aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg  >  2:3806781/1‑93 (MQ=33)
                                                                                          aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg  >  1:7898490/1‑93 (MQ=33)
                                                                                          |                                                                                            
CCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTG  >  NC_000913/1633933‑1634115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: