BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,740,514	T→C	R51R (CGT→CGC)	yiaK →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,740,514	0	T	C	99.8%	1763.9 / ‑5.6	540	R51R (CGT→CGC)	yiaK	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (215/324); total (215/325)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 552 total aligned reads displayed.

TGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGC  >  NC_000913/3740424‑3740603

    aGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcatt                                                                                     <  1:8409193/93‑1 (MQ=255)

     gCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcattc                                                                                    >  1:1118022/1‑93 (MQ=255)

      cGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcattca                                                                                   <  1:9859278/93‑1 (MQ=255)

       gAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaa                                                                                  >  2:7064144/1‑93 (MQ=255)

         aaCGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaaca                                                                                >  1:8456323/1‑93 (MQ=255)

          aCGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaacaa                                                                               <  2:9705572/93‑1 (MQ=255)

          aCGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaacaa                                                                               >  2:6736744/1‑93 (MQ=255)

           cGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCaacaac                                                                              <  2:4953232/93‑1 (MQ=255)

            ggCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACt                                                                             <  2:3673776/93‑1 (MQ=255)

            ggCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACt                                                                             >  1:7715725/1‑93 (MQ=255)

             gCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTg                                                                            <  2:2354809/93‑1 (MQ=255)

             gCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTg                                                                            >  2:6535315/1‑93 (MQ=255)

              cTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTgg                                                                           <  1:7869907/93‑1 (MQ=255)

              cTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTgg                                                                           <  1:6220804/93‑1 (MQ=255)

              cTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTgg                                                                           >  1:3882664/1‑93 (MQ=255)

                gACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGaa                                                                         <  2:4609717/93‑1 (MQ=255)

                   gCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAAc                                                                      <  1:3985265/93‑1 (MQ=255)

                      tgtgCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGc                                                                   <  2:4146469/93‑1 (MQ=255)

                      tgtgCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGc                                                                   >  1:5174009/1‑93 (MQ=255)

                        tgCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGa                                                                 <  2:4266635/93‑1 (MQ=255)

                        tgCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGa                                                                 <  2:8195319/93‑1 (MQ=255)

                         gCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGat                                                                >  1:2194590/1‑93 (MQ=255)

                         gCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGat                                                                <  2:7643505/93‑1 (MQ=255)

                         gCAGAGATGTTCGCCCGCACCACCGAACCCGGCGTTTATTCACACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGat                                                                <  2:1969301/93‑1 (MQ=255)

                          cAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGata                                                               >  1:7419142/1‑93 (MQ=255)

                           agagATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGatat                                                              <  2:2760793/93‑1 (MQ=255)

                            gagaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATc                                                             <  1:5056745/93‑1 (MQ=255)

                            gagaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATc                                                             <  2:3086911/93‑1 (MQ=255)

                              gaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCAt                                                           <  2:5034080/93‑1 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         <  2:4803929/93‑1 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         <  1:8624753/93‑1 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         <  1:6931790/93‑1 (MQ=255)

                                  ttCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCt                                                       >  1:5920555/1‑93 (MQ=255)

                                   tCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTg                                                      >  2:9192978/1‑93 (MQ=255)

                                    cGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGa                                                     <  2:5371888/93‑1 (MQ=255)

                                    cGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGa                                                     <  1:4324526/93‑1 (MQ=255)

                                      cccGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATg                                                   >  2:5213622/1‑93 (MQ=255)

                                        cGCACCACCGAATCCGGCGTTTATTCACACGGCGTTAATCGGTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcc                                                 >  2:1099477/1‑93 (MQ=255)

                                         gCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTTCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccc                                                <  2:3469940/93‑1 (MQ=255)

                                         gCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATTccc                                                >  1:5650661/1‑93 (MQ=255)

                                         gCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccc                                                >  2:5821014/1‑93 (MQ=255)

                                          caccacCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccca                                               <  2:2078237/93‑1 (MQ=255)

                                            ccaccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaac                                             >  1:2350159/1‑93 (MQ=255)

                                            ccaccGAATCCGGCGGTTATTCTCACGGCGTTACTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaac                                             <  2:6340998/93‑1 (MQ=255)

                                             caccGAATCCGGCGTTTATTCTCGCGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            <  1:8084554/93‑1 (MQ=255)

                                             caccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            <  2:1410764/93‑1 (MQ=255)

                                               ccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaaccca                                          >  1:8304733/1‑93 (MQ=255)

                                                cGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacccaa                                         <  2:5619924/93‑1 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        <  2:354936/93‑1 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        <  2:3335679/93‑1 (MQ=255)

                                                  aaTCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAAc                                       >  2:1186429/1‑93 (MQ=255)

                                                    tCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGt                                     >  1:445981/1‑93 (MQ=255)

                                                    tCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGt                                     >  2:7421439/1‑93 (MQ=255)

                                                     ccGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGta                                    <  1:6000247/93‑1 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   >  2:4043437/1‑93 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   >  2:3817403/1‑93 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   <  1:7142476/93‑1 (MQ=255)

                                                       ggCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtata                                  >  1:1317611/1‑93 (MQ=255)

                                                        gCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATaa                                 <  2:950236/93‑1 (MQ=255)

                                                        gCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATaa                                 <  2:9442632/93‑1 (MQ=255)

                                                         cGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAAc                                >  2:1684157/1‑93 (MQ=255)

                                                          gTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAAcc                               <  1:4923325/93‑1 (MQ=255)

                                                          gTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAAcc                               >  1:4712200/1‑93 (MQ=255)

                                                           tttATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCa                              >  1:3700183/1‑93 (MQ=255)

                                                           tttATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCa                              >  1:5421245/1‑93 (MQ=255)

                                                           tttATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATAATTCCTGATGCCCAACCCAAACGTATAACCa                              >  2:266560/1‑93 (MQ=255)

                                                            ttATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAg                             >  2:1880074/1‑93 (MQ=255)

                                                            ttATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAg                             <  2:7952202/93‑1 (MQ=255)

                                                            ttATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAg                             >  2:6150644/1‑93 (MQ=255)

                                                             tATTCTCACGTCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGc                            <  2:8690118/93‑1 (MQ=255)

                                                              aTTCTCACGGCGTTAATCGTTTCCCTCGCTTCCTTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGcc                           >  1:4539309/1‑93 (MQ=255)

                                                              aTTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGcc                           >  2:9944137/1‑93 (MQ=255)

                                                              aTTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGcc                           >  1:9402155/1‑93 (MQ=255)

                                                               ttCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCt                          >  1:1703599/1‑93 (MQ=255)

                                                                tctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCACAACCCAAACGTATAACCAGCCTc                         <  2:7239570/93‑1 (MQ=255)

                                                                 ctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCg                        >  2:7809514/1‑93 (MQ=255)

                                                                 ctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCg                        <  1:5277687/93‑1 (MQ=255)

                                                                  tcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCgg                       >  1:1928146/1‑93 (MQ=255)

                                                                  tcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCgg                       <  1:6779774/93‑1 (MQ=255)

                                                                   cACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgc                      <  2:9088512/93‑1 (MQ=255)

                                                                   cACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgc                      <  1:2538937/93‑1 (MQ=255)

                                                                   cACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgc                      <  1:3135993/93‑1 (MQ=255)

                                                                       gCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCaa                  <  1:6524057/93‑1 (MQ=255)

                                                                        cGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAATCGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAt                 <  1:7569481/93‑1 (MQ=255)

                                                                        cGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAt                 <  2:6919336/93‑1 (MQ=255)

                                                                          ttAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTg               >  1:8824232/1‑93 (MQ=255)

                                                                            taTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGaa             <  1:3420860/92‑1 (MQ=255)

                                                                            aaTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGaa             >  1:7289256/1‑93 (MQ=255)

                                                                             aTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAAc            >  1:9106823/1‑93 (MQ=255)

                                                                             aTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAAc            <  1:2094992/93‑1 (MQ=255)

                                                                             aTCGTTTCACTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAAc            <  1:1427395/93‑1 (MQ=255)

                                                                               cGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAg          <  1:4179914/93‑1 (MQ=255)

                                                                               cGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAg          <  1:9551552/93‑1 (MQ=255)

                                                                                 tttCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTg        >  2:589430/1‑93 (MQ=255)

                                                                                    cccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGa     <  2:4518465/93‑1 (MQ=255)

                                                                                    cccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGa     <  1:105301/93‑1 (MQ=255)

                                                                                    cccTCGCTTCATGCAACAACTGGAAAACTGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGa     <  1:27948/93‑1 (MQ=255)

                                                                                     ccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGAc    >  1:4852087/1‑93 (MQ=255)

                                                                                     ccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGAc    <  1:6411728/93‑1 (MQ=255)

                                                                                       tCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGc  <  1:9715608/93‑1 (MQ=255)

TGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGC  >  NC_000913/3740424‑3740603

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