Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	650,809	0	T	G	25.4%	144.2 / 11.8	59	T350P (ACC→CCC)	citC	citrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base T (16/28);  new base G (15/0);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 8.47e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT  >  NC_000913/650677‑650947                                                                                                                                     |                                                                                                                                            ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTGtttttg                                                                                                                                      >  2:1198444/1‑138 (MQ=255)          aTTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGcc                                                                                                                             <  1:1528099/140‑1 (MQ=255)              gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGt                                                                                                                         >  1:299707/1‑140 (MQ=255)              gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGt                                                                                                                         <  1:1332329/140‑1 (MQ=255)              gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGt                                                                                                                         >  1:1938168/1‑140 (MQ=255)                 ggTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCt                                                                                                                      >  1:1937537/1‑140 (MQ=255)                   tCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTgg                                                                                                                    <  1:677022/140‑1 (MQ=255)                      atatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCgg                                                                                                                 <  1:721288/140‑1 (MQ=255)                        atCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGa                                                                                                               <  2:1645164/140‑1 (MQ=255)                        atCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGa                                                                                                               <  2:1337813/140‑1 (MQ=255)                           cTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTg                                                                                                            <  1:142166/140‑1 (MQ=255)                            tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGc                                                                                                           >  1:1728008/1‑140 (MQ=255)                            tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTgg                                                                                                           >  2:75652/1‑139 (MQ=255)                              cGTATCGAGTGGGGCGATGCGTATCATCACGTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTc                                                                                                         >  1:86528/1‑140 (MQ=255)                              cGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTc                                                                                                         >  2:1628128/1‑140 (MQ=255)                                 aTCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGTCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTccaa                                                                                                      >  1:1145411/1‑137 (MQ=255)                                 aTCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAg                                                                                                     >  1:1101192/1‑140 (MQ=255)                                  tCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGc                                                                                                    <  1:1216974/140‑1 (MQ=255)                                    gAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAg                                                                                                  <  2:292940/140‑1 (MQ=255)                                     aGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGCTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑CAGCAgg                                                                                                 >  2:700512/1‑140 (MQ=255)                                       tGGGGCGATGCGTATCATCTCATCCCCAGACTCAAGGGTGCCTGCAACGGCTGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTTCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGtt                                                                                               <  1:1559986/140‑1 (MQ=255)                                        ggggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑CAACAgggtt                                                                                              >  2:1530345/1‑137 (MQ=255)                                         gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCt                                                                                             <  1:278813/140‑1 (MQ=255)                                         gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCt                                                                                             <  2:892448/140‑1 (MQ=255)                                          ggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑CAGCAGGGTCtt                                                                                            >  2:70466/1‑139 (MQ=255)                                               tGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCaaaa                                                                                       >  1:1416573/1‑139 (MQ=255)                                                   tATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAACAGGGTCTGCAAAAAATg                                                                                   >  1:775803/1‑140 (MQ=255)                                                      catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAg                                                                                <  1:236037/140‑1 (MQ=255)                                                      catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAg                                                                                <  1:1145996/140‑1 (MQ=255)                                                         cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGt                                                                             <  2:1101192/140‑1 (MQ=255)                                                         cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGt                                                                             <  2:578463/140‑1 (MQ=255)                                                         cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCCTGCGGGGGGGTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGGTCAAAGAAGGTT‑TTCAAAAAATTCAGccg                                                                             >  2:236037/1‑138 (MQ=255)                                                          acaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTg                                                                            <  2:1857438/140‑1 (MQ=255)                                                           caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑AAACAGGGTCTTCAAAAAAAGCAGCgggg                                                                           >  2:1539541/1‑137 (MQ=255)                                                              ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCCTTTTTCACCTGTTTTTCCTGGGGGGGGCCTTTTACCAGCTGCCGCGTCCTGGCGGGAAGGGTC‑CCTCAGGCTCTTCTAATAATGCCGCGGGAcgg                                                                        >  2:1145996/1‑138 (MQ=255)                                                                ccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCa                                                                      <  2:746609/140‑1 (MQ=255)                                                                ccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAGGGTTTGCAAAAAAAGCAACGTGAATGCa                                                                      >  2:360572/1‑140 (MQ=255)                                                                 cAGACTCAAGGGTGTCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTC‑AAGCAGGGTCCTCAAATAATGCAGCGGGACTGCAg                                                                     >  2:221903/1‑140 (MQ=255)                                                                         aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc                                                             <  1:1539541/140‑1 (MQ=255)                                                                          gggTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCCGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGTGGTC‑AAGCAGGTTCTTCCAAAAATGCAGCGTGGCTGCAGGGGCCAccg                                                            >  2:1444603/1‑140 (MQ=255)                                                                           ggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcgg                                                           <  1:1594823/140‑1 (MQ=255)                                                                             tGCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAGGGTTTTCTAATAATGCAGCGTGGCTGCCGGGGCCAgcgggg                                                         >  2:479737/1‑138 (MQ=255)                                                                             tGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCCTGTGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑AAGCAGGGTCCTCCAAAAAAGCAGCGGGAATGCCGGGGCCAgcggcg                                                         >  2:1461013/1‑140 (MQ=255)                                                                              gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑AAACAGGGTTTTCAAAAAATGCAGCGTGGCTGCAGGGACCAGCGgcgc                                                        >  2:1593792/1‑140 (MQ=255)                                                                              gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑CAGCAGGTTTTTCAAAAAATGCAGCGGGACTGCCGGGGCCAGCGgcgc                                                        >  2:721004/1‑140 (MQ=255)                                                                               ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAGGTTTTGCAAAAAAAGCAGCGTGACTGCAGGGGCCAGCGGcgcg                                                       >  1:1126706/1‑140 (MQ=255)                                                                               ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑AAGCAGGGTTTTCAAAAAATGCAGCGTGACTGCCGGGGCCCGCGGggcg                                                       >  2:1337876/1‑140 (MQ=255)                                                                                 tGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                                     <  2:299707/140‑1 (MQ=255)                                                                                      cggcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCg                                                <  1:1593792/140‑1 (MQ=255)                                                                                      cggcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCg                                                <  1:221903/140‑1 (MQ=255)                                                                                      cggcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAGGGTTTGCAAATAAAGCAGCGTGACTGCAGGGGCCAGCGGCGCGATAGCCg                                                >  2:1687191/1‑140 (MQ=255)                                                                                                          aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGTGCGATAGCCGTGAGATCGTTTTTCGCcagc                            <  2:1728008/140‑1 (MQ=255)                                                                                                          aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                            >  2:688220/1‑140 (MQ=255)                                                                                                              ttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGtt                        >  1:1619979/1‑140 (MQ=255)                                                                                                               tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTg                       <  1:688220/140‑1 (MQ=255)                                                                                                                     gTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAc                 <  1:1964519/140‑1 (MQ=255)                                                                                                                        tctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGa                      <  2:184246/132‑1 (MQ=255)                                                                                                                        tctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGa                      >  1:184246/1‑132 (MQ=255)                                                                                                                                gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCgg      >  2:634086/1‑140 (MQ=255)                                                                                                                                gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCgg      >  2:1349124/1‑140 (MQ=255)                                                                                                                                gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCgg      >  1:311179/1‑140 (MQ=255)                                                                                                                                 ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAc                 <  2:1429762/128‑1 (MQ=255)                                                                                                                                 ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAc                 >  1:1429762/1‑128 (MQ=255)                                                                                                                                 ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCggg     >  1:1635032/1‑139 (MQ=255)                                                                                                                                  gggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGat    >  1:811385/1‑140 (MQ=255)                                                                                                                                  gggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGat    >  1:547713/1‑140 (MQ=255)                                                                                                                                   ggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGata   >  2:1377567/1‑140 (MQ=255)                                                                                                                                    gTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGatat  >  1:1668534/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                            CTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT  >  NC_000913/650677‑650947

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