Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,871,978	A→G	R171R (CGT→CGC)	ispD ←	2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,871,978	0	A	G	100.0%	235.3 / NA	64	R171R (CGT→CGC)	ispD	2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/33);  total (31/33)

GGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTG  >  NC_000913/2871841‑2872113                                                                                                                                          |                                                                                                                                        gggCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTcgcg                                                                                                                                       >  1:775387/1‑140 (MQ=255)        tGACTTTAATGTTATCCGCACGGCCTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGaaa                                                                                                                                >  2:663434/1‑140 (MQ=255)              tAATGTTATCCGCACGGCCTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTgc                                                                                                                          >  2:37920/1‑140 (MQ=255)                   ttATCCGCACGGCATTCGACCAACTGAGGATGGAGTCCACAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGt                                                                                                                     <  2:291646/140‑1 (MQ=255)                     aTCCGCACGGCCTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCa                                                                                                                   >  2:6643/1‑140 (MQ=255)                         gCACGGCCTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAgcgc                                                                                                               <  2:775387/140‑1 (MQ=255)                               ccTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCa                                                                                                         >  1:398620/1‑140 (MQ=255)                                   cGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAg                                                                                                     <  2:297057/140‑1 (MQ=255)                                         aCTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTg                                                                                               >  2:302492/1‑140 (MQ=255)                                           tGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCg                                                                                             <  2:909319/140‑1 (MQ=255)                                             aGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGAt                                                                                           >  2:108790/1‑140 (MQ=255)                                              ggATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATc                                                                                          >  1:457126/1‑140 (MQ=255)                                              ggATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATc                                                                                          <  1:6643/140‑1 (MQ=255)                                                aTGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCaa                                                                                        <  2:341947/140‑1 (MQ=255)                                                 tGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAAc                                                                                       <  2:173420/140‑1 (MQ=255)                                                    aaTCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAg                                                                                                     >  1:991859/1‑123 (MQ=255)                                                    aaTCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAg                                                                                                     <  2:991859/123‑1 (MQ=255)                                                    aaTCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCaa                                                                                        <  1:41624/136‑1 (MQ=255)                                                    aaTCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCaa                                                                                        >  2:41624/1‑136 (MQ=255)                                                    aaTCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGt                                                                                    >  1:956870/1‑140 (MQ=255)                                                         gCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCaa                                                                                        <  2:875818/131‑1 (MQ=255)                                                         gCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCaa                                                                                        >  1:875818/1‑131 (MQ=255)                                                         gCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAg                                                                               <  1:108790/140‑1 (MQ=255)                                                          cAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGc                                                                              >  2:811651/1‑140 (MQ=255)                                                          cAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGc                                                                              >  1:58963/1‑140 (MQ=255)                                                                       ccGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTg                                                                 >  2:504247/1‑140 (MQ=255)                                                                         gAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGcc                                                               <  1:951237/140‑1 (MQ=255)                                                                         gAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGcc                                                               >  2:917704/1‑140 (MQ=255)                                                                          aGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGccc                                                              >  1:541528/1‑140 (MQ=255)                                                                               tcgtcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTc                                                         >  1:833914/1‑140 (MQ=255)                                                                                cgtcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCg                                                        <  2:541528/140‑1 (MQ=255)                                                                                cgtcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCg                                                        <  1:302492/140‑1 (MQ=255)                                                                                  tcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGc                                                      <  2:58963/140‑1 (MQ=255)                                                                                  tcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGc                                                      <  1:663434/140‑1 (MQ=255)                                                                                   cgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGTGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCa                                                     <  1:260541/140‑1 (MQ=255)                                                                                   cgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCa                                                     <  1:399166/140‑1 (MQ=255)                                                                                   cgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCa                                                     <  1:37920/140‑1 (MQ=255)                                                                                   cgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCa                                                     <  2:398620/140‑1 (MQ=255)                                                                                   cgGTAATAGTCGCGCCTCCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCa                                                     <  2:980402/140‑1 (MQ=255)                                                                                     gTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   <  2:314867/140‑1 (MQ=255)                                                                                     gTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   <  2:833914/140‑1 (MQ=255)                                                                                     gTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   <  1:504247/140‑1 (MQ=255)                                                                                     gTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   <  2:285647/140‑1 (MQ=255)                                                                                            tCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   >  2:615150/1‑133 (MQ=255)                                                                                            tCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   <  1:615150/133‑1 (MQ=255)                                                                                            tCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATa                                            >  2:922720/1‑140 (MQ=255)                                                                                            tCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATa                                            >  2:331491/1‑140 (MQ=255)                                                                                              gcgcCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   <  2:358057/131‑1 (MQ=255)                                                                                              gcgcCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACg                                                   >  1:358057/1‑131 (MQ=255)                                                                                               cgcCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTa                                         >  2:1027895/1‑140 (MQ=255)                                                                                                gcCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCa                                              <  2:437485/134‑1 (MQ=255)                                                                                                gcCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCa                                              >  1:437485/1‑134 (MQ=255)                                                                                                       tttAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCAc                                 <  1:811651/140‑1 (MQ=255)                                                                                                       tttAGAGCGCGCGTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCAc                                 <  1:922720/140‑1 (MQ=255)                                                                                                               gcgcgTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTgcgg                         >  1:1084966/1‑140 (MQ=255)                                                                                                               gcgcgTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTgcgg                         <  1:1082719/140‑1 (MQ=255)                                                                                                                cgcgTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTgcggc                        >  1:1074878/1‑140 (MQ=255)                                                                                                                cgcgTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTgcggc                        >  2:829963/1‑140 (MQ=255)                                                                                                                 gcgTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTgcggcg                       <  1:1027895/140‑1 (MQ=255)                                                                                                                 gcgTCAGACAGTCATGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTgcggcg                       <  1:720035/140‑1 (MQ=255)                                                                                                                               tGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGt         >  2:696395/1‑140 (MQ=255)                                                                                                                               tGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGt         >  2:354525/1‑140 (MQ=255)                                                                                                                               tGTAACAGCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGt         >  1:969385/1‑140 (MQ=255)                                                                                                                                      gCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTg  >  2:1032245/1‑140 (MQ=255)                                                                                                                                      gCTCGCGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCATTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGCTCCCCCCCGTGCGGCTg  >  2:980654/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        GGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCAACTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTG  >  NC_000913/2871841‑2872113

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