Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,410,406	T→C	S43P (TCT→CCT)	dusB →	tRNA‑dihydrouridine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,410,406	0	T	C	100.0%	246.5 / NA	70	S43P (TCT→CCT)	dusB	tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/34);  total (36/34)

CTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTG  >  NC_000913/3410272‑3410540                                                                                                                                       |                                                                                                                                       cTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGccttc                                                                                                                                    <  2:701167/139‑1 (MQ=255)    aaCTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTaa                                                                                                                                 >  2:427637/1‑139 (MQ=255)    aaCTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTaa                                                                                                                                 >  2:1118040/1‑139 (MQ=255)        cATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                             <  1:1081505/139‑1 (MQ=255)         aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAAccc                                                                                                                              <  2:468774/137‑1 (MQ=255)         aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAAccc                                                                                                                              >  1:468774/1‑137 (MQ=255)         aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                            >  2:1352874/1‑139 (MQ=255)         aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                            >  1:479565/1‑139 (MQ=255)                 cGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTCgg                                                                                                                    >  1:10532/1‑139 (MQ=255)                           cAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGccttc                                                                                                                                    >  1:1032171/1‑113 (MQ=255)                           cAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGccttc                                                                                                                                    <  2:1032171/113‑1 (MQ=255)                            aGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACaa                                                                                                         <  1:1118040/139‑1 (MQ=255)                              cTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAt                                                                                                       >  2:1109925/1‑139 (MQ=255)                              cTCAGAAATCGCCTGATCGCAGCGCCCATGACTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACATTATCCGAGATGATGCTTTCTAACCCACAGGTTTGGGAAAGCGTCAAAt                                                                                                       >  2:1444979/1‑139 (MQ=255)                                cAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                                     <  1:580417/139‑1 (MQ=255)                                cAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                                     <  1:1352874/139‑1 (MQ=255)                                   aaaTCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGt                                                                                                  >  2:1068524/1‑139 (MQ=255)                                     aTCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGttt                                                                                                <  1:1068524/139‑1 (MQ=255)                                        gCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACg                                                                                             >  1:626501/1‑139 (MQ=255)                                        gCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACg                                                                                             >  2:781969/1‑139 (MQ=255)                                          cTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                           <  1:1109925/139‑1 (MQ=255)                                             aTCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACCGGTTTGGGAAAGCGACAAATCTCGTTTACGGATgg                                                                                        >  1:1070503/1‑139 (MQ=255)                                             aTCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATgg                                                                                        >  1:44507/1‑139 (MQ=255)                                             aTCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGCATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGATTGGGAAAGCGACAAATCTCGTTTACGGATgg                                                                                        >  1:544948/1‑139 (MQ=255)                                                  aGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGcac                                                                                   >  1:64466/1‑139 (MQ=255)                                                       ccATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                                     >  2:792216/1‑116 (MQ=255)                                                       ccATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                                     <  1:792216/116‑1 (MQ=255)                                                        cATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATtgat                                                                             <  1:781969/139‑1 (MQ=255)                                                              tggcATTACAGATAAACCTTTTCTGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                       <  1:1444979/139‑1 (MQ=255)                                                              tggcATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                       <  2:1119253/139‑1 (MQ=255)                                                              tggcATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAACCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                       >  1:16189/1‑139 (MQ=255)                                                                 cATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTg                                                                                                                     <  1:260258/90‑1 (MQ=255)                                                                 cATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTg                                                                                                                     >  2:260258/1‑90 (MQ=255)                                                                     acagacagACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGca                                                                                    >  2:187506/1‑119 (MQ=255)                                                                     acagacagACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGca                                                                                    <  1:187506/119‑1 (MQ=255)                                                                     acagacagACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCAGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTCGGGCAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTc                                                                >  2:212259/1‑139 (MQ=255)                                                                        gacagacCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCa                                                             <  2:64466/139‑1 (MQ=255)                                                                          cagacCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                           >  1:6345/1‑107 (MQ=255)                                                                          cagacCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                           <  2:6345/107‑1 (MQ=255)                                                                             acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                           >  1:59422/1‑104 (MQ=255)                                                                             acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                           <  2:59422/104‑1 (MQ=255)                                                                             acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATtga                                                                              >  2:809476/1‑117 (MQ=255)                                                                             acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATtga                                                                              <  1:809476/117‑1 (MQ=255)                                                                             acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                       >  2:463617/1‑124 (MQ=255)                                                                             acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                       <  1:463617/124‑1 (MQ=255)                                                                                ttttCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGAATTCGCACCGTGCaa                                                     <  1:983050/139‑1 (MQ=255)                                                                                      gACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCt                                               <  1:186742/139‑1 (MQ=255)                                                                                           tgtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                          <  2:87671/139‑1 (MQ=255)                                                                                           tgtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                          <  1:32436/139‑1 (MQ=255)                                                                                            gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGc                                         >  1:1404348/1‑139 (MQ=255)                                                                                                    agaTGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGAt                                      <  1:5814/134‑1 (MQ=255)                                                                                                    agaTGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGAt                                      >  2:5814/1‑134 (MQ=255)                                                                                                       tGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCaa                                                     >  2:956294/1‑116 (MQ=255)                                                                                                       tGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaaga                              >  2:80447/1‑139 (MQ=255)                                                                                                       tGGGAGCCGGATTGACAGTAGCCGAGATGATGCCTTCTAACCCACAGGTTTGCGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCaa                                                     <  1:956294/116‑1 (MQ=255)                                                                                                        gggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaagaa                             <  1:427637/139‑1 (MQ=255)                                                                                                         ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGc                                         >  1:958198/1‑126 (MQ=255)                                                                                                         ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGc                                         <  2:958198/126‑1 (MQ=255)                                                                                                          gAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAAt                           <  2:1070503/139‑1 (MQ=255)                                                                                                             ccGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGc                        <  2:1404348/139‑1 (MQ=255)                                                                                                               ggATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAg                      >  1:1277076/1‑139 (MQ=255)                                                                                                                gATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGa                     <  1:165867/139‑1 (MQ=255)                                                                                                                 aTTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGAt                    >  1:1327869/1‑139 (MQ=255)                                                                                                                        gTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGTATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATgcag                <  1:1051031/136‑1 (MQ=255)                                                                                                                        gTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATgcag                >  2:1051031/1‑136 (MQ=255)                                                                                                                           tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                           >  2:833860/1‑58 (MQ=255)                                                                                                                           tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                           <  1:833860/58‑1 (MQ=255)                                                                                                                           tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACg            <  2:345740/137‑1 (MQ=255)                                                                                                                           tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACg            >  1:345740/1‑137 (MQ=255)                                                                                                                           tCCAAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTa          <  1:419760/139‑1 (MQ=255)                                                                                                                              gagaTGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTa       <  2:626501/139‑1 (MQ=255)                                                                                                                                   gatgCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTg  >  1:86193/1‑139 (MQ=255)                                                                                                                                   gatgCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTg  >  2:862726/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       CTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTG  >  NC_000913/3410272‑3410540

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: ‑