Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,573,528	A→G	100%	pseudogene (1115/1503 nt)	yjiT →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,528	0	A	G	100.0%	317.7 / ‑6.7	90	pseudogene (1115/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (35/54);  total (36/54)

AAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGC  >  NC_000913/4573432‑4573622                                                                                                 |                                                                                               agaGCTAAGTGAGAATTGGCAGGCAGAATCCAGCGGGAGGTTTGTCGTTAGTTGCTCATGCTTGAGGGATGATTTTTAGGAGCATTAAACTTGATGGCAg                                                                                             <  2:1976617/98‑1 (MQ=255)  aaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGt                                                                                            >  1:2085965/1‑100 (MQ=255)    gCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTg                                                                                           <  2:1803960/99‑1 (MQ=255)       aaGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                        >  2:1698359/1‑99 (MQ=255)         gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                     <  1:1855904/100‑1 (MQ=255)          tGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGGGCATTAAACTTGATGGCAGTGAAATTGc                                                                                    <  2:1126911/100‑1 (MQ=255)          tGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGc                                                                                    <  1:1003032/100‑1 (MQ=255)             gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCa                                                                                  <  1:782849/99‑1 (MQ=255)             gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCa                                                                                  >  1:134212/1‑99 (MQ=255)             gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                 >  2:1228718/1‑100 (MQ=255)             gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                 >  2:1373467/1‑100 (MQ=255)              aTTTTGCAGGCAAAATTCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGAGGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAtt                                                                                <  1:1755957/100‑1 (MQ=255)                ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTgg                                                                              <  1:32375/100‑1 (MQ=255)                ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTgg                                                                              <  1:2249399/100‑1 (MQ=255)                 tGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGt                                                                             <  1:1069214/100‑1 (MQ=255)                  ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGt                                                                             >  1:1898226/1‑99 (MQ=255)                  ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                            <  1:1081279/100‑1 (MQ=255)                    caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAg                                                                          <  2:672518/100‑1 (MQ=255)                     aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAgg                                                                         <  1:1802362/100‑1 (MQ=255)                      ggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGCGGGATGATTGTTGGGAGCATTAAAATTGATGGCAGTGAAATTGCCATTGGTGAgg                                                                         >  2:21796/1‑99 (MQ=255)                       gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                       <  1:695687/100‑1 (MQ=255)                        caAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGa                                                                           >  2:821400/1‑95 (MQ=255)                        caAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGc                                                                      >  2:2159799/1‑100 (MQ=255)                         aaaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGGAATTGCCATTGGTGAGGTGcn                                                                     <  2:1178906/100‑2 (MQ=255)                         aaaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGcc                                                                     <  2:463057/100‑1 (MQ=255)                         aaaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGcc                                                                     >  1:1104655/1‑100 (MQ=255)                          aaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGCATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCg                                                                    <  2:1578187/100‑1 (MQ=255)                           aaTCCATCGGCAGGTTTTTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGt                                                                   <  1:1568724/100‑1 (MQ=255)                           aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGt                                                                   <  2:1024558/100‑1 (MQ=255)                            aTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGtt                                                                  <  2:2004235/100‑1 (MQ=255)                                aTCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                              <  1:1395667/100‑1 (MQ=255)                                 tCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCt                                                             <  1:2092169/100‑1 (MQ=255)                                   ggCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTc                                                           >  1:2121682/1‑100 (MQ=255)                                   ggCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTc                                                           <  1:286071/100‑1 (MQ=255)                                      aGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATc                                                        >  2:688553/1‑100 (MQ=255)                                       gatctGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCg                                                       <  1:345080/96‑1 (MQ=255)                                        gTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGt                                                      <  2:252035/100‑1 (MQ=255)                                         tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGtt                                                     >  2:159613/1‑100 (MQ=255)                                         tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGtt                                                     >  2:1587520/1‑100 (MQ=255)                                           tGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGa                                                   <  1:2065334/100‑1 (MQ=255)                                           tGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGAGAGCCTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGa                                                   >  1:1579264/1‑100 (MQ=255)                                            gTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCgagatc                                                  >  2:345081/1‑95 (MQ=255)                                              cGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                <  1:389669/100‑1 (MQ=255)                                              cGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                >  2:893491/1‑100 (MQ=255)                                               gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                               <  1:1699559/100‑1 (MQ=255)                                               gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                               >  2:2135944/1‑100 (MQ=255)                                               gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGAAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                               >  2:960737/1‑100 (MQ=255)                                                ttAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAg                                              <  1:721109/100‑1 (MQ=255)                                                  aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                            <  2:191605/100‑1 (MQ=255)                                                  aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                            <  1:1601807/100‑1 (MQ=255)                                                     ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGCTGGCAGTGAAATTGCCATTGGTGAGGTGCGGTTAACCTTGATCGTTGATGCAGAt                                            >  2:869891/1‑97 (MQ=255)                                                     ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCaa                                         >  1:2042866/1‑100 (MQ=255)                                                     ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCaa                                         >  1:1754374/1‑100 (MQ=255)                                                       cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       >  1:8660/1‑100 (MQ=255)                                                       cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       <  1:565444/100‑1 (MQ=255)                                                        tCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATgg                                      <  1:1692069/100‑1 (MQ=255)                                                         cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                     <  2:267009/100‑1 (MQ=255)                                                          gTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCt                                    <  1:2186336/100‑1 (MQ=255)                                                           tGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTg                                   >  1:1260658/1‑100 (MQ=255)                                                            gCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGt                                  <  2:729711/100‑1 (MQ=255)                                                            gCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGt                                  <  2:1760028/100‑1 (MQ=255)                                                            gCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGt                                  >  1:510009/1‑100 (MQ=255)                                                            gCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGt                                  <  1:928754/100‑1 (MQ=255)                                                             cTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                 <  1:1446904/100‑1 (MQ=255)                                                             cTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                 >  1:971794/1‑100 (MQ=255)                                                             cTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                 >  1:939326/1‑100 (MQ=255)                                                               ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                               <  1:321552/100‑1 (MQ=255)                                                                gAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACa                              >  1:1500217/1‑100 (MQ=255)                                                                gAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACa                              <  1:1681961/100‑1 (MQ=255)                                                                 aGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAg                             >  1:858004/1‑100 (MQ=255)                                                                  gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                            <  1:1085038/100‑1 (MQ=255)                                                                  gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                            <  1:364054/100‑1 (MQ=255)                                                                  gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                            <  1:310202/100‑1 (MQ=255)                                                                    ggggatTGTTGGGAGCATTAAACTAGATGGAAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                            >  2:925341/4‑98 (MQ=255)                                                                    gatgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTTATGCAGATCAATGGCTGTTACAGGGa                          <  1:456383/100‑1 (MQ=255)                                                                    gatgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGa                          >  1:768237/1‑100 (MQ=255)                                                                       gatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAg                       <  2:2042852/100‑1 (MQ=255)                                                                        atTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAgg                      >  2:637137/1‑100 (MQ=255)                                                                          tgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGct                    <  1:1852846/100‑1 (MQ=255)                                                                           gttgGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGctt                   >  2:898074/1‑100 (MQ=255)                                                                             tgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttct                 <  2:1535252/100‑1 (MQ=255)                                                                              gggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt                <  1:2009288/100‑1 (MQ=255)                                                                               ggAGCATTAAACTTGATGGCAGTGAAATTGGAATTGGTGAGGTGCCGTTAACATTCGTCGTTGCTGCAGATCAATGGCTGTTACAGGGACAGGcttct                 >  2:1778468/1‑98 (MQ=255)                                                                               ggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTg               <  1:25660/100‑1 (MQ=255)                                                                               ggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTg               <  1:1693435/100‑1 (MQ=255)                                                                                  gCATTAAACTTGATGGCAGTGAAATTGCCATTGTTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg            <  1:414395/100‑1 (MQ=255)                                                                                   cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg            <  1:1728680/99‑1 (MQ=255)                                                                                     ttAAACTTGATGGCAGTGAAATTGCCATCGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtgt         >  2:356544/1‑100 (MQ=255)                                                                                        naCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCg      >  2:265529/2‑100 (MQ=255)                                                                                          cTTGATGGCAGTGAAATTGCCATTGGTGAGGTGTCGTTAATCTTCATCGTTGATGCAGATTAATGGCTGTTACTGGGACAGGCTTCTTGCAGTGTGc       <  2:558449/97‑1 (MQ=255)                                                                                            tgatggCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagc  >  2:1829094/1‑100 (MQ=255)                                                                                            tgatggCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagc  >  2:1924953/1‑100 (MQ=255)                                                                                                 |                                                                                               AAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGC  >  NC_000913/4573432‑4573622

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