BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,829,419	C→T	R446C (CGC→TGC)	yicH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,829,419	0	C	T	99.6%	831.8 / NA	232	R446C (CGC→TGC)	yicH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (87/144); minor base G (0/1); total (87/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

ACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC  >  NC_000913/3829329‑3829508

aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGc                                                                                         <  1:2978015/93‑1 (MQ=255)

   cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGGTTTTCGTTGGCCg                                                                                      >  2:2425151/1‑93 (MQ=255)

   cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCg                                                                                      >  2:4625693/1‑93 (MQ=255)

    tGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTTGTTGCCCGt                                                                                     >  2:3234695/1‑93 (MQ=255)

    tGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGt                                                                                     <  2:3035594/93‑1 (MQ=255)

       ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTTGTTGCCCGTCGc                                                                                  >  2:307767/1‑93 (MQ=255)

       ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGc                                                                                  >  1:128100/1‑93 (MQ=255)

       ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGGTTATGTTCGTTGCCCGTCGc                                                                                  >  2:1168280/1‑93 (MQ=255)

        tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTTATGTTCGTTGCCCGTCGCt                                                                                 <  1:1424503/93‑1 (MQ=255)

         aCCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTg                                                                                <  1:4254507/93‑1 (MQ=255)

          ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTgg                                                                               <  1:4326035/93‑1 (MQ=255)

            gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTTGTTGGCCGTCGCTGgcg                                                                             >  2:3225449/1‑93 (MQ=255)

            gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTCCCCGTCGCTGgcg                                                                             >  1:3006560/1‑93 (MQ=255)

             aTCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGgcgc                                                                            <  2:4317942/93‑1 (MQ=255)

             aTCATAAATGGGGCGTATGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGgcgc                                                                            <  1:2967117/93‑1 (MQ=255)

              tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCt                                                                           >  2:3396141/1‑93 (MQ=255)

                 tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGAc                                                                        >  1:2223131/1‑93 (MQ=255)

                 tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGAc                                                                        >  2:1755399/1‑93 (MQ=255)

                 aaaaGGGGGCGTCGGGGGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGAc                                                                        <  2:2928525/92‑1 (MQ=255)

                   aaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGAccc                                                                      >  1:3621144/1‑92 (MQ=255)

                     tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGAACGcc                                                                    >  1:1956472/1‑93 (MQ=255)

                     tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTATTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGcc                                                                    <  1:4230226/93‑1 (MQ=255)

                       gggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCaa                                                                  <  2:1468550/93‑1 (MQ=255)

                       gggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCaa                                                                  <  2:775074/93‑1 (MQ=255)

                        ggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAc                                                                 <  1:2125935/93‑1 (MQ=255)

                         gCGTCTGGAGTGGCTCGGCGAATCTGAGTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAca                                                                <  1:544079/93‑1 (MQ=255)

                         gCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGGCAAca                                                                >  1:950497/1‑93 (MQ=255)

                         gCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAca                                                                <  2:3571617/93‑1 (MQ=255)

                         gCGTCTGGAGTGGCTCGGCAAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAca                                                                <  2:3512861/93‑1 (MQ=255)

                           gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                              <  1:1675356/93‑1 (MQ=255)

                           gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                              >  1:3329696/1‑93 (MQ=255)

                           gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                              >  1:2567169/1‑93 (MQ=255)

                            tCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCCAcagca                                                             >  2:2512630/1‑93 (MQ=255)

                              tGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACagcagc                                                           >  2:1068238/1‑93 (MQ=255)

                               ggAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTTACCGCCAACAgcagca                                                          >  2:2769725/1‑93 (MQ=255)

                               ggAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAgcagca                                                          >  2:2890538/1‑93 (MQ=255)

                                gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGCCCGCCAACAGCAGCAc                                                         <  2:711061/93‑1 (MQ=255)

                                gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCAc                                                         <  1:1025930/93‑1 (MQ=255)

                                  gTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                                       <  2:1234738/93‑1 (MQ=255)

                                    ggCTCGCCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                                                     <  2:2317549/93‑1 (MQ=255)

                                       tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                                                  >  2:1082591/1‑93 (MQ=255)

                                        cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAata                                                 >  2:3987377/1‑93 (MQ=255)

                                         ggCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAatat                                                <  2:475422/93‑1 (MQ=255)

                                          gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAatata                                               >  2:2011983/1‑92 (MQ=255)

                                          gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                                               <  2:2670328/93‑1 (MQ=255)

                                          gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                                               <  1:2835545/93‑1 (MQ=255)

                                           cGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                                              >  2:1212896/1‑93 (MQ=255)

                                            gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                                             >  1:3208921/1‑93 (MQ=255)

                                            gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                                             <  1:1333353/93‑1 (MQ=255)

                                             aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                                            <  2:3758560/93‑1 (MQ=255)

                                             aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                                            <  1:1040538/93‑1 (MQ=255)

                                             aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                                            >  1:2080464/1‑93 (MQ=255)

                                             aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                                            >  1:2357116/1‑93 (MQ=255)

                                              aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCg                                           >  2:3549140/1‑93 (MQ=255)

                                                cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGaa                                         <  1:4182739/93‑1 (MQ=255)

                                                 tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                        <  1:649667/93‑1 (MQ=255)

                                                  gAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACt                                       <  1:944209/93‑1 (MQ=255)

                                                    aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                                     >  2:602751/1‑93 (MQ=255)

                                                    aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                                     >  2:163165/1‑93 (MQ=255)

                                                      gccgccgccGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                   <  1:1853690/93‑1 (MQ=255)

                                                      gccgccgccGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                   <  1:2938267/93‑1 (MQ=255)

                                                      gccgccgccGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                   <  1:459896/93‑1 (MQ=255)

                                                        cgccgccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc                                 <  1:4064599/93‑1 (MQ=255)

                                                        cgccgccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc                                 <  1:2485025/93‑1 (MQ=255)

                                                         gccgccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                <  1:1094427/93‑1 (MQ=255)

                                                           cgccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAtt                              <  2:1867503/93‑1 (MQ=255)

                                                             ccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                            <  1:3383030/93‑1 (MQ=255)

                                                             ccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                            <  2:2626007/93‑1 (MQ=255)

                                                             ccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                            <  1:4533769/93‑1 (MQ=255)

                                                             ccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                            <  2:3718216/93‑1 (MQ=255)

                                                             ccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                            >  2:1289253/1‑93 (MQ=255)

                                                              cgcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                           <  1:2695361/93‑1 (MQ=255)

                                                              cgcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                           <  2:3360385/93‑1 (MQ=255)

                                                               gcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                          <  2:2500235/93‑1 (MQ=255)

                                                               gcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                          >  1:3943636/1‑93 (MQ=255)

                                                               gcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                          >  1:3026203/1‑93 (MQ=255)

                                                                cgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTg                         >  1:760611/1‑93 (MQ=255)

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

                                                                    aTTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaa                     <  1:1630058/93‑1 (MQ=255)

                                                                    aTTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaa                     >  1:3742060/1‑93 (MQ=255)

                                                                     ttCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaaa                    >  1:2321920/1‑93 (MQ=255)

                                                                     ttCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaaa                    <  1:850892/93‑1 (MQ=255)

                                                                     ttCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaaa                    <  2:1600199/93‑1 (MQ=255)

                                                                      tCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGCACTGAGTGCATTTACTGAAAAAg                   <  1:402580/93‑1 (MQ=255)

                                                                      tCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAg                   <  1:1157292/93‑1 (MQ=255)

                                                                      tCAATCGTGTTGATGTTCGTTGCCCGTCGCTAGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAg                   <  1:2524531/93‑1 (MQ=255)

                                                                        aaTCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGc                 <  2:1938573/93‑1 (MQ=255)

                                                                         aTCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCa                <  2:2220915/93‑1 (MQ=255)

                                                                          tCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAt               >  1:3900723/1‑93 (MQ=255)

                                                                           cGTGTTGATGTTCGTTGCCCGTCGCTGGTGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAtt              <  2:3854179/93‑1 (MQ=255)

                                                                            gtgtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTc             <  2:2956589/93‑1 (MQ=255)

                                                                            gtgtTGATGTTCGTTGCCCGGCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAtgc             <  2:2988201/93‑3 (MQ=255)

                                                                             tgtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCt            <  1:1262568/93‑1 (MQ=255)

                                                                              gtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTg           >  1:4381762/1‑93 (MQ=255)

                                                                               ttGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTgg          >  2:93614/1‑93 (MQ=255)

                                                                               ttGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTgg          >  1:1797309/1‑93 (MQ=255)

                                                                                 gATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGaa        <  2:2139057/93‑1 (MQ=255)

                                                                                   tGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGc      <  2:2090118/93‑1 (MQ=255)

                                                                                   tGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGc      >  1:891375/1‑93 (MQ=255)

                                                                                      tcgttgCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCAc   <  2:4430886/93‑1 (MQ=255)

ACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC  >  NC_000913/3829329‑3829508

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