Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,717,671	0	A	C	63.6%	35.5 / 47.8	44	T286T (ACT→ACG)	yfiF	putative methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (5/11);  new base C (15/13);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

AAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGAT  >  NC_000913/2717525‑2717811                                                                                                                                                   |                                                                                                                                             aaaaCGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGttt                                                                                                                                          <  2:28365/151‑1 (MQ=255)             gCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGtt                                                                                                                              <  1:453726/151‑1 (MQ=255)             gCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGt                                                                                                                               <  1:645097/150‑1 (MQ=255)               gACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTa                                                                                                                            <  2:206353/151‑1 (MQ=255)                 cAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTAcc                                                                                                                          >  1:228119/1‑151 (MQ=255)                  agagATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTAccc                                                                                                                         >  1:188791/1‑151 (MQ=255)                   gagaTATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCt                                                                                                                        <  2:571877/151‑1 (MQ=255)                      atatTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCt                                                                                                                     >  2:257393/1‑151 (MQ=255)                          tCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGc                                                                                                                  <  2:911088/150‑1 (MQ=255)                                    aaCGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTc                                                                                                       <  2:292858/151‑1 (MQ=255)                                    aaCGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTc                                                                                                       <  1:405271/151‑1 (MQ=255)                                    aaCGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTc                                                                                                       <  1:76693/151‑1 (MQ=255)                                              gTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGgtg                                                                                              <  1:257834/150‑1 (MQ=255)                                                aCCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTaa                                                                                           <  2:603688/151‑1 (MQ=255)                                                                      tCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATcc                                                                      >  1:254271/1‑150 (MQ=255)                                                                       cGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCa                                                                     <  1:312530/150‑1 (MQ=255)                                                                       cGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCa                                                                     <  2:766662/150‑1 (MQ=255)                                                                                   gTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAAt                                                        >  2:302324/1‑151 (MQ=255)                                                                                     gCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAAtgt                                                      <  1:32522/151‑1 (MQ=255)                                                                                       ggCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAAtgttg                                                    <  1:785403/151‑1 (MQ=255)                                                                                          aTCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCg                                                 >  1:591321/1‑151 (MQ=255)                                                                                           tCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGc                                                >  1:825457/1‑151 (MQ=255)                                                                                            ccGGTAACCCTTCATATTCCTGACCCAGTACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGc                                                <  2:135935/150‑1 (MQ=255)                                                                                            ccGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGcc                                               <  2:800124/151‑1 (MQ=255)                                                                                              ggTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGcc                                               >  2:35418/1‑149 (MQ=255)                                                                                               gTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGt                                            >  1:919939/1‑151 (MQ=255)                                                                                                tAACCCTTCATATTCCTGACCCAGCACCAGTGCCATTTTCGCGGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAg                                             <  1:480800/149‑1 (MQ=255)                                                                                                   cccTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAAt                                         <  2:878499/150‑1 (MQ=255)                                                                                                    ccTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCg                                       >  1:854746/1‑151 (MQ=255)                                                                                                     cTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCgg                                      >  2:572670/1‑151 (MQ=255)                                                                                                       tCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGc                                     >  2:118173/1‑150 (MQ=255)                                                                                                       tCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCt                                    >  1:186114/1‑151 (MQ=255)                                                                                                              ccTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGt                              <  2:235210/150‑1 (MQ=255)                                                                                                               cTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGc                            >  2:82746/1‑151 (MQ=255)                                                                                                                      aGCACCAGTACCATTTTCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTgcgc                     <  1:704915/151‑1 (MQ=255)                                                                                                                      aGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTgcgc                     <  2:753135/151‑1 (MQ=255)                                                                                                                           cAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCt                >  1:452986/1‑151 (MQ=255)                                                                                                                           cAGTACCATTTTCGCTGGCAGACTAGTCTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGTTGTAACCTGCCTGACGGAAATCATCCAGCGCGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCTCCGCCt                <  1:369561/151‑1 (MQ=255)                                                                                                                                     ttCGCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTAc      >  1:766188/1‑151 (MQ=255)                                                                                                                                     ttCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTAc      >  1:14167/1‑151 (MQ=255)                                                                                                                                        gCTGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACgg    >  2:936567/1‑150 (MQ=255)                                                                                                                                        gCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGa   <  1:543863/151‑1 (MQ=255)                                                                                                                                         cTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGa   >  1:460612/1‑150 (MQ=255)                                                                                                                                          tGGCAGACTCGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGAt  >  2:538128/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                             AAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATCCGGTAACCCTTCATATTCCTGACCCAGCACCAGTACCATTTTCGCTGGCAGACTAGTTTTGAACAGCGGTTTACCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGAT  >  NC_000913/2717525‑2717811

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