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

CTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT  >  NC_000913/650670‑650947                                                                                                                                            |                                                                                                                                            ctcAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGt                                                                                                                                             <  2:565252/140‑1 (MQ=255)    aaCGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCtt                                                                                                                                          >  1:896158/1‑140 (MQ=255)     aCGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGGTGGGTTATTTTCATTTTTTACCTGTTTTTCATGCGGGGGGCttt                                                                                                                                         >  2:1345006/1‑140 (MQ=255)     aCGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCttt                                                                                                                                         >  2:1434966/1‑140 (MQ=255)        ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGa                                                                                                                                      >  2:944484/1‑140 (MQ=255)            tATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAg                                                                                                                                  <  2:896158/140‑1 (MQ=255)                   ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCcgc                                                                                                                           >  1:586995/1‑140 (MQ=255)                   ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCcgc                                                                                                                           >  1:1022973/1‑140 (MQ=255)                   ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCCGGTTTTTTTTCATTTTTTAACCGGCTCTCATGCGGGGGTCATTTTGCAAGCTGcccc                                                                                                                           >  2:610064/1‑138 (MQ=255)                     gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGt                                                                                                                         >  2:1259316/1‑140 (MQ=255)                     gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTTACGAGCTGCCGCGt                                                                                                                         >  1:1544875/1‑140 (MQ=255)                     gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCcgcgg                                                                                                                         >  1:491058/1‑139 (MQ=255)                               atCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGa                                                                                                               <  1:1541651/140‑1 (MQ=255)                               atCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGa                                                                                                               <  1:1434966/140‑1 (MQ=255)                                    gCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCt                                                                                                          <  1:1345006/140‑1 (MQ=255)                                        aTCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTC‑aaa                                                                                                     >  1:582162/1‑139 (MQ=255)                                         gcGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGc                                                                                                    <  1:55091/139‑1 (MQ=255)                                          cGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCa                                                                                                   >  1:1566724/1‑140 (MQ=255)                                           gAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAg                                                                                                  <  1:944484/140‑1 (MQ=255)                                           gAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAg                                                                                                  <  1:586621/140‑1 (MQ=255)                                           gAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTCTTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCCC‑cagcag                                                                                                  >  2:586621/1‑140 (MQ=255)                                            aGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAgg                                                                                                 >  2:522340/1‑140 (MQ=255)                                             gTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGt                                                                                                <  1:580476/140‑1 (MQ=255)                                                     aTGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAGGGTCTGCaac                                                                                        >  2:55091/1‑139 (MQ=255)                                                          tATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATg                                                                                   <  2:1508841/140‑1 (MQ=255)                                                           atcatcACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAACAGGGTCTGCAAATAATGc                                                                                  >  2:1107397/1‑140 (MQ=255)                                                            tcatcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTC‑AAGACGGGTTTGCAAATAATGCa                                                                                 >  1:1525357/1‑140 (MQ=255)                                                             catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAg                                                                                <  2:11656/140‑1 (MQ=255)                                                                cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑AAGCAGGGTTTGCCAATAATTCAGCgg                                                                             >  2:570956/1‑139 (MQ=255)                                                                  caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                           <  1:1107397/140‑1 (MQ=255)                                                                     ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCCC‑AAGCAGGGTTTGCAAAAAAATCAGCGGGACTg                                                                        >  2:1579843/1‑140 (MQ=255)                                                                       ccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAACAGGGTTTTCAAATAATTCAGCGCGGCTGcc                                                                      >  2:175452/1‑139 (MQ=255)                                                                            cTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGAc                                                                 <  1:1240286/140‑1 (MQ=255)                                                                                     gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGGCCAGCGgcgc                                                        >  2:1554123/1‑140 (MQ=255)                                                                                      ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑CAACAGGGTTTTCAAAAAATGCCGCGGGGCTGCAGGGGCCAGCGGcgcg                                                       >  2:1524847/1‑140 (MQ=255)                                                                                      ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTC‑AAGCAGGGTTTTCAAATAAATCAGCGTGGCTGCAGGGGCCAGCGGggcg                                                       >  2:835543/1‑140 (MQ=255)                                                                                      ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCATTTTCTCATGCAGGGGACTTTTGACGCGCCGCCCCGTCCCGGCGGGGGTGCGCAAAACAGATTTTCCAAA‑AATCGCGCGCGACTCCTGGGGCCAGCGgcgcc                                                       >  2:1117335/1‑139 (MQ=255)                                                                                        tGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                                     <  1:1579843/140‑1 (MQ=255)                                                                                         gCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTC‑AAGCAGGGTTTGCAAATAAAGCAGCGCGGCTGCAGGGGCCAGCGGCGCGATa                                                    >  2:466309/1‑140 (MQ=255)                                                                                             cggcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTC‑AAGCAGGGTTTGCAAATAATGCAGCGGGACTGTAGGGGCCCGCGGCGCGATAGCCg                                                >  2:1197554/1‑140 (MQ=255)                                                                                                             tttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGc                                <  1:1117335/140‑1 (MQ=255)                                                                                                             tttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGc                                <  1:1259316/140‑1 (MQ=255)                                                                                                             tttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGc                                <  1:835543/140‑1 (MQ=255)                                                                                                              ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGTAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                               <  1:1197554/140‑1 (MQ=255)                                                                                                              ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                               <  1:175452/140‑1 (MQ=255)                                                                                                              ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                               <  2:1097870/140‑1 (MQ=255)                                                                                                                 aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                            >  2:972439/1‑140 (MQ=255)                                                                                                                 aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                            >  2:552664/1‑140 (MQ=255)                                                                                                                 aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                            >  2:1089477/1‑140 (MQ=255)                                                                                                                 aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                            >  2:764518/1‑140 (MQ=255)                                                                                                                  tttttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGTAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGGTAGCCGTGAGATCGTTTTTCGCcagca                           >  1:875305/1‑140 (MQ=255)                                                                                                                  tttttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagca                           <  2:1566724/140‑1 (MQ=255)                                                                                                                   ttttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCGTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagcag                          <  1:570956/140‑1 (MQ=255)                                                                                                                      tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTg                       >  2:610696/1‑140 (MQ=255)                                                                                                                      tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTg                       >  2:929344/1‑140 (MQ=255)                                                                                                                             tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc                >  2:1241266/1‑140 (MQ=255)                                                                                                                               tctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCg              >  2:763877/1‑140 (MQ=255)                                                                                                                                ctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCgg             <  1:552664/140‑1 (MQ=255)                                                                                                                                    tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg         >  1:827740/1‑140 (MQ=255)                                                                                                                                    tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg         >  1:380289/1‑140 (MQ=255)                                                                                                                                       gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCgg      >  2:970533/1‑140 (MQ=255)                                                                                                                                        ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATCCAGCGCGACTTCAGGGACCAACGGACCGATAGCCGTGAGAGCGTTTTTCGCCCGCAAATGAAGGAGCCGGGAAGAGGa     >  2:1541485/1‑140 (MQ=255)                                                                                                                                        ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGa     >  2:994285/1‑140 (MQ=255)                                                                                                                                         gggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGat    <  2:586995/140‑1 (MQ=255)                                                                                                                                           gTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGatat  >  1:1580444/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                            CTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTC‑AAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT  >  NC_000913/650670‑650947

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