Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,208,051	G→T	intergenic (‑519/‑96)	purH ← / → rrsE	bifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,208,051	0	G	T	100.0%	249.9 / NA	72	intergenic (‑519/‑96)	purH/rrsE	bifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (37/35);  total (37/35)

AAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTC  >  NC_000913/4207918‑4208128                                                                                                                                        |                                                                              aaaaTCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCt                                                                         >  2:42187/1‑140 (MQ=255)          aaaTTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGc                                                                >  1:1995786/1‑140 (MQ=255)            aTTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCaa                                                              <  2:756311/140‑1 (MQ=255)            aTTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCaa                                                              <  2:606855/140‑1 (MQ=255)               cAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGAc                                                           <  1:549833/140‑1 (MQ=255)               cAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGAc                                                           >  1:1313469/1‑140 (MQ=255)                 gggTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGa                                                         <  1:1347709/140‑1 (MQ=255)                  ggTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGaa                                                        <  2:728540/140‑1 (MQ=255)                  ggTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGaa                                                        <  2:488147/140‑1 (MQ=255)                    ttGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGaaaa                                                      >  2:2006630/1‑140 (MQ=255)                     tGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGaaaaa                                                     <  2:1152658/140‑1 (MQ=255)                       aCTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATg                                                   >  1:1686666/1‑140 (MQ=255)                       aCTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATg                                                   >  2:2281368/1‑140 (MQ=255)                       aCTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATg                                                   >  2:362123/1‑140 (MQ=255)                       aCTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATg                                                   >  2:1148338/1‑140 (MQ=255)                            gAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                              <  1:1323252/140‑1 (MQ=255)                            gAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                              <  1:1265181/140‑1 (MQ=255)                            gAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                              <  2:1987512/140‑1 (MQ=255)                            gAAAGAGGAAAGCGTAATATACGC‑‑CCTCGCGACAGTGCGCTGAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCa                                            <  2:1245166/140‑1 (MQ=11)                              aaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCa                                            <  1:2069236/140‑1 (MQ=255)                              aaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCa                                            <  2:1686666/140‑1 (MQ=255)                                 agGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGt                                         <  1:2006630/140‑1 (MQ=255)                                  ggAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtc                                        >  1:325415/1‑140 (MQ=255)                                   gAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                       >  1:2577039/1‑140 (MQ=255)                                   gAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                       <  2:2577039/140‑1 (MQ=255)                                    aaaGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtctc                                      <  2:1995786/140‑1 (MQ=255)                                      aGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCaa                                    >  2:1009077/1‑140 (MQ=255)                                       gCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAag                                   >  1:938827/1‑140 (MQ=255)                                       gCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAag                                   >  1:1147654/1‑140 (MQ=255)                                          tAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGt                                <  1:1009077/140‑1 (MQ=255)                                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTg                               <  2:1775649/140‑1 (MQ=255)                                            atatACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                              <  1:2278684/140‑1 (MQ=255)                                            atatACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                              <  1:1148338/140‑1 (MQ=255)                                                 cGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACg                         <  1:2281368/140‑1 (MQ=255)                                                 cGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACg                         <  2:1485025/140‑1 (MQ=255)                                                  gCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAacac                          >  2:991935/1‑140 (MQ=255)                                                  gCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACg                         <  2:2035580/139‑1 (MQ=255)                                                  gCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACg                         >  1:2035580/1‑139 (MQ=255)                                                  gCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGt                        >  1:1129492/1‑140 (MQ=255)                                                  gCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGt                        >  1:1129493/1‑140 (MQ=255)                                                   ccACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTa                       >  1:2091374/1‑140 (MQ=255)                                                   ccACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTa                       >  2:1680523/1‑140 (MQ=255)                                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAac                            <  1:248373/131‑1 (MQ=255)                                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAac                            >  2:248373/1‑131 (MQ=255)                                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTc                   >  2:893302/1‑140 (MQ=255)                                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTc                   >  2:1716220/1‑140 (MQ=255)                                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTc                   >  1:1291249/1‑140 (MQ=255)                                                       atCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAAtcc                   >  1:1701309/2‑138 (MQ=255)                                                         cgcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCAt                 >  2:1969846/1‑140 (MQ=255)                                                           cgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa               <  2:1584487/140‑1 (MQ=255)                                                           cgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa               >  1:1204528/1‑140 (MQ=255)                                                           cgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa               <  1:1716220/140‑1 (MQ=255)                                                           cgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa               >  1:911838/1‑140 (MQ=255)                                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGTACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATGCATTACg             <  2:2494428/140‑1 (MQ=255)                                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             <  1:1803051/140‑1 (MQ=255)                                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             <  2:2159598/140‑1 (MQ=255)                                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             <  2:1129493/140‑1 (MQ=255)                                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             <  1:2158374/140‑1 (MQ=255)                                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             <  1:1488120/140‑1 (MQ=255)                                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             <  2:1129492/140‑1 (MQ=255)                                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt         >  2:583016/1‑140 (MQ=38)                                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt         >  2:1323091/1‑140 (MQ=38)                                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt         >  1:848865/1‑140 (MQ=38)                                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt         >  2:1985065/1‑140 (MQ=38)                                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt         >  2:1816771/1‑140 (MQ=38)                                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt         >  1:2216574/1‑140 (MQ=38)                                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt         <  1:42187/140‑1 (MQ=38)                                                                   cgcTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATGCCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt       >  1:445975/1‑140 (MQ=255)                                                                    gcTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             <  2:884038/133‑1 (MQ=25)                                                                    gcTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg             >  1:884038/1‑133 (MQ=25)                                                                      tAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAt    >  1:819430/1‑140 (MQ=18)                                                                        aaGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACTGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTc  >  1:1873119/1‑140 (MQ=17)                                                                                                                                        |                                                                              AAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATC‑‑TGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTC  >  NC_000913/4207918‑4208128

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