BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,879,747	C→T	A125T (GCC→ACC)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,879,747	0	C	T	98.9%	1050.3 / NA	281	A125T (GCC→ACC)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand): ref base C (1/0); major base T (109/167); minor base G (3/0); total (114/167)
Fisher's exact test for biased strand distribution p-value = 6.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

GGTAACCGCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCAC  >  NC_000913/3879649‑3879843

ggTAACCGCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGgtgt                                                                                                <  1:5105525/101‑1 (MQ=255)

  tAACCGCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGttt                                                                                              <  1:5105360/101‑1 (MQ=255)

      cGCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTtacta                                                                                          <  1:5105629/101‑1 (MQ=255)

      cGCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTtacta                                                                                          <  1:5105467/101‑1 (MQ=255)

       gCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTtactac                                                                                         >  1:5105914/1‑101 (MQ=255)

       gCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTtactac                                                                                         <  1:5105390/101‑1 (MQ=255)

       gCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTtactac                                                                                         <  2:5105332/101‑1 (MQ=255)

        ccAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTAcc                                                                                        <  1:5105432/101‑1 (MQ=255)

          aGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGa                                                                                      >  2:5105921/1‑101 (MQ=255)

          aGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGa                                                                                      <  1:5105519/101‑1 (MQ=255)

           gCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGaa                                                                                     <  1:5105408/101‑1 (MQ=255)

           gCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGaa                                                                                     <  1:5105587/101‑1 (MQ=255)

            cGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGaaa                                                                                    >  1:5105922/1‑101 (MQ=255)

             ggggCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAac                                                                                   <  1:5105297/101‑1 (MQ=255)

             ggggCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAac                                                                                   <  2:5105372/101‑1 (MQ=255)

                 ccTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCa                                                                               >  1:5105926/1‑101 (MQ=255)

                   tGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAAc                                                                             <  2:5105559/101‑1 (MQ=255)

                   tGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAAc                                                                             <  2:5105598/101‑1 (MQ=255)

                   tGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAAc                                                                             >  1:5105930/1‑101 (MQ=255)

                   tGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAAc                                                                             <  1:5105793/101‑1 (MQ=255)

                     cGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGc                                                                           <  1:5105505/101‑1 (MQ=255)

                        tACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCa                                                                               >  2:5105934/1‑94 (MQ=255)

                        tACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGt                                                                        <  1:5105574/101‑1 (MQ=255)

                          cacCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGc                                                                      >  1:5105937/1‑101 (MQ=255)

                          cacCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGc                                                                      <  2:5105611/101‑1 (MQ=255)

                           acCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCa                                                                     <  1:5105659/101‑1 (MQ=255)

                             cGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGa                                                                   >  2:5105942/1‑101 (MQ=255)

                              gtgtTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGCTTTTGCGACAGGGTGTTTACTACCGAAac                                                                                   <  2:5105446/84‑1 (MQ=255)

                               tgtTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGaaa                                                                                    >  1:5105946/1‑82 (MQ=255)

                                gtTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGAccg                                                                <  1:5105513/101‑1 (MQ=255)

                                 ttCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGAccgc                                                               >  2:5105947/1‑101 (MQ=255)

                                   cGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTTCGACAGGGTGTTtactac                                                                                         >  2:5105951/1‑73 (MQ=255)

                                   cGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGAccgccg                                                             <  1:5105546/101‑1 (MQ=255)

                                    gTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCgg                                                            <  1:5105610/101‑1 (MQ=255)

                                     tAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGa                                                           >  2:5105958/1‑101 (MQ=255)

                                     tAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGa                                                           <  2:5105329/101‑1 (MQ=255)

                                     tAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACAGAAACACCAACGCCGTGCAGACCGCCGGa                                                           >  1:5105955/1‑101 (MQ=255)

                                      aGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGac                                                          <  2:5105492/101‑1 (MQ=255)

                                       gATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGaca                                                         <  2:5105460/101‑1 (MQ=255)

                                       gATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGaca                                                         <  1:5105570/101‑1 (MQ=255)

                                       gATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGac                                                          >  1:5105961/1‑100 (MQ=255)

                                         tCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACt                                                       <  1:5105473/101‑1 (MQ=255)

                                          cTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACtt                                                      <  1:5105742/101‑1 (MQ=255)

                                          cTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACtt                                                      >  1:5105965/1‑101 (MQ=255)

                                           tGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCaa                                                                              <  2:5105687/76‑1 (MQ=255)

                                             aCGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTtat                                                   >  2:5105969/1‑101 (MQ=255)

                                              cGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCCCCGGACACTTtata                                                  >  1:5105971/1‑101 (MQ=255)

                                               ggTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAg                                                 <  1:5105498/101‑1 (MQ=255)

                                                gTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAgg                                                <  2:5105631/101‑1 (MQ=255)

                                                    aTTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGtt                                            <  1:5105482/101‑1 (MQ=255)

                                                      tttACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAgg                                                >  2:5105978/1‑95 (MQ=255)

                                                      tttACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTAt                                          >  2:5105975/1‑101 (MQ=255)

                                                       ttACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGGGTTTACTACCg                                                                                       >  1:5105979/1‑55 (MQ=255)

                                                        tACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCg                                        <  2:5105699/101‑1 (MQ=255)

                                                         aCCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGt                                       <  2:5105584/101‑1 (MQ=255)

                                                          cccTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGt                                       <  1:5105209/100‑1 (MQ=255)

                                                          cccTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGt                                       <  2:5105437/100‑1 (MQ=255)

                                                          cccTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTc                                      <  2:5105548/101‑1 (MQ=255)

                                                           ccTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCa                                     >  2:5105983/1‑101 (MQ=255)

                                                           ccTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCa                                     <  2:5105375/101‑1 (MQ=255)

                                                             tCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCaaa                                   >  2:5105986/1‑101 (MQ=255)

                                                              cgcgCTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAAt                                  <  2:5105651/101‑1 (MQ=255)

                                                               gcgcTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAAtt                                 <  2:5105481/101‑1 (MQ=255)

                                                                cgcTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGAc                                                                  >  2:5105990/1‑67 (MQ=255)

                                                                 gcTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTa                               >  1:5105991/1‑101 (MQ=255)

                                                                  cTGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTAc                              <  1:5105596/101‑1 (MQ=255)

                                                                   tGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGAccgcc                                                              >  2:5105997/1‑68 (MQ=255)

                                                                   tGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTAcc                             <  2:5105530/101‑1 (MQ=255)

                                                                   tGGATAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTAcc                             >  1:5105996/1‑101 (MQ=255)

                                                                      aTAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGcc                          >  1:5106000/1‑101 (MQ=255)

                                                                       tAACCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCt                         <  1:5105680/101‑1 (MQ=255)

                                                                        aaCCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTg                        >  2:5106005/1‑101 (MQ=255)

                                                                        aaCAAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTg                        >  2:5106003/1‑101 (MQ=255)

                                                                         aCCAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTGTAGGAGTTATCGTCAAATTTACCGCCTGc                       <  2:5105388/101‑1 (MQ=255)

                                                                          ccAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTg                                                                       <  2:5105515/52‑1 (MQ=255)

                                                                           cAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGt                     <  1:5105635/101‑1 (MQ=255)

                                                                           cAGCTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGt                     >  1:5106009/1‑101 (MQ=255)

                                                                              cTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCa                  <  2:5105401/101‑1 (MQ=255)

                                                                              cTCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCa                  <  1:5105710/101‑1 (MQ=255)

                                                                               tCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTtata                                                  <  2:5105269/68‑1 (MQ=255)

                                                                               tCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAg                 <  1:5105810/101‑1 (MQ=255)

                                                                               tCCAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAg                 <  1:5105560/101‑1 (MQ=255)

                                                                                ccAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGa                >  1:5106013/1‑101 (MQ=255)

                                                                                 cAGTTTTTGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGaa               >  2:5106014/1‑101 (MQ=255)

                                                                                    tttttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACgg            >  2:5106020/1‑101 (MQ=255)

                                                                                    tttttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACgg            <  2:5105472/101‑1 (MQ=255)

                                                                                    tttttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACgg            >  1:5106019/1‑101 (MQ=255)

                                                                                     ttttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGt           <  2:5105886/101‑1 (MQ=255)

                                                                                     ttttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGt           <  2:5105665/101‑1 (MQ=255)

                                                                                      tttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTc          <  2:5105806/101‑1 (MQ=255)

                                                                                      tttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTc          >  2:5106023/1‑101 (MQ=255)

                                                                                       ttGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCa         <  2:5105734/101‑1 (MQ=255)

                                                                                        tGCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCAt        <  1:5105620/101‑1 (MQ=255)

                                                                                         gCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATg       <  2:5105701/101‑1 (MQ=255)

                                                                                         gCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATg       >  1:5106027/1‑101 (MQ=255)

                                                                                         gCGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATg       <  1:5105427/101‑1 (MQ=255)

                                                                                          cGACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGa      >  2:5106028/1‑101 (MQ=255)

                                                                                           gACAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGAt     >  1:5106031/1‑101 (MQ=255)

                                                                                             cAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCa   <  2:5105623/101‑1 (MQ=255)

                                                                                             cAGGGTGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCa   >  1:5106035/1‑101 (MQ=255)

GGTAACCGCCAGCGGGGCCTGCGGTACACCGTGTTCGTAGATCTGACGGTGAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCAC  >  NC_000913/3879649‑3879843

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