BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,140,841	A→G	V917A (GTG→GCG)	rne ←	RNase E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,140,841	0	A	G	99.5%	607.4 / ‑4.9	196	V917A (GTG→GCG)	rne	RNase E
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (97/98); total (97/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

Only 100 of 208 total aligned reads displayed.

ACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACT  >  NC_000913/1140751‑1140930

accGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGAATCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACt                                                                                         >  1:1970453/1‑93 (MQ=255)

  cGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTc                                                                                       <  1:4858536/93‑1 (MQ=255)

  cGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTc                                                                                       >  2:3998606/1‑93 (MQ=255)

   ggTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCa                                                                                      >  1:5418176/1‑93 (MQ=255)

    gTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAg                                                                                     <  1:1217444/93‑1 (MQ=255)

    gTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAg                                                                                     >  2:4996384/1‑93 (MQ=255)

    gTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAg                                                                                     >  1:6792812/1‑93 (MQ=255)

    gTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAg                                                                                     >  1:7474329/1‑93 (MQ=255)

     ttCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAgg                                                                                    >  2:7632400/1‑93 (MQ=255)

     ttCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAgg                                                                                    <  2:1165175/93‑1 (MQ=255)

      tCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGa                                                                                   <  2:1641303/93‑1 (MQ=255)

      tCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGa                                                                                   <  2:6207865/93‑1 (MQ=255)

      tCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGa                                                                                   <  1:4161353/93‑1 (MQ=255)

       ttGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGTCTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGAt                                                                                  <  1:6051989/92‑1 (MQ=255)

        tGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGGTACCGCGGCAGCGATCGCTTCAAGATg                                                                                 >  1:2390296/1‑93 (MQ=255)

         gCTTGTTCTGTAGCTTCCTGGGCTACGGTTACATCAGACTCGGTAATCGCCTGCGGCTTTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGc                                                                                <  2:2163213/93‑1 (MQ=255)

         gCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGc                                                                                >  1:4335998/1‑93 (MQ=255)

          cTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcg                                                                               <  1:2763381/93‑1 (MQ=255)

          cTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcg                                                                               >  2:4377466/1‑93 (MQ=255)

          cTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcg                                                                               <  1:400712/93‑1 (MQ=255)

           ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcgt                                                                              >  1:7455860/1‑93 (MQ=255)

           ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcgt                                                                              >  2:5066352/1‑93 (MQ=255)

            tGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcgtc                                                                             >  2:5709223/1‑93 (MQ=255)

            tGTACTGCAACTTCCTGGGATGCGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTCCAGGATAcgtc                                                                             <  2:5238355/93‑1 (MQ=255)

              ttCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcgtcgt                                                                           >  2:1979484/1‑93 (MQ=255)

              ttCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGcgtcgt                                                                           >  1:7814860/1‑93 (MQ=255)

               tCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGtt                                                                          <  2:2010230/93‑1 (MQ=255)

                 tGCAACTTCCTGGGCTACGGCAACATCAGGCTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTc                                                                        >  2:1342955/1‑93 (MQ=255)

                 tGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTc                                                                        <  2:7775921/93‑1 (MQ=255)

                 tGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTc                                                                        >  1:4810463/1‑93 (MQ=255)

                   cAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCaa                                                                      <  2:4656348/93‑1 (MQ=255)

                    aaCTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAAc                                                                     >  1:1816244/1‑93 (MQ=255)

                    aaCTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAAc                                                                     >  2:6310937/1‑93 (MQ=255)

                      cctCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCa                                                                   <  2:5246953/91‑1 (MQ=255)

                      cTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCa                                                                   >  1:4629305/1‑93 (MQ=255)

                      cTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCa                                                                   <  1:60180/93‑1 (MQ=255)

                          cTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTc                                                               <  2:2811334/93‑1 (MQ=255)

                           tGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTTAACCACTTcc                                                              >  1:1990151/1‑93 (MQ=255)

                              gCTCCGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGt                                                           <  2:1655691/93‑1 (MQ=255)

                              gCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGt                                                           >  2:2991366/1‑93 (MQ=255)

                               cTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGtt                                                          <  1:3271578/93‑1 (MQ=255)

                                 aCGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTgg                                                        >  2:3396272/1‑92 (MQ=255)

                                 aCGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGc                                                        <  2:1819423/93‑1 (MQ=255)

                                   ggCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCgg                                                      >  2:5462478/1‑93 (MQ=255)

                                   ggCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCgg                                                      <  2:7657114/93‑1 (MQ=255)

                                   ggCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCgg                                                      <  1:5516323/93‑1 (MQ=255)

                                    gCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGt                                                     >  1:3582327/1‑93 (MQ=255)

                                     cAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGtt                                                    <  1:3697673/93‑1 (MQ=255)

                                     cAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGtt                                                    <  1:5086881/93‑1 (MQ=255)

                                      aaCATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTc                                                   <  2:936039/93‑1 (MQ=255)

                                      aaCATCAGACTCGGTAATCACCTGCGGCTGTTCAGCTACCGCGGCAGCGATCGCTTCACCATGCGTCGTTTCAACCTCTTCCCGCTGCGGTTc                                                   >  2:7556526/1‑93 (MQ=255)

                                       aCATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCg                                                  >  1:2121187/1‑93 (MQ=255)

                                        cATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGTATGCCTCGTTTCAACAACTTCCTGTTGCGGTTCga                                                 >  1:2646433/1‑92 (MQ=255)

                                        cATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCgg                                                 >  1:6879440/1‑93 (MQ=255)

                                          tCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCa                                               >  2:6861490/1‑93 (MQ=255)

                                          tCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCa                                               >  2:6744358/1‑93 (MQ=255)

                                           cAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCaa                                              <  2:6001147/93‑1 (MQ=255)

                                            aGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAAc                                             <  1:3175080/93‑1 (MQ=255)

                                            aGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAAc                                             >  2:3180979/1‑93 (MQ=255)

                                                tCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGa                                         >  1:6893655/1‑93 (MQ=255)

                                                 cGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGaa                                        <  2:6599476/93‑1 (MQ=255)

                                                 cGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGaa                                        <  2:6251191/93‑1 (MQ=255)

                                                  ggTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAAc                                       >  2:6908618/1‑93 (MQ=255)

                                                  ggTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAAc                                       >  1:6745733/1‑93 (MQ=255)

                                                    tAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACgg                                     >  2:6037567/1‑93 (MQ=255)

                                                     aaTCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg                                    >  2:5948072/1‑93 (MQ=255)

                                                     aaTCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg                                    <  1:4131923/93‑1 (MQ=255)

                                                      aTCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACgggg                                   >  1:6553736/1‑93 (MQ=255)

                                                       tCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGACCggggg                                  >  1:4030410/1‑93 (MQ=255)

                                                       tCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggggg                                  >  1:222841/1‑93 (MQ=255)

                                                        cACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGc                                 <  2:407346/93‑1 (MQ=255)

                                                        cACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGc                                 >  2:4830798/1‑93 (MQ=255)

                                                        cACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGc                                 <  1:4290782/93‑1 (MQ=255)

                                                           cTGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTc                              <  1:3974113/93‑1 (MQ=255)

                                                            tGCGGCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCt                             >  1:7241162/1‑93 (MQ=255)

                                                              cGGCTGTACAGTTACCGAGGCAGCAATCGATTCAGAATGCATCGTTTCATCCTCTTCCGGTTGCGGTTAGGCAACATGAACGGGGGCTTCTaa                           <  2:7012044/93‑2 (MQ=255)

                                                                gCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGa                         >  2:514769/1‑93 (MQ=255)

                                                                gCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGa                         >  2:4732896/1‑93 (MQ=255)

                                                                gCTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGa                         >  2:2425298/1‑93 (MQ=255)

                                                                 cTGTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGac                        >  1:616752/1‑93 (MQ=255)

                                                                   gTTCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGacac                      >  2:7333797/1‑93 (MQ=255)

                                                                    ttCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACAcc                     <  2:4168056/93‑1 (MQ=255)

                                                                     tCAGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCg                    >  1:1458096/1‑93 (MQ=255)

                                                                       aGTTCCCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGc                  <  1:5618743/93‑1 (MQ=255)

                                                                       aGTTACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGc                  <  1:4561714/93‑1 (MQ=255)

                                                                          tACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCAc               >  2:588158/1‑93 (MQ=255)

                                                                          tACCGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCAc               >  1:3925685/1‑93 (MQ=255)

                                                                            ccGCGGCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCActt             <  2:2217938/93‑1 (MQ=255)

                                                                                ggCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACttcttc         <  1:5756732/93‑1 (MQ=255)

                                                                                ggCAGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACttcttc         >  1:2834041/1‑93 (MQ=255)

                                                                                  cAGCGACCGCTCCAGGATGCGTCGTTTCAACCACTCCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTcaa       <  2:6432636/93‑1 (MQ=255)

                                                                                   aGCGATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTcaac      <  1:763422/93‑1 (MQ=255)

                                                                                      gATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCaacaac   <  1:3458698/93‑1 (MQ=255)

                                                                                      gATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCaacaac   >  1:6980867/1‑93 (MQ=255)

                                                                                      gATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCaacaac   >  1:6587306/1‑93 (MQ=255)

                                                                                      gATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCaacaac   <  1:4987526/93‑1 (MQ=255)

                                                                                      gATCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCaacaac   <  1:5885457/93‑1 (MQ=255)

                                                                                       aTCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACt  >  2:1478230/1‑93 (MQ=255)

                                                                                       aTCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACt  <  2:3509037/93‑1 (MQ=255)

                                                                                       aTCGCTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACt  <  2:3898969/93‑1 (MQ=255)

ACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACT  >  NC_000913/1140751‑1140930

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