BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	484,877	G→A	15.6%	T248I (ACC→ATC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	484,877	0	G	A	15.6%	394.5 / 47.5	158	T248I (ACC→ATC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base G (69/63); new base A (4/21); total (73/85)
Fisher's exact test for biased strand distribution p-value = 8.54e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

GGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACAT  >  NC_000913/484778‑484976

gggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTgg                                                                                                     >  1:1391921/1‑100 (MQ=255)

gggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTgg                                                                                                     <  1:154920/100‑1 (MQ=255)

   aaGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGa                                                                                                  <  1:247650/100‑1 (MQ=255)

     gATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATc                                                                                                <  2:506812/100‑1 (MQ=255)

     gATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATc                                                                                                >  1:1075178/1‑100 (MQ=255)

         gcgcGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTg                                                                                            <  2:2388969/100‑1 (MQ=255)

         gcgcGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTg                                                                                            <  2:370770/100‑1 (MQ=255)

         gcgcGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTg                                                                                            <  1:769667/100‑1 (MQ=255)

           gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGac                                                                                          <  1:2599331/100‑1 (MQ=255)

            cgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacg                                                                                         <  2:1667315/100‑2 (MQ=255)

             gTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGacac                                                                                        >  1:199108/1‑100 (MQ=255)

             gTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacgg                                                                                        <  2:2520110/100‑3 (MQ=255)

             gTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacgg                                                                                        <  2:3292608/100‑3 (MQ=255)

             gTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacgg                                                                                        <  2:1470058/100‑3 (MQ=255)

             gTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacgg                                                                                        <  1:1261304/100‑3 (MQ=255)

              tAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacggc                                                                                       <  2:1923887/100‑4 (MQ=255)

                gggTGATAGGCCCAGTGGTCTGATCAACGGTAATTTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacggcat                                                                                     <  2:1198855/100‑6 (MQ=255)

                gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACttt                                                                                     >  2:2709203/1‑100 (MQ=255)

                gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacggcat                                                                                     <  1:3538361/100‑6 (MQ=255)

                gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacggc                                                                                       <  1:1881976/98‑4 (MQ=255)

                 ggtggtAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTg                                                                                    >  2:425760/1‑100 (MQ=255)

                 ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTg                                                                                    >  1:2560316/1‑100 (MQ=255)

                    gATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCt                                                                                 <  1:2392952/100‑1 (MQ=255)

                    gATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCt                                                                                 >  1:2493879/1‑100 (MQ=255)

                    gATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacggcattaag                                                                                 <  1:1170644/100‑10 (MQ=37)

                      tAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCttt                                                                               >  2:3085287/1‑100 (MQ=255)

                      tAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCttt                                                                               >  1:343711/1‑100 (MQ=255)

                      tAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCttt                                                                               >  2:8862/1‑100 (MQ=255)

                      tAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTACTGGGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCtt                                                                                >  2:710240/1‑99 (MQ=255)

                       aGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTg                                                                              <  2:1426330/100‑1 (MQ=255)

                        gACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                             <  1:2802645/100‑1 (MQ=255)

                        gACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                             >  2:3783602/1‑100 (MQ=255)

                        gACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                             <  1:3372639/100‑1 (MQ=255)

                        gACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                             <  1:1451678/100‑1 (MQ=255)

                        gACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                             >  1:1413724/1‑100 (MQ=255)

                         aCCCAGTGGTCTGATCAGCGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGcc                                                                            <  1:716241/100‑1 (MQ=255)

                         aCCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGcc                                                                            <  1:55857/100‑1 (MQ=255)

                           ccAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGt                                                                          <  2:3127481/100‑1 (MQ=255)

                           ccAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGt                                                                          <  2:582916/100‑1 (MQ=255)

                            cAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGtt                                                                         >  2:2026239/1‑100 (MQ=255)

                                 gTCTGATCAACGGTGACTTCAGAGAATTCCATCGTACCGTCCTGCGGTAACGGGTGGCCGTCACTGGGGAGCAGTGACACTGTGGCTTTGCCGTTctct                                                                     <  2:3031345/99‑1 (MQ=255)

                                 gTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTctct                                                                     >  1:1881616/1‑99 (MQ=255)

                                  tCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTg                                                                   <  1:2785515/100‑1 (MQ=255)

                                  tCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTg                                                                   <  1:2195214/100‑1 (MQ=255)

                                     gATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGttt                                                                >  2:2964054/1‑100 (MQ=255)

                                     gATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGttt                                                                >  1:810354/1‑100 (MQ=255)

                                        cAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAg                                                             >  1:1312770/1‑100 (MQ=255)

                                        cAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAg                                                             <  2:1986374/100‑1 (MQ=255)

                                         aaCGGTAACCTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGc                                                            >  2:2495928/1‑100 (MQ=255)

                                          aCGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCg                                                           >  2:1305120/1‑100 (MQ=255)

                                           cGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGt                                                          >  1:371252/1‑100 (MQ=255)

                                                  gTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACAGTTTGGATTTGCCGTACTCTTGTTTCAGCGTGCCATTc                                                   >  2:2579516/1‑100 (MQ=255)

                                                    cAGAGCATTCCTTCTTACGGTCCTGCGTGAACTTATTGCCTTCTCTGGTGTTCATTTTCACTTTGGCTTTGCCGTTCTTTTTTTTCTGCGTGCCATTCgt                                                 <  2:3235649/100‑2 (MQ=255)

                                                      gagaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCa                                               <  2:1807496/100‑1 (MQ=255)

                                                        gaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGt                                             <  1:740132/100‑1 (MQ=255)

                                                          aTTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTc                                           <  2:2984340/100‑1 (MQ=255)

                                                           ttCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTcc                                          <  1:2014142/100‑1 (MQ=255)

                                                           ttCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTcc                                          <  1:2005420/100‑1 (MQ=255)

                                                           ttCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTcc                                          >  1:426823/1‑100 (MQ=255)

                                                             ccAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTg                                        >  1:981908/1‑100 (MQ=255)

                                                             ccAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTg                                        >  2:643123/1‑100 (MQ=255)

                                                                    aCCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAgg                                 <  2:3756918/100‑1 (MQ=255)

                                                                     ccGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGgc                                <  1:2390994/100‑1 (MQ=255)

                                                                     ccGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGgc                                >  1:294229/1‑100 (MQ=255)

                                                                      cGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGgcg                               <  1:3372330/100‑1 (MQ=255)

                                                                           tGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGa                          >  2:3590203/1‑100 (MQ=255)

                                                                            gCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGaa                         >  2:1714936/1‑100 (MQ=255)

                                                                            gCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGaa                         >  2:222067/1‑100 (MQ=255)

                                                                            gCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTCCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGaa                         >  2:2109705/1‑100 (MQ=255)

                                                                             cGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAg                        >  1:2668430/1‑100 (MQ=255)

                                                                             cGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAg                        >  2:1557985/1‑100 (MQ=255)

                                                                                gAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCg                     <  1:1923829/100‑1 (MQ=255)

                                                                                 aaCTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGt                    >  2:1764852/1‑100 (MQ=255)

                                                                                  aCTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGtt                   >  2:3352385/1‑100 (MQ=255)

                                                                                   cTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCTTTCTCTTGTGTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGtg                   >  2:2756291/1‑98 (MQ=255)

                                                                                   cTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTg                  >  1:959728/1‑100 (MQ=255)

                                                                                   cTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTg                  >  1:1541870/1‑100 (MQ=255)

                                                                                   cTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTg                  >  1:2891141/1‑100 (MQ=255)

                                                                                   cTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTg                  >  2:3228940/1‑100 (MQ=255)

                                                                                   agccaaagtgTCACTGATCACCAGTGACGCTTTGGCTTTGCCGTTCTCTTGTTTCATCGTGCCATTCGCCAGGTCCTGTTTCCGGCGCAGGAAGTCGTTg                  >  2:1620207/10‑100 (MQ=25)

                                                                                    ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGc                 <  2:691593/100‑1 (MQ=255)

                                                                                    ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGc                 >  1:1124026/1‑100 (MQ=255)

                                                                                     tAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCt                >  1:3795958/1‑100 (MQ=255)

                                                                                       aTGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTgg              <  1:1352239/100‑1 (MQ=255)

                                                                                        tGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGa             >  1:1664965/1‑100 (MQ=255)

                                                                                        aagtgTCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGa             >  2:3042271/5‑100 (MQ=255)

                                                                                         gCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGAc            >  1:2050831/1‑100 (MQ=255)

                                                                                          ccGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACt           <  2:1413725/100‑1 (MQ=255)

                                                                                             tCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTggg        <  1:287294/100‑1 (MQ=255)

                                                                                             tCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTggg        >  1:1868863/1‑100 (MQ=255)

                                                                                             tCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTagg        >  2:187578/1‑100 (MQ=255)

                                                                                              cACTGGTTATCAGTGACACTTTAGCTTTGCCGTTCTCTTGTTTCAGCGTGCTATTTGCCATTACCCGTTTCTGGCTCAGGGATTCTTTGCTgg              <  2:2200279/93‑1 (MQ=255)

                                                                                               aCTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTc      <  2:3500706/100‑1 (MQ=255)

                                                                                               aCTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGGTGCTGGACTGGGTc      <  2:775738/100‑1 (MQ=255)

                                                                                                cTGGTGATCAGTGACACTTTGGCTTTGCCGGTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTca     >  2:3221894/1‑100 (MQ=255)

                                                                                                  ggTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca   >  1:650884/1‑100 (MQ=255)

                                                                                                  ggTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca   <  1:3817988/100‑1 (MQ=255)

                                                                                                  ggTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca   >  2:3347510/1‑100 (MQ=255)

                                                                                                  ggTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca   >  1:1222198/1‑100 (MQ=255)

                                                                                                   gTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACAt  <  1:2037063/100‑1 (MQ=255)

GGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACAT  >  NC_000913/484778‑484976

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