Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,327,013	Δ1 bp	coding (162/876 nt)	rluB →	23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,327,013	0	T	.	100.0%	269.7 / NA	65	coding (162/876 nt)	rluB	23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (34/31);  total (34/31)

CGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCG  >  NC_000913/1326880‑1327147                                                                                                                                      |                                                                                                                                       cgcgTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGa                                                                                                                                 >  2:365097/1‑140 (MQ=255) cgcgTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGGGTGAGTGTTGTTGGCAAAATTGCTAAACTTGGAGTTCGTGTTGAAGTTACGCCTGGATTGAAAATCCGTATCGATGG‑CACCTGa                                                                                                                                 >  2:306860/1‑140 (MQ=255)  gcgTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGAt                                                                                                                                <  2:1052457/140‑1 (MQ=255)    gTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGAttt                                                                                                                              <  2:1103368/140‑1 (MQ=255)     tGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTc                                                                                                                             <  1:1075165/140‑1 (MQ=255)     tGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTc                                                                                                                             >  2:633131/1‑140 (MQ=255)           ccACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACg                                                                                                                       <  1:478660/140‑1 (MQ=255)                   ctcGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCg                                                                                                               <  1:257928/140‑1 (MQ=255)                   ctcGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCg                                                                                                               <  2:589685/140‑1 (MQ=255)                         gTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAAc                                                                                                         >  2:466860/1‑140 (MQ=255)                           gAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACaa                                                                                                       >  2:512835/1‑140 (MQ=255)                                        ttattGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTc                                                                                          <  2:1324374/140‑1 (MQ=255)                                        ttattGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTc                                                                                          <  1:278353/140‑1 (MQ=255)                                          attGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTg                                                                                        <  2:114714/140‑1 (MQ=255)                                            tGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGc                                                                                      >  2:85059/1‑140 (MQ=255)                                            tGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGc                                                                                      <  2:1466413/140‑1 (MQ=255)                                            tGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGc                                                                                      <  1:466860/140‑1 (MQ=255)                                             gAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGcc                                                                                     <  1:633131/140‑1 (MQ=255)                                                gCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTAt                                                                                  >  1:184230/1‑140 (MQ=255)                                                   ggTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGcc                                                                                     >  2:1510206/1‑134 (MQ=255)                                                   ggTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGcc                                                                                     <  1:1510206/134‑1 (MQ=255)                                                     tCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACaa                                                                             >  2:276915/1‑140 (MQ=255)                                                     tCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACaa                                                                             >  2:876900/1‑140 (MQ=255)                                                     tCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACaa                                                                             >  1:478073/1‑140 (MQ=255)                                                        tgtgAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAAcc                                                                          >  1:1108221/1‑140 (MQ=255)                                                          tgAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCgg                                                                        <  1:876900/140‑1 (MQ=255)                                                           gAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGa                                                                       >  1:1079025/1‑140 (MQ=255)                                                            aGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGaa                                                                      <  1:365097/140‑1 (MQ=255)                                                              tgtTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAgg                                                                    >  1:797693/1‑140 (MQ=255)                                                              tgtTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAgg                                                                    >  1:282572/1‑140 (MQ=255)                                                              tgtTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATCTGTCGCGTTCTGGCCTATTACAAACCGGAAgg                                                                    >  1:511644/1‑140 (MQ=255)                                                               gtTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGt                                                                   >  2:530067/1‑140 (MQ=255)                                                               gtTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGt                                                                   <  2:1170279/140‑1 (MQ=255)                                                                  gATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAg                                                                <  2:138947/140‑1 (MQ=255)                                                                  gATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAg                                                                >  2:1277420/1‑140 (MQ=255)                                                                  gATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAg                                                                <  1:276915/140‑1 (MQ=255)                                                                  gATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAg                                                                <  2:797693/140‑1 (MQ=255)                                                                  gATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAg                                                                <  2:153119/140‑1 (MQ=255)                                                                   aTGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCgg                                                                        >  2:1461736/1‑131 (MQ=255)                                                                   aTGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCgg                                                                        <  1:1461736/131‑1 (MQ=255)                                                                   aTGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTg                                                                  <  2:752484/137‑1 (MQ=255)                                                                   aTGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTg                                                                  >  1:752484/1‑137 (MQ=255)                                                                   aTGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGt                                                               >  1:282264/1‑140 (MQ=255)                                                                   aTGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGt                                                               >  2:543559/1‑140 (MQ=255)                                                                   aTGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGt                                                               >  2:1361027/1‑140 (MQ=255)                                                                    atgCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGtt                                                              >  1:1429244/3‑140 (MQ=255)                                                                      gCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTtgt                                                            >  1:370259/1‑140 (MQ=255)                                                                      gCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTtgt                                                            >  2:1496250/1‑140 (MQ=255)                                                                      gCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTtgt                                                            <  2:282264/140‑1 (MQ=255)                                                                        aaaaTTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTgtgt                                                          <  2:370259/140‑1 (MQ=255)                                                                        aaaaTTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTgtgt                                                          <  1:530067/140‑1 (MQ=255)                                                                               cTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTa                                                   >  1:833417/1‑140 (MQ=255)                                                                               cTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTa                                                   >  1:851995/1‑140 (MQ=255)                                                                                tAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTaa                                                  <  2:1079025/140‑1 (MQ=255)                                                                                               tgtTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACg                                   >  2:789256/1‑140 (MQ=255)                                                                                                    aaGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCaa                              <  2:282572/140‑1 (MQ=255)                                                                                                    aaGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCaa                              <  1:1361027/140‑1 (MQ=255)                                                                                                    aaGTTACCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCaa                              <  1:1221015/140‑1 (MQ=255)                                                                                                         aCCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGgtg                         >  1:445358/1‑140 (MQ=255)                                                                                                         aCCCCTGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGgtg                         >  2:706923/1‑140 (MQ=255)                                                                                                              tGGATTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGa                    >  1:1005313/1‑140 (MQ=255)                                                                                                                 aTTGAAAATCCGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCg                 <  2:1446071/140‑1 (MQ=255)                                                                                                                          ccGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCaaa        >  2:1172712/1‑140 (MQ=255)                                                                                                                           cGTATCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCaaaa       >  1:322256/1‑140 (MQ=255)                                                                                                                               tCGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACATACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGc   >  2:1349085/1‑140 (MQ=255)                                                                                                                                cGATGG‑CACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCg  <  1:512835/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       CGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAAGCCGGTCGTGTGAGTGTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCG  >  NC_000913/1326880‑1327147

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