BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,321,028	A→G	R228R (CGT→CGC)	mrsA ←	MrsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,321,028	0	A	G	99.7%	1073.6 / NA	345	R228R (CGT→CGC)	mrsA	MrsA
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (179/165); total (179/166)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 378 total aligned reads displayed.

ATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGC  >  NC_000913/3320938‑3321117

aTCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGTAATACCGAGATCCGCTTTTTCAGCCAGCAcgcg                                                                                         >  1:7627661/1‑93 (MQ=255)

aTCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCAcgcg                                                                                         <  2:1254810/93‑1 (MQ=255)

aTCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCAcgcg                                                                                         >  1:4799313/1‑93 (MQ=255)

aTCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCAcgcg                                                                                         >  1:4441773/1‑93 (MQ=255)

aTCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCAcgcg                                                                                         >  2:4937671/1‑93 (MQ=255)

 tCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGa                                                                                        <  1:6169282/93‑1 (MQ=255)

 tCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGa                                                                                        <  2:3910360/93‑1 (MQ=255)

   gACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGc                                                                                      <  1:7181261/93‑1 (MQ=255)

   gACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATACGCTTTTTCAGCCAGCACGCGAGc                                                                                      >  1:3313975/1‑93 (MQ=255)

    aCTTTATTGCCTTGATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGcc                                                                                     <  1:5801218/93‑1 (MQ=255)

    aCTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGcc                                                                                     <  2:5438459/93‑1 (MQ=255)

    aCTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGcc                                                                                     >  2:1371601/1‑93 (MQ=255)

     cTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCt                                                                                    >  2:7366278/1‑93 (MQ=255)

      tttATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTg                                                                                   >  2:2360242/1‑93 (MQ=255)

      tttATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTg                                                                                   >  1:7730239/1‑93 (MQ=255)

         aTTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAg                                                                                >  1:326389/1‑93 (MQ=255)

         aTTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAg                                                                                >  1:1607599/1‑93 (MQ=255)

              cTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAgcgcgc                                                                           >  1:6859236/1‑93 (MQ=255)

               ttCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGcgcgcg                                                                          >  1:5880903/1‑93 (MQ=255)

               ttCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGcgcgcg                                                                          >  2:7579192/1‑93 (MQ=255)

                 cATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGaa                                                                        <  2:2314443/93‑1 (MQ=255)

                  aTGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGCAc                                                                       >  2:1062115/1‑93 (MQ=255)

                   tGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACg                                                                      <  2:5276156/93‑1 (MQ=255)

                   tGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTCGAGCGCGCGAACg                                                                      >  1:368374/1‑93 (MQ=255)

                     gTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTc                                                                    <  1:4227057/93‑1 (MQ=255)

                     gTCAACCATAAGCACGAGATCGCCATCGCCGTGGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTc                                                                    <  1:121193/93‑1 (MQ=255)

                      tCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCg                                                                   >  1:2876046/1‑93 (MQ=255)

                      tCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCg                                                                   <  2:1946094/93‑1 (MQ=255)

                      tCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTGTTTCAGCCAGCACGCGAGCCTGGAGGGCGCGAACGGCg                                                                   >  1:82629/1‑93 (MQ=255)

                        aaCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCggg                                                                 >  2:6957530/1‑92 (MQ=255)

                        aaCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGt                                                                 >  2:1329183/1‑93 (MQ=255)

                         aCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTa                                                                <  2:949669/93‑1 (MQ=255)

                         aCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTa                                                                >  1:2149463/1‑93 (MQ=255)

                         aCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTa                                                                >  1:4130763/1‑93 (MQ=255)

                          ccATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAg                                                               >  2:1027373/1‑93 (MQ=255)

                            aTAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGcc                                                             <  1:3931956/93‑1 (MQ=255)

                            aTAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGcc                                                             <  1:1675133/93‑1 (MQ=255)

                             tAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGccc                                                            >  1:4977671/1‑93 (MQ=255)

                             tAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGccc                                                            >  2:2242605/1‑93 (MQ=255)

                             tAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGccc                                                            >  1:1331991/1‑93 (MQ=255)

                              aaTCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGcccc                                                           <  2:6347214/93‑1 (MQ=255)

                                 cACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCAc                                                        <  1:2357541/93‑1 (MQ=255)

                                 cACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCAc                                                        >  2:7462950/1‑93 (MQ=255)

                                  ccgcgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACt                                                       <  2:2863640/92‑1 (MQ=255)

                                  aCGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACt                                                       >  2:2118798/1‑93 (MQ=255)

                                   cgcgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACtt                                                      >  2:1668404/1‑93 (MQ=255)

                                   cgcgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACtc                                                      >  1:190691/1‑92 (MQ=255)

                                    gcgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTc                                                     <  1:6636518/93‑1 (MQ=255)

                                     cgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCg                                                    <  1:7533908/93‑1 (MQ=255)

                                       atcgccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGc                                                  >  1:1987487/1‑93 (MQ=255)

                                        tcgccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCGCGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCa                                                 <  2:81847/93‑1 (MQ=255)

                                           ccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTg                                              <  1:2567515/93‑1 (MQ=255)

                                           ccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCAAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTg                                              <  2:1426163/93‑1 (MQ=255)

                                            catcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGa                                             <  2:4490979/93‑1 (MQ=255)

                                             atcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGAt                                            >  2:2906892/1‑93 (MQ=255)

                                             atcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGAt                                            <  2:1501163/93‑1 (MQ=255)

                                             atcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGAt                                            >  2:5757539/1‑93 (MQ=255)

                                              tcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATg                                           >  1:4156458/1‑93 (MQ=255)

                                                  cGTCGAAGGCAATACCGCGAGCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTa                                       <  2:3379966/93‑1 (MQ=255)

                                                   gTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTac                                      <  2:6469404/93‑1 (MQ=255)

                                                   gTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTac                                      >  2:2198680/1‑93 (MQ=255)

                                                   gTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGGAGCCCCCACTTCGGCATTGATGTTTac                                      >  1:4311788/1‑93 (MQ=255)

                                                    tCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTaca                                     <  2:2580480/93‑1 (MQ=255)

                                                     cGAAGACAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTacac                                    <  2:3853630/93‑1 (MQ=255)

                                                      gAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACAcc                                   <  1:2977444/93‑1 (MQ=255)

                                                      gAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACAcc                                   <  1:4730146/93‑1 (MQ=255)

                                                      gAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACAcc                                   <  1:568402/93‑1 (MQ=255)

                                                       aaGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCg                                  >  2:5977850/1‑93 (MQ=255)

                                                       aaGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCg                                  >  2:7303059/1‑93 (MQ=255)

                                                        aGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGt                                 >  2:4458036/1‑93 (MQ=255)

                                                         ggCAATGCCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGtt                                <  1:1774366/93‑1 (MQ=255)

                                                         ggCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGtt                                <  1:5324102/93‑1 (MQ=255)

                                                         ggCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGtt                                <  2:667998/93‑1 (MQ=255)

                                                          gCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGttt                               <  2:160192/93‑1 (MQ=255)

                                                          gCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGttt                               <  2:62828/93‑1 (MQ=255)

                                                           cAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTg                              <  2:3064937/93‑1 (MQ=255)

                                                           cAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTg                              <  2:7793601/93‑1 (MQ=255)

                                                           cAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTg                              <  2:4114351/93‑1 (MQ=255)

                                                           cAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTg                              <  1:3802912/93‑1 (MQ=255)

                                                             aTACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGc                            <  1:7274615/93‑1 (MQ=255)

                                                               aCCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTc                          >  2:7060019/1‑93 (MQ=255)

                                                                ccGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTca                         <  1:6396830/93‑1 (MQ=255)

                                                                  gagaTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTcaca                       >  2:2503950/1‑93 (MQ=255)

                                                                    gaTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAAc                     >  2:4338476/1‑93 (MQ=255)

                                                                      tCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCg                   <  1:1042120/93‑1 (MQ=255)

                                                                        cGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGAt                 >  1:3551826/1‑93 (MQ=255)

                                                                         gCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATa                >  1:5966260/1‑93 (MQ=255)

                                                                         gCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATa                >  1:4385480/1‑93 (MQ=255)

                                                                         gCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATa                <  2:1727137/93‑1 (MQ=255)

                                                                          cTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAg               <  2:5541044/93‑1 (MQ=255)

                                                                           tttttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGc              >  1:3154077/1‑93 (MQ=255)

                                                                           tttttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGc              <  1:2931704/93‑1 (MQ=255)

                                                                              ttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGAt           <  2:3233515/93‑1 (MQ=255)

                                                                               tCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATa          >  1:3597263/1‑93 (MQ=255)

                                                                               tCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATa          >  2:3788862/1‑93 (MQ=255)

                                                                                cagccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATaa         <  2:6114114/93‑1 (MQ=255)

                                                                                 agccagtACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAAc        <  1:2190374/93‑1 (MQ=255)

                                                                                  gccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACg       >  1:2702671/1‑93 (MQ=255)

                                                                                   ccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGt      >  2:6757120/1‑93 (MQ=255)

                                                                                       cacgCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTAGCTCACAACCGATAGCGATAACGTTCGc  >  1:1096151/1‑93 (MQ=255)

ATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGC  >  NC_000913/3320938‑3321117

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