BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,472,351	A→G	H36R (CAC→CGC)	ydbC →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,472,351	0	A	G	98.5%	400.1 / ‑3.5	129	H36R (CAC→CGC)	ydbC	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (46/81); total (47/82)
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.32e-01

Only 100 of 140 total aligned reads displayed.

TACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGG‑CCC‑CCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATC  >  NC_000913/1472262‑1472440

tACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGG‑CCC‑CCCACGAGATCGCCGCGt                                                                                        >  1:4344238/1‑93 (MQ=255)

 aCTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGATTTTGGG‑CCC‑CCCACGAGATCGCcgcgc                                                                                        >  2:3000554/1‑92 (MQ=255)

 aCTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCcgcgc                                                                                        >  2:362395/1‑92 (MQ=255)

  ctctCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTc                                                                                       >  2:4397904/1‑93 (MQ=255)

  ctctCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTc                                                                                       <  1:4312180/93‑1 (MQ=255)

   tctcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCg                                                                                      >  1:2156770/1‑93 (MQ=255)

    ctcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGc                                                                                     <  1:4113670/93‑1 (MQ=255)

     tcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCt                                                                                    >  1:1780677/1‑93 (MQ=255)

     tcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCt                                                                                    >  1:2715724/1‑93 (MQ=255)

     tcGGTACAAAATCCGTTAACCGTCTTGGGTATGGCGCGCTGCAACTGGCAGGTCCTCGAG‑TTTTTGG‑CCCTCCGACGAGATCGCCGCGTCGcc                                                                                    >  1:3091983/1‑92 (MQ=255)

       ggTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCtat                                                                                  >  2:3122593/1‑93 (MQ=255)

         tACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATaa                                                                                >  1:1279987/1‑93 (MQ=255)

         tACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATaa                                                                                <  1:4414416/93‑1 (MQ=255)

          aCAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAAc                                                                               <  1:4337309/93‑1 (MQ=255)

            aaaaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCCCGAGATCGCCGGTTCGCTATACCCg                                                                             >  1:2901386/1‑93 (MQ=255)

            aaaaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCg                                                                             <  1:3053367/93‑1 (MQ=255)

               aTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGtgc                                                                          <  1:1310450/93‑1 (MQ=255)

               aTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGtgc                                                                          <  2:1493436/93‑1 (MQ=255)

               aTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGtgc                                                                          <  1:4680344/93‑1 (MQ=255)

               aTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGtgc                                                                          <  1:3172507/93‑1 (MQ=255)

                tCCGTTATCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCTCTATAACCGtgct                                                                         <  2:1123744/93‑1 (MQ=255)

                tCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGtgct                                                                         <  2:3449712/93‑1 (MQ=255)

                accGTTAACCGTCTTGGTTATGGCGCGATGCAGCTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGtgct                                                                         <  2:3414910/92‑1 (MQ=255)

                 ccGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑GTTTTGGTCCC‑CCCACGAGATCGCCGCGTCGCTATAAACGtgctg                                                                        >  1:1909429/1‑93 (MQ=255)

                   gTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCg                                                                      >  1:3302301/1‑93 (MQ=255)

                    ttAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGt                                                                     <  2:717850/93‑1 (MQ=255)

                    ttAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGt                                                                     <  1:3301521/93‑1 (MQ=255)

                    gtAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGt                                                                     <  1:4188325/92‑1 (MQ=255)

                      aaCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGa                                                                   >  1:778731/1‑93 (MQ=255)

                      aaCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGa                                                                   <  1:1418390/93‑1 (MQ=255)

                       aCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAg                                                                  <  2:4012596/93‑1 (MQ=255)

                        ccGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAgg                                                                 <  2:2151705/93‑1 (MQ=255)

                          gTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgcg                                                               >  1:2290815/1‑93 (MQ=255)

                          gTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgcg                                                               >  1:3616653/1‑93 (MQ=255)

                           tCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgcgc                                                              >  1:649494/1‑93 (MQ=255)

                           tCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgcgc                                                              <  1:865624/93‑1 (MQ=255)

                           tCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgcgc                                                              <  1:3946531/93‑1 (MQ=255)

                           tCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATAGCCGCGTCGCTATAACCGTGCTGCGTGAGgcga                                                              >  1:62395/1‑92 (MQ=255)

                            cTTGGTTATGGCGCGATGCAACTGGTAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATTGCCGCGTCGCTATAACCGTGCGGTTTGAGGGGCt                                                             >  1:3539244/1‑93 (MQ=255)

                            cTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCt                                                             <  2:1773443/93‑1 (MQ=255)

                            cTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCt                                                             <  2:1395909/93‑1 (MQ=255)

                            cTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCt                                                             <  2:4212788/93‑1 (MQ=255)

                              tGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTgg                                                           >  1:1863423/1‑93 (MQ=255)

                                gTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGACGCTATAACCGTGCTGCGTGAGGCGCTGGCa                                                         >  2:1052382/1‑93 (MQ=255)

                                   aTGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTg                                                      >  2:4508171/1‑93 (MQ=255)

                                    tGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTgg                                                     <  1:3571964/93‑1 (MQ=255)

                                     ggCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTggg                                                    <  2:1557373/93‑1 (MQ=255)

                                      gcgcGATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                                                   <  2:2663415/93‑1 (MQ=255)

                                       cgcgATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCg                                                  <  2:1658807/93‑1 (MQ=255)

                                         cgATGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTc                                                <  2:2405698/93‑1 (MQ=255)

                                          gATGCAACTGGCAGGTCCTGGCG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCa                                               >  2:1524484/1‑93 (MQ=255)

                                           aTGCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCaa                                              >  2:404419/1‑93 (MQ=255)

                                             gCAACTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATc                                            >  1:1213829/1‑93 (MQ=255)

                                             gCAACTGGCAGGTCCTGGAG‑TTTTTGG‑CCC‑CCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCa                                           >  2:2371519/1‑93 (MQ=255)

                                               aaCTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCat                                          >  2:3671357/1‑93 (MQ=255)

                                                 cTGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCatat                                        <  1:3244230/93‑1 (MQ=255)

                                                  tGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATAtt                                       <  1:1189430/93‑1 (MQ=255)

                                                  tGGCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATAtt                                       <  1:2578241/93‑1 (MQ=255)

                                                   ggCAGGTCCTGGAG‑TTTTTGGCCCC‑CCCCCGAGATCGCCGCATCGCTATAACCGTGCCGCGTGAGGCGCTGGCATTAGGCGTCAATCATATTg                                      <  2:1357757/93‑1 (MQ=255)

                                                     cAGGTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGAt                                    <  1:3694182/93‑1 (MQ=255)

                                                       ggTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATAc                                  >  1:3795430/1‑93 (MQ=255)

                                                        gTCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATAcc                                 <  2:1317472/93‑1 (MQ=255)

                                                         tCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGTGTCTCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa                                <  2:3862803/93‑1 (MQ=255)

                                                         tCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa                                <  2:154494/93‑1 (MQ=255)

                                                         tCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa                                <  2:154495/93‑1 (MQ=255)

                                                         tCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa                                <  2:3204456/93‑1 (MQ=255)

                                                         tCCTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCACATTGATACCa                                >  2:3948315/1‑93 (MQ=255)

                                                           cTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc                              <  1:1018887/93‑1 (MQ=255)

                                                           cTGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc                              <  1:2230790/93‑1 (MQ=255)

                                                            tGGAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCg                             >  2:297580/1‑93 (MQ=255)

                                                             ggAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGa                            >  2:4192070/1‑93 (MQ=255)

                                                             ggAG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGa                            >  2:4192069/1‑93 (MQ=255)

                                                               aG‑TTTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACt                          <  2:730318/93‑1 (MQ=255)

                                                               aG‑TCTTTGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACt                          <  1:4654548/93‑1 (MQ=255)

                                                                  tttttGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACttt                        <  1:4069105/93‑1 (MQ=255)

                                                                  ttttgGC‑CCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACtttt                       <  1:856999/93‑1 (MQ=255)

                                                                    tttGGCCCC‑CCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTa                      <  1:1189737/93‑1 (MQ=255)

                                                                      tggcccc‑cccACGTGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATg                    <  2:4611995/89‑1 (MQ=255)

                                                                      tggcccc‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATg                    <  1:3842622/89‑1 (MQ=255)

                                                                      tggcccc‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATg                    <  1:4173114/89‑1 (MQ=255)

                                                                       ggcccc‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATgg                   >  2:467672/4‑93 (MQ=255)

                                                                        gcccc‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGt                  <  1:2346067/91‑1 (MQ=255)

                                                                        gcccc‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGt                  <  1:4638269/91‑1 (MQ=255)

                                                                         cccc‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTc                 <  1:2025479/92‑1 (MQ=255)

                                                                           cc‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCg               >  2:1236962/1‑93 (MQ=255)

                                                                            c‑cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGc              >  2:2648022/1‑93 (MQ=255)

                                                                              cccACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCa             >  2:2677999/1‑93 (MQ=255)

                                                                                 aCGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGt          >  2:3487246/1‑93 (MQ=255)

                                                                                  cGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTc         <  1:3689822/93‑1 (MQ=255)

                                                                                   gagaTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCa        <  1:1495166/93‑1 (MQ=255)

                                                                                   gagaTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCa        <  1:1495165/93‑1 (MQ=255)

                                                                                    agaTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCAc       <  2:455095/93‑1 (MQ=255)

                                                                                    agaTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCAc       >  2:2924670/1‑93 (MQ=255)

                                                                                     gaTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCAcc      <  2:1768311/93‑1 (MQ=255)

                                                                                     gaTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCAcc      <  1:3297660/93‑1 (MQ=255)

                                                                                      aTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCAGCa     >  2:2219485/1‑93 (MQ=255)

                                                                                      aTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCa     <  1:1427669/93‑1 (MQ=255)

                                                                                       tCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCaa    >  2:4150636/1‑93 (MQ=255)

                                                                                       tCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCaa    <  1:416348/93‑1 (MQ=255)

                                                                                       tCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCaa    <  2:2273802/93‑1 (MQ=255)

TACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAG‑TTTTTGG‑CCC‑CCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATC  >  NC_000913/1472262‑1472440

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