Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,490,055	T→G	100%	E3D (GAA→GAC)	yjgR ←	DUF853 family protein with NTPase fold

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,490,055	0	T	G	100.0%	322.2 / ‑6.7	97	E3D (GAA→GAC)	yjgR	DUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base T (2/0);  major base G (41/52);  minor base A (2/0);  total (45/52)

TACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGC  >  NC_000913/4489958‑4490147                                                                                                  |                                                                                             tACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTc                                                                                            <  2:951759/100‑1 (MQ=255)    ccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACt                                                                                         >  2:1910672/1‑100 (MQ=255)     cGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACt                                                                                         <  1:83196/99‑1 (MQ=255)     cGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTc                                                                                        <  1:1442621/100‑1 (MQ=255)      gTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCa                                                                                       <  1:2219993/100‑1 (MQ=255)      gTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCa                                                                                       <  1:2195769/100‑1 (MQ=255)      gTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCa                                                                                       <  1:2121853/100‑1 (MQ=255)         cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                    >  2:1178561/1‑100 (MQ=255)          ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                   <  2:2254527/100‑1 (MQ=255)          ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                   >  1:2168165/1‑100 (MQ=255)          ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                   <  2:33735/100‑1 (MQ=255)          ccGTTGCGCAGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                   <  1:903124/100‑1 (MQ=255)           cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCTGCAGTAAAAACAGTTCTGTGTCCGGCGTGCAGGCAATTAACAGGGGGTCACTCATgaga                                                                                  <  2:1491734/100‑1 (MQ=255)           cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgaga                                                                                  <  1:1759209/100‑1 (MQ=255)              tGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGatat                                                                               >  1:2377874/1‑100 (MQ=255)                 gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                             <  1:1083856/99‑1 (MQ=255)                   cttccgatctcccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                          <  1:2188407/90‑1 (MQ=255)                   cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                          <  1:493771/100‑1 (MQ=255)                   cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                          <  1:730252/100‑1 (MQ=255)                   cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctc                                                                           <  2:962647/99‑1 (MQ=255)                    ggTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGAAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGTGTCACTCATGAGATATCCTCCa                                                                         >  2:1944162/1‑100 (MQ=255)                     gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                        <  1:648883/100‑1 (MQ=255)                     gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                        <  2:1288846/100‑1 (MQ=255)                       aaTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGtt                                                                      >  2:1370288/1‑100 (MQ=255)                       aaTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGtt                                                                      >  2:1674295/1‑100 (MQ=255)                        aTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGttt                                                                     >  1:538358/1‑100 (MQ=255)                          cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                   <  1:2381618/100‑1 (MQ=255)                           aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcct                                                                  >  1:2377306/1‑100 (MQ=255)                            gCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTTCGGCGTGCGTGCAATTAACAGGGGGTCACTCATGAGTTACCCTCCAGTTTcctg                                                                 <  2:527515/100‑1 (MQ=255)                             cccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAgatcggaag                                                                >  2:2188381/1‑92 (MQ=255)                             cccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgc                                                                <  2:2230478/100‑1 (MQ=255)                             cccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgc                                                                <  2:578/100‑1 (MQ=255)                             cccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGATCACTCATGAGATATCCTCCAGTTTcctg                                                                 >  2:5525/1‑99 (MQ=255)                             cccATGACGGTTAGCCATTCCAGGCAGTAAAAACAGTTCTGTGTACGGCGGGCGGGCAATTAACAGGGGTTCACTCATGAGATATTCTCAAATTTcctg                                                                 >  2:406034/1‑99 (MQ=255)                                 tGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcct                                                              >  1:94694/1‑98 (MQ=255)                                   aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTTTGTCCGGCGTGCGGGCAATTAACAGGGGGTCATTCATGAGATATCCTCCAGTTTCCTGCCTGGa                                                           <  2:913039/99‑1 (MQ=255)                                   aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                          <  1:1235857/100‑1 (MQ=255)                                   aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                          <  2:1876132/100‑1 (MQ=255)                                    cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCCATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTTCCTGGAg                                                          >  1:184044/1‑99 (MQ=255)                                    cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                         >  1:266389/1‑100 (MQ=255)                                    cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGCGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                         <  2:1480555/100‑1 (MQ=255)                                      gTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtat                                                       >  1:1360841/1‑100 (MQ=255)                                       ttAGCCATTCCCTGCAGTAAAAACAGTTCTGTCTACGGCGTTCGGGCAATTTACAGGGGGTCACTCATGAGATTTCCTCCAGTTTCCTGCCTGGAGtata                                                      <  2:779037/100‑1 (MQ=255)                                       ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                      <  1:1928450/100‑1 (MQ=255)                                        tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                     <  2:1370348/100‑1 (MQ=255)                                         aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                    <  2:1624662/100‑1 (MQ=255)                                          gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                   <  2:1409512/100‑1 (MQ=255)                                          gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                   >  2:2011577/1‑100 (MQ=255)                                          gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                   <  1:1164148/100‑1 (MQ=255)                                          gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                   <  1:878256/100‑1 (MQ=255)                                           ctaTTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                  <  1:1236843/98‑1 (MQ=255)                                            cATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCaa                                                 <  1:1519203/100‑1 (MQ=255)                                              ttCCAGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGTCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGCAGGATGGGCaa                                                 >  2:856576/1‑98 (MQ=255)                                                cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAAc                                                >  2:2206114/1‑97 (MQ=255)                                                cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                             >  1:300474/1‑100 (MQ=255)                                                 ccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTcc                                            <  1:452226/100‑1 (MQ=255)                                                    gCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                         <  2:979878/100‑1 (MQ=255)                                                    gCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                         <  1:1235942/100‑1 (MQ=255)                                                     cAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCa                                        <  1:1427074/100‑1 (MQ=255)                                                      actAAAAACAGTTCTTTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                         >  2:1141337/3‑98 (MQ=255)                                                       gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                      >  2:50508/1‑100 (MQ=255)                                                        tAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACTGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAag                                     >  2:1285988/1‑100 (MQ=255)                                                         aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCTTCGAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGGAAaga                                    <  2:94692/100‑1 (MQ=255)                                                         aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                    <  2:167610/100‑1 (MQ=255)                                                         aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                    >  2:599141/1‑100 (MQ=255)                                                            aaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGaa                                 >  2:1002607/1‑100 (MQ=255)                                                             aCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAAt                                >  1:23885/1‑100 (MQ=255)                                                                gTTCTGTGTCCGGTGTGCGGGCAAGTAACAGGGGGTGACTCATGAGAGATCCGCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                             <  1:985455/100‑1 (MQ=255)                                                                gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                             >  2:1354502/1‑100 (MQ=255)                                                                   cTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtgt                          >  1:1986783/1‑100 (MQ=255)                                                                   cTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtgt                          >  2:292455/1‑100 (MQ=255)                                                                    tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTa                         <  2:271911/100‑1 (MQ=255)                                                                    tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGAAAAGAGAATACgtgt                          >  2:569228/1‑99 (MQ=255)                                                                     gtgtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTaa                        >  2:735143/1‑100 (MQ=255)                                                                      tgtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAAt                       <  2:184041/100‑1 (MQ=255)                                                                       gtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATg                      <  2:845636/100‑1 (MQ=255)                                                                       gtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATg                      >  2:248965/1‑100 (MQ=255)                                                                        tCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGa                     >  2:646627/1‑100 (MQ=255)                                                                         ccGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACTTCCGGCAAAGAGAATACGTGTAATg                      <  2:815714/98‑1 (MQ=255)                                                                          cGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAag                   <  2:834949/100‑1 (MQ=255)                                                                           ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                  <  1:1674315/100‑1 (MQ=255)                                                                              gTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCACCGTCCGGCAAAGAGAATACGTGTAATGAAGAGaa               <  2:672756/100‑1 (MQ=255)                                                                              gTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaa               <  1:394274/100‑1 (MQ=255)                                                                              gTGCGGGCAATAAACAGGGGGTCACTCATGAGATATCCTCCACTTTCCTGCCTGGAGTATAGGCAACGTCCTGCAAAGAGAATACGTGTAATGAaga                  <  1:444600/97‑1 (MQ=255)                                                                               tGCGGGCAATTAACAGAGAATCACTCATGAGATATACTCCAGTTTAATGCCTGGAGTATAGGCAACGTCAGGCAAAGAGAATACGTGTAATGAAGAGaaa              >  2:1517375/1‑100 (MQ=255)                                                                                 cGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAAt            >  2:1763528/1‑100 (MQ=255)                                                                                   ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa          >  1:617266/1‑100 (MQ=255)                                                                                   ggCAATTAACAGGGGGACACTCATGAGATATCCTCCAGTTTCCTACCTGGAGTATAGGCAACGTCCGCCAAAGAGAATACGTGTAATGAAGAGAAAATa           >  2:491543/1‑99 (MQ=255)                                                                                   cgCAATTAACAGGGGGACACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATTAAGAGAAAAt            >  2:867435/2‑98 (MQ=255)                                                                                     cAATTAACAGGGGGTCACTCATGAGATATCCTCCGGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc        >  1:782031/1‑100 (MQ=255)                                                                                     cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc        >  2:2195312/1‑100 (MQ=255)                                                                                      aaTTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCt       <  2:1137304/100‑1 (MQ=255)                                                                                        ttAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa     >  2:1964536/1‑100 (MQ=255)                                                                                         tAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaaa    >  1:1662318/1‑100 (MQ=255)                                                                                         tAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaaa    >  2:254761/1‑100 (MQ=255)                                                                                         tAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAAAGAAAATAACCaaaa    >  1:195979/1‑100 (MQ=255)                                                                                           aCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGc  >  1:1567560/1‑100 (MQ=255)                                                                                            catggggTCACGCATGAGATATCCTCCAGTTTCCTGCATGGAGTATAGGCAACGTCCGGCAAAGCGAATACGTGTAATGCAGCGAAAATATCCTaccg   <  2:2246124/95‑4 (MQ=255)                                                                                                  |                                                                                             TACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGC  >  NC_000913/4489958‑4490147

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