BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	979,652	T→C	G1368G (GGT→GGC)	mukB →	MukB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	979,652	0	T	C	83.2%	528.5 / 99.6	244	G1368G (GGT→GGC)	mukB	MukB
Reads supporting (aligned to +/- strand): ref base T (23/18); new base C (92/111); total (115/129)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 246 total aligned reads displayed.

TGGATTACCGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTAAAGATATCTCTCC  >  NC_000913/979561‑979744

  gATTACCGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAg                                                                                           <  1:4494934/93‑1 (MQ=255)

  gATTACCGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAg                                                                                           <  1:1423461/93‑1 (MQ=255)

    ttACCGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGc                                                                                         <  1:2178859/93‑1 (MQ=255)

     tACCGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCg                                                                                        <  2:4700279/93‑1 (MQ=255)

       ccGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGAt                                                                                      <  1:688178/93‑1 (MQ=255)

        cGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGAtt                                                                                     <  1:2762527/93‑1 (MQ=255)

         gCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTg                                                                                    <  2:6994360/93‑1 (MQ=255)

         gCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTg                                                                                    <  1:5913542/93‑1 (MQ=255)

         gCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTg                                                                                    >  1:773285/1‑93 (MQ=255)

         gCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTg                                                                                    >  1:3410165/1‑93 (MQ=255)

           aaCTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGt                                                                                  <  1:1839047/93‑1 (MQ=255)

            aCTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTa                                                                                 >  2:1291606/1‑93 (MQ=255)

            aCTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCACTGTCGACCGGTGAGGCGATTGGTa                                                                                 <  1:3124904/93‑1 (MQ=255)

             cTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTAc                                                                                >  1:2814871/1‑93 (MQ=255)

              tATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTAcc                                                                               >  1:3827932/1‑93 (MQ=255)

              tATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTAcc                                                                               >  1:960077/1‑93 (MQ=255)

              tATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTAcc                                                                               <  2:1937201/93‑1 (MQ=255)

                  tGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTa                                                                           <  1:5584099/93‑1 (MQ=255)

                  tGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTa                                                                           >  1:6520839/1‑93 (MQ=255)

                   ggAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTAt                                                                          <  2:6841577/93‑1 (MQ=255)

                     aaaTGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGCCCGGCGAGGCGATTGGTACCGGTATGt                                                                        <  1:6018528/93‑1 (MQ=255)

                     aaaTGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGt                                                                        <  2:2941154/93‑1 (MQ=255)

                     aaaTGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGt                                                                        >  1:5377764/1‑93 (MQ=255)

                      aaTGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTc                                                                       >  1:3734871/1‑93 (MQ=255)

                        tGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATTTCGa                                                                     >  2:3237317/1‑93 (MQ=255)

                        tGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGa                                                                     >  2:6646925/1‑93 (MQ=255)

                        tGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGa                                                                     >  1:7421676/1‑93 (MQ=255)

                           aaGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTc                                                                  <  2:4346294/93‑1 (MQ=255)

                            aGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCt                                                                 >  1:2368246/1‑93 (MQ=255)

                             gTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTg                                                                >  2:1574/1‑93 (MQ=255)

                             gTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTg                                                                >  1:6552056/1‑93 (MQ=255)

                             gTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTg                                                                >  2:7275080/1‑93 (MQ=255)

                              ttGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTgg                                                               >  2:1676976/1‑93 (MQ=255)

                              ttGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTgg                                                               <  1:302785/93‑1 (MQ=255)

                              ttGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTgg                                                               <  1:1021258/93‑1 (MQ=255)

                                gAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTg                                                             >  1:2206172/1‑93 (MQ=255)

                                  ggTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAt                                                           >  2:5393648/1‑93 (MQ=255)

                                  ggTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAt                                                           >  1:4889563/1‑93 (MQ=255)

                                    ttAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAtgg                                                         >  1:6630172/1‑93 (MQ=255)

                                    ttAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAtgg                                                         >  2:934637/1‑93 (MQ=255)

                                    ttAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAtgg                                                         <  1:4696072/93‑1 (MQ=255)

                                      aaCCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAtggtg                                                       <  1:1196605/93‑1 (MQ=255)

                                      aaCCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAtggtg                                                       <  2:711212/93‑1 (MQ=255)

                                       aCCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGAtggtgg                                                      >  1:2223331/1‑93 (MQ=255)

                                         cGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTa                                                    <  1:994632/93‑1 (MQ=255)

                                         cGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTa                                                    >  2:2520988/1‑93 (MQ=255)

                                         cGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTa                                                    >  2:3441642/1‑93 (MQ=255)

                                          gTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTAc                                                   >  2:1047688/1‑93 (MQ=255)

                                          gTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTAc                                                   >  2:3870267/1‑93 (MQ=255)

                                           tGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACa                                                  >  2:6020811/1‑93 (MQ=255)

                                            ggTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACaa                                                 <  2:5001181/93‑1 (MQ=255)

                                             gTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACaaa                                                <  2:4284701/93‑1 (MQ=255)

                                             gTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACaaa                                                <  2:4871351/93‑1 (MQ=255)

                                             gTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACaaa                                                >  1:3087841/1‑93 (MQ=255)

                                             gTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACaaa                                                <  1:1738954/93‑1 (MQ=255)

                                             gTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACaaa                                                <  1:1091233/93‑1 (MQ=255)

                                              ttCCGATGGCTGGCTGCGGCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACaaa                                                <  1:6185425/93‑1 (MQ=255)

                                              ttCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAg                                               >  1:3197838/1‑93 (MQ=255)

                                                ccGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCt                                             <  2:1145671/93‑1 (MQ=255)

                                                  gATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTgg                                           >  1:2911590/1‑93 (MQ=255)

                                                  gATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTgg                                           <  1:1608735/93‑1 (MQ=255)

                                                  gATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTgg                                           <  2:6820880/93‑1 (MQ=255)

                                                  gATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATATCGATTCTGGTGATGGTGGTACAAAGCTgg                                           >  2:51256/1‑93 (MQ=255)

                                                   aTGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTggg                                          <  1:541383/93‑1 (MQ=255)

                                                   aTGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTggg                                          >  2:3732265/1‑93 (MQ=255)

                                                    tggctggcTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGa                                         <  1:2493396/93‑1 (MQ=255)

                                                    tggctggcTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGa                                         >  2:2846671/1‑93 (MQ=255)

                                                      gctggctgC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAg                                       >  1:1543989/1‑93 (MQ=255)

                                                       ctggctgC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGa                                      <  1:2553981/93‑1 (MQ=255)

                                                       ctggctgC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGa                                      <  2:1815079/93‑1 (MQ=255)

                                                        tggctgC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGAt                                     <  1:209860/93‑1 (MQ=255)

                                                        tggctgC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGAt                                     <  1:1327438/93‑1 (MQ=255)

                                                         ggctgC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATg                                    <  2:5620146/93‑1 (MQ=255)

                                                          gctgC‑GCGCAGAGTCTGGTGCATTGTTGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGa                                   >  1:546225/1‑93 (MQ=255)

                                                            tgC‑GCGCAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAAt                                 <  1:4226147/93‑1 (MQ=255)

                                                             gc‑gcgcAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAAtc                                <  1:4521002/93‑1 (MQ=255)

                                                                gcgcAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAAtctc                              <  2:1868284/93‑1 (MQ=255)

                                                                 cgcAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcg                             >  1:5851969/1‑93 (MQ=255)

                                                                  tcAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcgc                            <  1:5202203/92‑1 (MQ=255)

                                                                  gcAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcgc                            <  1:719400/93‑1 (MQ=255)

                                                                  gcAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcgc                            <  1:6776844/93‑1 (MQ=255)

                                                                  gcAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcgc                            <  1:363996/93‑1 (MQ=255)

                                                                   cAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcgcc                           <  2:1970081/93‑1 (MQ=255)

                                                                   cAGAGTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcgcc                           <  2:6360935/93‑1 (MQ=255)

                                                                     gagTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTcgccgc                         >  1:4417145/1‑93 (MQ=255)

                                                                      agTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCgccgcc                        <  1:5618915/93‑1 (MQ=255)

                                                                       gTCTGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCt                       >  1:6761019/1‑93 (MQ=255)

                                                                          tGGTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTgcg                    >  1:810685/1‑93 (MQ=255)

                                                                           ggTGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTgcgc                   <  1:4127303/93‑1 (MQ=255)

                                                                             tGCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCgg                 <  2:2220763/93‑1 (MQ=255)

                                                                              gCATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGt                >  2:5643771/1‑93 (MQ=255)

                                                                               cATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTTTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTa               >  1:6863283/1‑93 (MQ=255)

                                                                               cATTGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTa               >  1:7321318/1‑93 (MQ=255)

                                                                                 ttGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTaaa             <  1:495050/93‑1 (MQ=255)

                                                                                 ttGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTaaa             <  2:2721469/93‑1 (MQ=255)

                                                                                 ttGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTaaa             >  1:4643371/1‑93 (MQ=255)

                                                                                 ttGTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTaaa             >  1:3922254/1‑93 (MQ=255)

                                                                                   gTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTAAAGa           <  2:6182483/93‑1 (MQ=255)

                                                                                   gTCGACCGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTAAAGa           >  1:5743539/1‑93 (MQ=255)

                                                                                         cGGCGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCGCGCCGCCTGCGCGGTAAAGATAtctc     <  1:6314894/93‑1 (MQ=255)

TGGATTACCGCAACTATCTGGAAATGGAAGTTGAGGTTAACCGTGGTTCCGATGGCTGGCTGC‑GCGCAGAGTCTGGTGCATTGTCGACCGGTGAGGCGATTGGTACCGGTATGTCGATTCTGGTGATGGTGGTACAAAGCTGGGAAGATGAATCTCGCCGCCTGCGCGGTAAAGATATCTCTCC  >  NC_000913/979561‑979744

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