BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,208,197	A→C	T360P (ACC→CCC)	b2122 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,208,197	0	A	C	97.5%	389.9 / ‑1.8	123	T360P (ACC→CCC)	b2122	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base C (67/52); total (70/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 138 total aligned reads displayed.

GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCT  >  NC_000913/2208107‑2208285

ggAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATCTTGAGCATTATGGTTTACGCATCTGGTTTccc                                                                                        >  2:2997777/1‑93 (MQ=255)

ggAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTccc                                                                                        >  1:2182259/1‑93 (MQ=255)

ggAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGGTTccc                                                                                        >  2:3284110/1‑93 (MQ=255)

ggAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGCTTACGCATCTGGTTTAcc                                                                                        >  2:3284106/1‑93 (MQ=255)

  agaAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGa                                                                                      >  2:1975685/1‑93 (MQ=255)

   gaAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGAt                                                                                     <  2:1232153/93‑1 (MQ=255)

   gaAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGAt                                                                                     <  2:1650973/93‑1 (MQ=255)

     aGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATcc                                                                                   <  2:3785563/93‑1 (MQ=255)

     aGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATcc                                                                                   <  1:3914731/93‑1 (MQ=255)

       cAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCg                                                                                 <  2:2286721/93‑1 (MQ=255)

        aaaGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTtcccgatcccga                                                                                >  2:4422786/1‑93 (MQ=255)

        aaaGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTtcccgatcccga                                                                                >  2:4303827/1‑93 (MQ=255)

         aaGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGac                                                                               >  1:3560684/1‑93 (MQ=255)

           gCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGacac                                                                             <  1:2259838/93‑1 (MQ=255)

           gCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGacac                                                                             >  1:1642796/1‑93 (MQ=255)

            cTTGATCATCTGCGGCTGGTGTCTTTAGGTATTCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGCTTCCCGATACCGACAcc                                                                            >  1:1863327/1‑93 (MQ=255)

            cTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACAcc                                                                            >  2:2924301/1‑93 (MQ=255)

               gATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCGCGACACCGGc                                                                         >  1:1541294/1‑93 (MQ=255)

               gATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGc                                                                         >  2:1464649/1‑93 (MQ=255)

               gATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGc                                                                         >  1:3443664/1‑93 (MQ=255)

                atcatcTGCGTCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCa                                                                        <  1:605556/93‑1 (MQ=255)

                atcatcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGATCATTATGGCTTACGCATCTGGTTTCCCGATCTCGAAACCGGcg                                                                        >  2:628690/1‑92 (MQ=255)

                 tcatcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAg                                                                       <  2:3293112/93‑1 (MQ=255)

                   atcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTa                                                                     <  1:31357/93‑1 (MQ=255)

                   atcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTa                                                                     >  2:3256042/1‑93 (MQ=255)

                    tcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAt                                                                    >  1:3479538/1‑93 (MQ=255)

                    tcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAt                                                                    >  2:3025961/1‑93 (MQ=255)

                     cTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtt                                                                   <  1:823756/93‑1 (MQ=255)

                       gCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtttt                                                                 >  2:4051306/1‑93 (MQ=255)

                       gCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtttt                                                                 <  1:1533638/93‑1 (MQ=255)

                       gCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtttt                                                                 <  1:3886885/93‑1 (MQ=255)

                        cGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTg                                                                <  2:3155432/93‑1 (MQ=255)

                         ggCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGc                                                               <  1:3131857/93‑1 (MQ=255)

                         ggCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGc                                                               >  1:562088/1‑93 (MQ=255)

                          gCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCa                                                              <  2:1933171/93‑1 (MQ=255)

                           cTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTTCCGATCCCGACACCGGGAGttttttttac                                                             >  2:4255524/1‑89 (MQ=255)

                              gtgtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCtt                                                          <  1:426433/93‑1 (MQ=255)

                              gtgtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCtt                                                          >  2:158647/1‑93 (MQ=255)

                               tgtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCttt                                                         >  2:2956821/1‑93 (MQ=255)

                                gtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTc                                                        <  1:2168817/93‑1 (MQ=255)

                                gtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTc                                                        >  2:3207921/1‑93 (MQ=255)

                                 tCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcg                                                       <  1:4010567/93‑1 (MQ=255)

                                   tttAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgcg                                                     >  1:4153441/1‑93 (MQ=255)

                                   tttAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgcg                                                     >  2:4193242/1‑93 (MQ=255)

                                   tttAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgcg                                                     >  1:2545722/1‑93 (MQ=255)

                                   tttAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgcg                                                     >  2:2370668/1‑93 (MQ=255)

                                   tttAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgcg                                                     >  1:503811/1‑93 (MQ=255)

                                       ggTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGt                                                 >  1:2675504/1‑93 (MQ=255)

                                       ggTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGt                                                 >  2:3494099/1‑93 (MQ=255)

                                       ggTAAGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGt                                                 >  1:2369379/1‑93 (MQ=255)

                                        gTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGtt                                                <  1:289166/93‑1 (MQ=255)

                                        gTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGtt                                                <  2:3488374/93‑1 (MQ=255)

                                         tATGCGTTGTTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTg                                               >  2:3143576/1‑93 (MQ=255)

                                         tATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTg                                               <  1:3766892/93‑1 (MQ=255)

                                          aTGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTgg                                              >  1:1756485/1‑93 (MQ=255)

                                            gCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGcc                                            <  2:3453243/93‑1 (MQ=255)

                                              gTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCcgc                                          >  1:4156408/1‑93 (MQ=255)

                                                 gCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGaa                                       <  2:595480/93‑1 (MQ=255)

                                                   tGGCAGGATATTGAGCATTATGGTGGACGCATCTGGTGTCCCGAACCCGACTCCGGCAGGATTTTGCATCTTTCGCGCAGTCGTTCGCGAAGt                                     >  2:4226078/1‑93 (MQ=255)

                                                    ggCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTg                                    <  2:4588076/93‑1 (MQ=255)

                                                    ggCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTg                                    >  1:2827124/1‑93 (MQ=255)

                                                     gCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGa                                   >  2:1570964/1‑93 (MQ=255)

                                                       aGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAAc                                 <  1:3035307/93‑1 (MQ=255)

                                                         gATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAg                               <  2:4514116/93‑1 (MQ=255)

                                                           tatTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGa                             <  1:4534620/93‑1 (MQ=255)

                                                           tatTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGa                             >  2:1044258/1‑93 (MQ=255)

                                                            atTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGaa                            >  2:2606635/1‑93 (MQ=255)

                                                                aGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACt                        >  2:1957311/1‑93 (MQ=255)

                                                                  cATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCa                      <  2:3777832/93‑1 (MQ=255)

                                                                   aTTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAc                     <  1:4203667/93‑1 (MQ=255)

                                                                   aTTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAc                     >  2:3684749/1‑93 (MQ=255)

                                                                    ttATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                    >  2:4383110/1‑93 (MQ=255)

                                                                    ttATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                    >  1:1978188/1‑93 (MQ=255)

                                                                    ttATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                    >  2:2969006/1‑93 (MQ=255)

                                                                    atataGTTCACGTATCTGGTTGCCCTATCCCGACACCGGCAGTATTTTGCATCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                    <  2:1713574/92‑1 (MQ=255)

                                                                     tATGGTTTACGCATCTGGCTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAATTGAACAGGAAAACTCACCg                   >  1:4686810/1‑93 (MQ=255)

                                                                       tGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGc                 <  2:3910559/93‑1 (MQ=255)

                                                                         gTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAg               <  1:3416970/93‑1 (MQ=255)

                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGc              >  1:3997107/1‑93 (MQ=255)

                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGc              <  1:4249500/93‑1 (MQ=255)

                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGc              >  2:131779/1‑93 (MQ=255)

                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGc              <  2:3066929/93‑1 (MQ=255)

                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGc              <  1:3855782/93‑1 (MQ=255)

                                                                            tACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTa            <  1:4551575/93‑1 (MQ=255)

                                                                            tACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTa            <  1:1393408/93‑1 (MQ=255)

                                                                             aCGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAc           >  2:451055/1‑93 (MQ=255)

                                                                             aCGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAc           <  2:3429604/93‑1 (MQ=255)

                                                                             aCGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAc           >  2:4269491/1‑93 (MQ=255)

                                                                             aCGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAc           >  1:1167618/1‑93 (MQ=255)

                                                                              cGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcg          <  1:2852175/93‑1 (MQ=255)

                                                                               gCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcgc         >  1:2731433/1‑93 (MQ=255)

                                                                               gCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcgc         >  1:663462/1‑93 (MQ=255)

                                                                                 ctCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGt       <  2:124064/92‑1 (MQ=255)

                                                                                  tCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTc      >  1:2463354/1‑93 (MQ=255)

                                                                                    tGGTTTCTCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGGAGTGACCAGGAAAACTCACCGGCAGCTACGCGTCgg    >  1:1001381/1‑93 (MQ=255)

                                                                                     ggTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGc   <  1:2705709/93‑1 (MQ=255)

                                                                                     ggTTTCCCGATCCCGACACCGGCAGTATTTCGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGc   <  1:601095/93‑1 (MQ=255)

                                                                                      gTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCt  <  2:4269547/93‑1 (MQ=255)

                                                                                      gTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCt  >  2:3497065/1‑93 (MQ=255)

                                                                                      gTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCt  <  1:3870122/93‑1 (MQ=255)

GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCT  >  NC_000913/2208107‑2208285

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