BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,692	C→T	G81S (GGC→AGC)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,692	0	C	T	99.5%	750.9 / ‑5.6	207	G81S (GGC→AGC)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (104/102); total (105/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01

Only 100 of 219 total aligned reads displayed.

GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781

gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  1:3977781/1‑93 (MQ=255)

gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  1:5166315/1‑93 (MQ=255)

gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  2:246116/1‑93 (MQ=255)

 ctctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgcc                                                                                        <  1:3335134/93‑1 (MQ=255)

 ctctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgcc                                                                                        >  1:6671973/1‑93 (MQ=255)

  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  1:5699177/1‑93 (MQ=255)

  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       <  2:1767381/93‑1 (MQ=255)

  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       <  2:1436738/93‑1 (MQ=255)

    ttCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGctgcctgc                                                                                     >  1:7720976/1‑93 (MQ=255)

      ccGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCaa                                                                                   <  2:1394798/93‑1 (MQ=255)

        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 >  2:4268801/1‑93 (MQ=255)

        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 <  1:3069862/93‑1 (MQ=255)

        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 >  1:2807616/1‑93 (MQ=255)

          cTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCg                                                                               >  1:4632420/1‑93 (MQ=255)

          cTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCg                                                                               >  1:6459286/1‑93 (MQ=255)

           ttGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGc                                                                              >  2:4813397/1‑93 (MQ=255)

            tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             <  1:4373946/93‑1 (MQ=255)

              cTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCTATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGa                                                                           >  1:2771742/1‑93 (MQ=255)

               ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          <  1:1563817/93‑1 (MQ=255)

                tttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAcc                                                                         <  2:5360658/93‑1 (MQ=255)

                  tCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCgg                                                                       <  1:4029366/93‑1 (MQ=255)

                  tCGCCTCCAGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCgg                                                                       <  1:1282764/93‑1 (MQ=255)

                     ccTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGaa                                                                    <  1:2134104/93‑1 (MQ=255)

                     ccTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGaa                                                                    >  2:711867/1‑93 (MQ=255)

                      cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAAt                                                                   <  2:919301/93‑1 (MQ=255)

                      cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAAt                                                                   <  1:3281386/93‑1 (MQ=255)

                       tCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATc                                                                  <  2:5207028/93‑1 (MQ=255)

                       tCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATc                                                                  >  2:785866/1‑93 (MQ=255)

                        ccGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCg                                                                 >  1:1488628/1‑93 (MQ=255)

                         cGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGa                                                                >  2:1054142/1‑93 (MQ=255)

                          ggTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaa                                                               <  1:1388421/93‑1 (MQ=255)

                           gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              <  2:1595564/93‑1 (MQ=255)

                           gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              <  2:5435246/93‑1 (MQ=255)

                            tGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAAt                                                             <  1:2099273/93‑1 (MQ=255)

                             gAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATg                                                            >  2:6302779/1‑93 (MQ=255)

                              aatcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGc                                                           >  1:5950326/1‑93 (MQ=255)

                               atcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGcc                                                          >  1:7914753/1‑93 (MQ=255)

                               atcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGcc                                                          <  1:2416829/93‑1 (MQ=255)

                                tcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCa                                                         >  1:2150222/1‑93 (MQ=255)

                                  aatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaaa                                                       >  1:412103/1‑93 (MQ=255)

                                  aatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaaa                                                       <  2:6159985/93‑1 (MQ=255)

                                  aatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaaa                                                       <  1:2638562/93‑1 (MQ=255)

                                   atcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAAt                                                      <  2:7471620/93‑1 (MQ=255)

                                    tcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATg                                                     >  2:5469936/1‑93 (MQ=255)

                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  1:4726425/1‑93 (MQ=255)

                                       cTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt                                                  >  1:6070072/1‑93 (MQ=255)

                                       cTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt                                                  >  2:4089077/1‑93 (MQ=255)

                                       cTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt                                                  <  2:3853209/93‑1 (MQ=255)

                                       cTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt                                                  <  2:5897246/93‑1 (MQ=255)

                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 >  1:6249510/1‑93 (MQ=255)

                                         ttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATcc                                                >  1:7439879/1‑93 (MQ=255)

                                           cTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCAt                                              <  2:1879756/93‑1 (MQ=255)

                                            tGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATc                                             <  2:2284033/93‑1 (MQ=255)

                                              ttGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCAc                                           <  2:2329352/93‑1 (MQ=255)

                                              ttGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCAc                                           <  1:2248011/93‑1 (MQ=255)

                                                 cGGAAAGAGCATCCGACAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAg                                        <  1:7000036/93‑1 (MQ=255)

                                                  ggAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGt                                       >  2:2266711/1‑93 (MQ=255)

                                                  ggAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGt                                       <  2:6755792/93‑1 (MQ=255)

                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGtt                                      >  2:3042268/1‑93 (MQ=255)

                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGtt                                      <  1:5780995/93‑1 (MQ=255)

                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGtt                                      >  1:7065662/1‑93 (MQ=255)

                                                      agagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGcc                                   >  1:6316285/1‑93 (MQ=255)

                                                       gagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGccc                                  >  1:1541171/1‑93 (MQ=255)

                                                       gagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGccc                                  >  2:1724827/1‑93 (MQ=255)

                                                       gagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGccc                                  <  2:6239853/93‑1 (MQ=255)

                                                         gCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCgg                                >  2:6886214/1‑93 (MQ=255)

                                                            tCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGt                             >  1:5679370/1‑93 (MQ=255)

                                                            tCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGt                             >  2:2579907/1‑93 (MQ=255)

                                                            tCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGt                             >  1:7206231/1‑93 (MQ=255)

                                                             ccGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTa                            >  1:529524/1‑93 (MQ=255)

                                                             ccGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTa                            >  2:2792059/1‑93 (MQ=255)

                                                               gAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAAt                          >  2:3717581/1‑93 (MQ=255)

                                                               gAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAAt                          >  2:3546477/1‑93 (MQ=255)

                                                                  aaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       >  2:1281964/1‑93 (MQ=255)

                                                                   aCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGc                      <  1:345950/93‑1 (MQ=255)

                                                                     cGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGcca                    >  2:1103184/1‑93 (MQ=255)

                                                                      gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   <  1:1697998/93‑1 (MQ=255)

                                                                       aTAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccacc                  >  1:5046730/1‑93 (MQ=255)

                                                                        tAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccacca                 <  1:3742424/93‑1 (MQ=255)

                                                                           aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              <  1:1048811/93‑1 (MQ=255)

                                                                           aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              <  2:7720975/93‑1 (MQ=255)

                                                                             gcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatg            >  2:4413208/1‑93 (MQ=255)

                                                                             gcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatg            >  2:4041384/1‑93 (MQ=255)

                                                                              cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           <  2:842784/93‑1 (MQ=255)

                                                                              cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           <  1:336920/93‑1 (MQ=255)

                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTCATACGCCACCATCatgat          >  1:5967486/1‑93 (MQ=255)

                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          >  1:4662006/1‑93 (MQ=255)

                                                                                acgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAGTGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgatg         <  2:7258551/93‑1 (MQ=255)

                                                                                 cgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATgg        <  2:6599283/93‑1 (MQ=255)

                                                                                  gaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGc       >  2:5074546/1‑93 (MQ=255)

                                                                                   aTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCa      <  2:1206345/93‑1 (MQ=255)

                                                                                   aTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCa      >  2:2419142/1‑93 (MQ=255)

                                                                                    ttAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaa     <  2:794359/93‑1 (MQ=255)

                                                                                    ttAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaa     <  2:5281403/93‑1 (MQ=255)

                                                                                     tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  1:2958348/93‑1 (MQ=255)

                                                                                     tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  1:5579828/93‑1 (MQ=255)

                                                                                     tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  1:6894625/93‑1 (MQ=255)

                                                                                       aGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATa  >  2:634302/1‑93 (MQ=255)

                                                                                       aGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATa  <  1:1970330/93‑1 (MQ=255)

                                                                                       aGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATa  <  2:4338810/93‑1 (MQ=255)

GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`