Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,621,542	G→A	D71N (GAC→AAC)	ydeE →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,542	0	G	A	100.0%	239.5 / NA	73	D71N (GAC→AAC)	ydeE	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/43);  total (30/43)

CTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATT  >  NC_000913/1621410‑1621662                                                                                                                                     |                                                                                                                         cTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATTGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    <  1:706141/139‑1 (MQ=255) cTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAagaaga                                                                                                                    <  1:1045352/139‑1 (MQ=255) cTGCCATTTATCACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    <  2:1167330/139‑1 (MQ=255)  tGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaa                                                                                                                   <  2:1567777/139‑1 (MQ=255)   gCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                  <  1:413071/139‑1 (MQ=255)      aTTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGc                                                                                                               <  1:1541274/139‑1 (MQ=255)              ccATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATACGATTACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGtt                                                                                                       >  1:1567656/1‑139 (MQ=255)                        gAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAAtt                                                                                             >  2:550965/1‑139 (MQ=255)                         aGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTa                                                                                            <  1:550965/139‑1 (MQ=255)                         aGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTa                                                                                            >  1:508910/1‑139 (MQ=255)                          atCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaa                                                                                                                   <  2:156179/114‑1 (MQ=255)                          atCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaa                                                                                                                   >  1:156179/2‑115 (MQ=255)                           tCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                <  1:497836/117‑1 (MQ=255)                           tCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                >  2:497836/1‑117 (MQ=255)                            cGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCg                                                                                         >  2:1409029/1‑139 (MQ=255)                                         tGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTAATGGCAATTACCGCCTTCGCCAg                                                                               >  2:901654/1‑136 (MQ=255)                                         tGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTAATGGCAATTACCGCCTTCGCCAg                                                                               <  1:901654/136‑1 (MQ=255)                                          gAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGtt                                                                                                       >  1:1713110/1‑111 (MQ=255)                                          gAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGtt                                                                                                       <  2:1713110/111‑1 (MQ=255)                                           aGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    <  1:382314/97‑1 (MQ=255)                                           aGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    >  2:382314/1‑97 (MQ=255)                                           aGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTatat                                                                                                          >  2:178669/1‑107 (MQ=255)                                           aGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTatat                                                                                                          <  1:178669/107‑1 (MQ=255)                                            gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGttt                                                                         >  2:1528243/1‑139 (MQ=255)                                            gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGttt                                                                         >  1:1031788/1‑139 (MQ=255)                                              gtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                       <  2:1567656/139‑1 (MQ=255)                                                cGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTAtt                                                                     <  2:582594/139‑1 (MQ=255)                                                cGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTAtt                                                                     <  2:959869/139‑1 (MQ=255)                                                 gATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAAc                                                                                                                        >  1:1483864/1‑87 (MQ=255)                                                 gATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTg                                                                    >  1:1665987/1‑139 (MQ=255)                                                  aTCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                   <  1:513558/139‑1 (MQ=255)                                                  aTCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                   <  2:630217/139‑1 (MQ=255)                                                    cTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                 <  2:508910/139‑1 (MQ=255)                                                    cTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                 <  1:1138669/139‑1 (MQ=255)                                                     tAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                <  1:44099/91‑1 (MQ=255)                                                     tAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                >  2:44099/1‑91 (MQ=255)                                                       aTCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTa                                                              <  2:471383/139‑1 (MQ=255)                                                             tATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAg                                                        >  1:1085095/1‑139 (MQ=255)                                                              aTGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGt                                                       <  1:1467514/139‑1 (MQ=255)                                                                 cGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGaa                                                    <  2:1435272/139‑1 (MQ=255)                                                                      aCAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                 <  2:1707342/121‑1 (MQ=255)                                                                      aCAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                 >  1:1707342/1‑121 (MQ=255)                                                                       cAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCAGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGt                                              <  2:1085095/139‑1 (MQ=255)                                                                          ttGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGAc                                           >  1:1366776/1‑139 (MQ=255)                                                                               cTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTgg                                      <  2:1722938/139‑1 (MQ=255)                                                                                  aCTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTtg                                   >  2:370505/1‑139 (MQ=255)                                                                                     aTTCGCGTCGTTTTTAGCCTCGGTTTTTGTATCCTGGCGGATAAGTTCAACAAGAATCGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGc                                <  2:1156691/139‑1 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCtat                                                                                                            >  2:1151265/1‑56 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCtat                                                                                                            <  1:1151265/56‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAg                                                                               <  2:403992/83‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAg                                                                               >  1:403992/1‑83 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCAt                                                                >  2:976444/1‑98 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCAt                                                                <  1:976444/98‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACg                                               <  1:894259/115‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACg                                               >  2:894259/1‑115 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                       <  2:159214/139‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                       >  1:1636786/1‑139 (MQ=255)                                                                                                ttttAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGcc                     <  2:531718/139‑1 (MQ=255)                                                                                                     gCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCAt                <  1:1528243/139‑1 (MQ=255)                                                                                                        tCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGt                                              <  2:1585323/106‑1 (MQ=255)                                                                                                        tCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGt                                              >  1:1585323/1‑106 (MQ=255)                                                                                                               tGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGcc      >  2:1408066/1‑139 (MQ=255)                                                                                                                  tATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTAt   <  2:1031788/139‑1 (MQ=255)                                                                                                                   aTCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTAtt  <  1:1408066/139‑1 (MQ=255)                                                                                                                   aTCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCAGTAACTGCGCCTAtt  <  1:1409029/139‑1 (MQ=255)                                                                                                                    tCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTa                                                                                                      <  1:1592879/38‑1 (MQ=38)                                                                                                                    tCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTa                                                                                                      >  2:1592879/1‑38 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAg                                                                               >  2:1418629/1‑59 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAg                                                                               <  1:1418629/59‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGtt                                    >  2:1683092/1‑102 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGtt                                    <  1:1683092/102‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAG‑AACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                 <  2:1262382/72‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAG‑AACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                 >  1:1262382/1‑72 (MQ=255)                                                                                                                                     |                                                                                                                         CTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATT  >  NC_000913/1621410‑1621662

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