Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,648,042	T→C	intergenic (‑254/‑486)	dinQ ← / → arsR	UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,042	0	T	C	100.0%	325.5 / NA	89	intergenic (‑254/‑486)	dinQ/arsR	UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (51/38);  total (51/38)

TGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647947‑3648138                                                                                                |                                                                                                  tGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATgg                                                                                               <  1:1726574/100‑1 (MQ=255)  gcgcATTGCCCGGGGGAAAATGCCCCCGAGCTAACCGGGTAATCATACAGTACATGTAGGGTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATgg                                                                                               >  2:1308333/1‑99 (MQ=255)  gcgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGt                                                                                              >  2:995558/1‑100 (MQ=255)  gcgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGt                                                                                              >  2:1393315/1‑100 (MQ=255)  gcgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAATACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGt                                                                                              >  2:1309646/1‑100 (MQ=255)     cATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTAt                                                                                            <  2:1942642/99‑1 (MQ=255)        tGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAg                                                                                         <  2:1586263/99‑1 (MQ=255)         gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                       >  1:899660/1‑100 (MQ=255)           cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAtt                                                                                     >  1:425783/1‑100 (MQ=255)           cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAtt                                                                                     >  1:787723/1‑100 (MQ=255)           cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAtt                                                                                     <  1:1133376/100‑1 (MQ=255)            agagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTg                                                                                    <  2:1392301/100‑1 (MQ=255)              agTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGAt                                                                                  >  2:1865699/1‑100 (MQ=255)              agTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGAt                                                                                  <  2:1879577/100‑1 (MQ=255)               gTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATg                                                                                 >  2:1233259/1‑100 (MQ=255)                 gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                               >  2:1423790/1‑100 (MQ=255)                 gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                               >  2:716811/1‑100 (MQ=255)                  aaaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                               <  2:608051/99‑1 (MQ=255)                   aaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGtt                                                                             <  2:348084/100‑1 (MQ=255)                    aaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTTCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTa                                                                            <  2:1536004/100‑1 (MQ=255)                     aTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAg                                                                           <  2:820170/100‑1 (MQ=255)                       gCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgat                                                                         >  2:1761521/1‑100 (MQ=255)                       gCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAga                                                                          <  1:1393324/99‑1 (MQ=255)                         cTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatga                                                                       >  1:1563441/1‑100 (MQ=255)                           cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                     <  2:1502424/100‑1 (MQ=255)                               gCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTa                                                                 <  1:356169/100‑1 (MQ=255)                               gCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTa                                                                 >  1:1681777/1‑100 (MQ=255)                                 tAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtc                                                               >  1:1636885/1‑100 (MQ=255)                                  aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                              >  1:1558983/1‑100 (MQ=255)                                    cTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCa                                                            >  2:1033070/1‑100 (MQ=255)                                    cTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCa                                                            <  1:1645478/100‑1 (MQ=255)                                    cTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCa                                                            >  2:1171122/1‑100 (MQ=255)                                     tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAc                                                           >  1:941032/1‑100 (MQ=255)                                        aTAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTcc                                                        <  2:662570/100‑1 (MQ=255)                                           aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGt                                                     >  2:1473026/1‑100 (MQ=255)                                           aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGt                                                     >  1:1717068/1‑100 (MQ=255)                                            tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc                                                    >  2:1618266/1‑100 (MQ=255)                                             cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                   <  2:662154/100‑1 (MQ=255)                                              aTACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCag                                                  >  1:614059/1‑100 (MQ=255)                                               tACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCaga                                                 <  2:853859/100‑1 (MQ=255)                                                aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                                <  2:787720/100‑1 (MQ=255)                                                aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                                <  1:1761535/100‑1 (MQ=255)                                                  aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcc                                              >  1:61568/1‑100 (MQ=255)                                                   gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcca                                             >  2:825799/1‑100 (MQ=255)                                                   gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcca                                             >  2:1908917/1‑100 (MQ=255)                                                      cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                          >  1:742850/1‑100 (MQ=255)                                                      cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                          >  1:1630771/1‑100 (MQ=255)                                                      cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                          >  1:15441/1‑100 (MQ=255)                                                        tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                        <  1:722099/100‑1 (MQ=255)                                                        tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                        >  2:437069/1‑100 (MQ=255)                                                         gCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACt                                       >  1:891202/1‑100 (MQ=255)                                                          cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTc                                      <  1:1830017/100‑1 (MQ=255)                                                          cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTc                                      <  1:1828209/100‑1 (MQ=255)                                                            ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                    >  1:2011007/1‑100 (MQ=255)                                                            ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                    <  2:1584548/100‑1 (MQ=255)                                                              ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGCCTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAggg                                  >  1:84423/1‑100 (MQ=255)                                                               tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                                 >  1:1236531/1‑100 (MQ=255)                                                               tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                                 <  1:944623/100‑1 (MQ=255)                                                               atctaaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                                 <  2:1563815/96‑1 (MQ=255)                                                                ataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCt                                <  1:1233265/100‑1 (MQ=255)                                                                ataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCt                                >  1:1165707/1‑100 (MQ=255)                                                                 taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCt                                >  2:69229/1‑99 (MQ=255)                                                                  aaaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTGAGAGCCACCAACTCAGGGCTgg                              >  2:266666/1‑100 (MQ=255)                                                                  aaaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTgg                              <  2:1636039/100‑1 (MQ=255)                                                                  aaaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTAATAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTgg                              >  2:1091043/1‑100 (MQ=255)                                                                   aaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGcagat                             >  1:1563825/1‑96 (MQ=255)                                                                    aaCCAGCACTTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaa                            <  1:716813/100‑1 (MQ=255)                                                                    aaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCACCTCAGGGCTGGaa                            >  2:1877918/1‑100 (MQ=255)                                                                    aaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaa                            <  1:702991/100‑1 (MQ=255)                                                                      ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                          <  2:425782/100‑1 (MQ=255)                                                                       cAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGt                         <  1:1877935/100‑1 (MQ=255)                                                                            cGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑aaaaa                   >  1:56688/1‑100 (MQ=255)                                                                               ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAAcc                >  2:177403/1‑100 (MQ=255)                                                                               ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAAcc                <  2:1391980/100‑1 (MQ=255)                                                                               ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAAcc                <  1:995562/100‑1 (MQ=255)                                                                                cTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCg               <  2:1558973/100‑1 (MQ=255)                                                                                 ttGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGa              >  1:1689404/1‑100 (MQ=255)                                                                                 ttGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGa              <  1:50606/100‑1 (MQ=255)                                                                                  tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCCGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGAc             <  2:1092702/100‑1 (MQ=255)                                                                                   gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACg            <  2:742847/100‑1 (MQ=255)                                                                                     aaTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCa          >  2:157036/1‑100 (MQ=255)                                                                                     aaTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCa          >  1:533267/1‑100 (MQ=255)                                                                                     aaTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCa          >  1:895811/1‑100 (MQ=255)                                                                                       tAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCaaa        <  2:899657/100‑1 (MQ=255)                                                                                         gTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGt      >  2:540560/1‑100 (MQ=255)                                                                                         gTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGt      >  2:941375/1‑100 (MQ=255)                                                                                          tttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGActttgcgt      >  1:1697425/1‑93 (MQ=25)                                                                                          tttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTc     >  1:339947/1‑100 (MQ=255)                                                                                             cAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGt  >  1:1142908/1‑100 (MQ=255)                                                                                                |                                                                                                  TGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647947‑3648138

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