Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,744,443	0	C	A	72.7%	43.6 / 20.4	33	V77F (GTT→TTT)	apt	adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (10/14);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCT  >  CP000730/1744307‑1744579                                                                                                                                         |                                                                                                                                         gCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGACATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGt                                                                                                                                      >  1:252407/1‑141 (MQ=255)    aataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCa                                                                                                                                   <  1:342066/141‑1 (MQ=255)     ataaATCATCTGTAATTAACACACGTTGCCCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTCCTAACAGGTGCaa                                                                                                                                  >  2:444583/1‑141 (MQ=255)     ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCaa                                                                                                                                  >  1:390856/1‑141 (MQ=255)          tcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCa                                                                                                                             <  1:325493/141‑1 (MQ=255)                gTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATcccc                                                                                                                       <  2:295395/141‑1 (MQ=255)                   aTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATTACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCAtt                                                                                                                    >  2:46899/1‑141 (MQ=255)                      aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGaa                                                                                                                 <  2:336388/141‑1 (MQ=255)                      aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGaa                                                                                                                 <  2:99803/141‑1 (MQ=255)                                gACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAg                                                                                                       <  2:71230/141‑1 (MQ=255)                                   cTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGAc                                                                                                    >  1:197238/1‑141 (MQ=255)                                            ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGa                                                                                           <  1:392218/141‑1 (MQ=255)                                                         cATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                                              >  1:82086/1‑141 (MQ=255)                                                              ttAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAg                                                                         <  2:239586/141‑1 (MQ=255)                                                                    cATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTa                                                                   >  1:31543/1‑141 (MQ=255)                                                                       ttGTACCATATTCTAGGTCATACTCATAACGAATTACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                                                <  1:46899/141‑1 (MQ=255)                                                                        tttACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAc                                                               <  1:225073/139‑1 (MQ=255)                                                                            ccATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGata                                                           >  1:165212/1‑141 (MQ=255)                                                                            ccATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGata                                                           <  2:205264/141‑1 (MQ=255)                                                                                atTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATATAGCCACGCGCTTCAGGTCCTACAACGATATCaa                                                       <  2:120683/141‑1 (MQ=255)                                                                                  tCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACa                                                     >  1:326210/1‑141 (MQ=255)                                                                                          aTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTCCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc                                             <  1:444583/141‑1 (MQ=255)                                                                                            aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtt                                           <  2:390856/141‑1 (MQ=255)                                                                                                       aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                                >  2:1275/1‑141 (MQ=255)                                                                                                         gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCTACATCTCTGTCTTTTGCGTATTCTAc                              >  2:122055/1‑141 (MQ=255)                                                                                                         gACTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc                              >  2:29160/1‑141 (MQ=255)                                                                                                             tCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt                          >  1:213348/1‑141 (MQ=255)                                                                                                                 gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCGCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTAt                      <  1:117245/141‑1 (MQ=255)                                                                                                                     tAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt                  <  2:197238/141‑1 (MQ=255)                                                                                                                       aTTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTg                <  2:31543/141‑1 (MQ=255)                                                                                                                                    tCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGc   <  2:382221/141‑1 (MQ=255)                                                                                                                                    tCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGc   <  2:378561/141‑1 (MQ=255)                                                                                                                                     cTAACAGGTGCAAAGCCAATCCCCATCGAATAAGCTACAGGACACCCAACGATAAAGCCCCGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCt  >  2:12314/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                                         GCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCT  >  CP000730/1744307‑1744579

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