Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	459,256	T→C	100%	D123D (GAT→GAC)	lon →	DNA‑binding ATP‑dependent protease La

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	459,256	0	T	C	100.0%	274.5 / ‑6.7	83	D123D (GAT→GAC)	lon	DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (34/47);  total (35/48)

TGGTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGTCG  >  NC_000913/459159‑459352                                                                                                  |                                                                                                 tGGTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATtgacga                                                                                                >  2:1501880/1‑100 (MQ=255)  ggTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAg                                                                                               >  1:2839637/1‑100 (MQ=255)    tCGAGGTGTTACCGCTCGCGCGTATTTCTCCGCTGTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTCACGAGc                                                                                              >  2:2668622/1‑99 (MQ=255)     cGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCgg                                                                                            <  1:260081/100‑1 (MQ=255)      gAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCggg                                                                                           <  2:2010255/100‑1 (MQ=255)       aGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGa                                                                                          <  1:2919568/100‑1 (MQ=255)        ggggTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGaa                                                                                         >  1:1061269/1‑100 (MQ=255)        cgggTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGCCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCCTTGACGAGc                                                                                              >  2:673834/2‑95 (MQ=255)           ggtACAGCTGGCGTGTATTTATGAGCTATCTGACAATGGCGAACACTTTTCTGCGAAGGCGCAGTATTTGGATTTGCTGACCATAGACGAGCggg                                                                                           <  2:659634/93‑1 (MQ=255)           gTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACAATTGACGAGCGGGAACAg                                                                                      >  2:2393072/1‑100 (MQ=255)            ttACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAgg                                                                                     <  1:868366/100‑1 (MQ=255)            ttACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAgg                                                                                     <  2:1049487/100‑1 (MQ=255)            ttACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAgg                                                                                     >  1:1053417/1‑100 (MQ=255)            ttACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAgg                                                                                     <  1:1754504/100‑1 (MQ=255)               cAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAg                                                                                  <  2:2239608/100‑1 (MQ=255)                 gcgcgcgcGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTg                                                                                >  2:710287/1‑100 (MQ=255)                 gcgcgcgcGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTg                                                                                >  2:1024933/1‑100 (MQ=255)                  cgcgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTg                                                                                >  1:279898/1‑99 (MQ=255)                    cgcgcgTATTTCTGCGCTCTTTTACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATCGACGAGCGGGAACAGGAAGTGCTg                                                                             <  1:2712078/100‑1 (MQ=255)                    cgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCt                                                                              >  1:1293400/1‑99 (MQ=255)                       gcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGt                                                                           <  1:2451873/99‑1 (MQ=255)                          tATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGt                                                                       <  1:2435510/100‑1 (MQ=255)                           aTTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCTTTGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTa                                                                      <  1:46556/100‑1 (MQ=255)                           aTTTCTGCGCTCTCTGACAATGGCGAACACATTTCTGCGAAGGCGCAGTATCTGGAGACGCCGAACATTGACGAGCGGGAACAGGAAAGGCTGGTGCGTa                                                                      >  2:643101/1‑100 (MQ=255)                            tttCTGCGCTATCTGACAATGGCGAACACTTTTCTGCGAATGCGGAGTATCTGGAGTCGCTGACCATTGACGGGCGGGAACAGGAAGTGCTGGTGCGTAc                                                                     <  2:2377335/100‑1 (MQ=255)                              tCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTg                                                                   >  2:2929683/1‑100 (MQ=255)                               cTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGc                                                                  <  2:2117275/100‑1 (MQ=255)                                  cgcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAAt                                                               >  2:1426816/1‑100 (MQ=255)                                  cgcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAAt                                                               >  2:1426897/1‑100 (MQ=255)                                   gcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCaa                                                                <  2:1841001/98‑1 (MQ=255)                                   gcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATc                                                              <  2:2884674/100‑1 (MQ=255)                                   gcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATc                                                              <  2:1747053/100‑1 (MQ=255)                                    ctctctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCa                                                             <  1:2237220/100‑1 (MQ=255)                                    ctctctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGAAAGTGCTGGTGCGTACTGCAATc                                                              >  1:1805930/1‑99 (MQ=255)                                     tctctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAg                                                            <  1:785219/100‑1 (MQ=255)                                      ctctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGc                                                           <  1:858241/100‑1 (MQ=255)                                        ctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCa                                                         >  1:933694/1‑100 (MQ=255)                                         tGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAg                                                        >  2:2359082/1‑100 (MQ=255)                                         tGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAg                                                        <  1:300395/100‑1 (MQ=255)                                         tGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAg                                                        <  2:48078/100‑1 (MQ=255)                                             aaTGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCg                                                    >  2:2643415/1‑100 (MQ=255)                                             aaTGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCg                                                    >  2:564381/1‑100 (MQ=255)                                               tGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGaa                                                  <  1:1195315/100‑1 (MQ=255)                                                ggCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAg                                                 <  1:1400551/100‑1 (MQ=255)                                                ggCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAg                                                 >  2:1460272/1‑100 (MQ=255)                                                    aaCACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTa                                             <  1:128863/100‑1 (MQ=255)                                                    aaCACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTa                                             <  2:2488084/100‑1 (MQ=255)                                                    aaCACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTa                                             <  2:801855/100‑1 (MQ=255)                                                    aaCACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTa                                             <  2:1007516/100‑1 (MQ=255)                                                     acacTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCTTTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTa                                             <  1:761181/99‑1 (MQ=255)                                                      cacTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACa                                           <  1:361350/100‑1 (MQ=255)                                                      cacTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACa                                           <  2:350267/100‑1 (MQ=255)                                                       acTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACAt                                          <  1:625251/100‑1 (MQ=255)                                                        cTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGATCGAAGGCTACATc                                         >  1:763824/1‑100 (MQ=255)                                                        cTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTACTGGAACGAACTGCAATCAGCCAGTTCGAAGGCTACATc                                         >  1:110240/1‑100 (MQ=255)                                                        cTTTTCTGCGAAGGCGGAGAATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATc                                         >  2:469977/1‑100 (MQ=255)                                                        cTTTTCTGCGAAGGCAGAGAATCTAGAGTCGCCGACCCTTGAAGAGCGGGAACAGAAAGTGCTAGTGAGTAGTGACATCAGCGAATTCGAAAGCTAc                                            >  2:631494/1‑97 (MQ=255)                                                            tCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGc                                     <  1:806548/100‑1 (MQ=255)                                                             cTGCGAGGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGGGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCt                                    >  2:2120266/1‑100 (MQ=255)                                                             cTGCGAAGGCGGAGTATCTGGAGTCGTCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCt                                    <  2:151316/100‑1 (MQ=255)                                                             cTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCt                                    >  2:1548409/1‑100 (MQ=255)                                                             cTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCt                                    <  2:684356/100‑1 (MQ=255)                                                                cGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGaa                                 <  1:1447090/100‑1 (MQ=255)                                                                cGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGaa                                 >  2:97972/1‑100 (MQ=255)                                                                 gAAGGCGGAGTATCTGGAGTCGCTGACCATTGACGGGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAAc                                <  1:1468892/100‑1 (MQ=255)                                                                          gTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATc                       <  2:2441618/100‑1 (MQ=255)                                                                            aTCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATccc                     <  2:110239/100‑1 (MQ=255)                                                                             tCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCcca                    <  2:2929341/100‑1 (MQ=255)                                                                             tCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCcca                    >  2:528843/1‑100 (MQ=255)                                                                              cTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccac                   <  1:710291/100‑1 (MQ=255)                                                                               tGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacc                  <  2:343793/100‑1 (MQ=255)                                                                               tGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacc                  <  1:3012351/100‑1 (MQ=255)                                                                               tGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGATAAAAAAAAAAATCccacc                  <  2:955601/100‑1 (MQ=255)                                                                                ggAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacca                 >  2:2964280/1‑100 (MQ=255)                                                                                  aGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGa               >  2:2281452/1‑100 (MQ=255)                                                                                    tCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAg             <  1:1768735/100‑1 (MQ=255)                                                                                       ccGACCATTGACGAGCGGGAACAGGAATTGCTTGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAg             <  2:627370/97‑1 (MQ=255)                                                                                         gACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTg        <  2:2839612/100‑1 (MQ=255)                                                                                          aCCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTg        <  1:741010/99‑1 (MQ=255)                                                                                           ccATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGAc      >  2:1052244/1‑100 (MQ=255)                                                                                            cATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACg     >  1:1621145/1‑100 (MQ=255)                                                                                             aTTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCACTCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACg     >  1:642158/1‑99 (MQ=255)                                                                                             aTTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGt    >  2:2558394/1‑100 (MQ=255)                                                                                              ttGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCAAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCACACCAGAAGTGCTGACATc   >  2:503845/1‑100 (MQ=255)                                                                                               tgacgaGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGTCg  >  2:723129/1‑100 (MQ=255)                                                                                                  |                                                                                                 TGGTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGTCG  >  NC_000913/459159‑459352

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