BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,235,117	C→T	A219A (GCG→GCA)	mglC ←	MglC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,235,117	0	C	T	87.1%	689.0 / 68.9	233	A219A (GCG→GCA)	mglC	MglC
Reads supporting (aligned to +/- strand): ref base C (13/17); new base T (103/100); total (116/117)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 245 total aligned reads displayed.

GCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGC  >  NC_000913/2235025‑2235206

gCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGcc                                                                                           >  1:3687529/1‑93 (MQ=255)

  ttCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTtgctgc                                                                                         >  2:6155303/1‑93 (MQ=255)

  ttCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTtgctgc                                                                                         <  1:3985914/93‑1 (MQ=255)

   tCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTgctgct                                                                                        >  1:1901471/1‑93 (MQ=255)

   tCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTgctgct                                                                                        <  2:1435086/93‑1 (MQ=255)

    cTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCtt                                                                                       <  1:5205968/93‑1 (MQ=255)

    cTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCtt                                                                                       >  1:785868/1‑93 (MQ=255)

     tAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATAGTTTTGCTGCTTc                                                                                      <  2:4524411/93‑1 (MQ=255)

     tAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTc                                                                                      >  1:2114230/1‑93 (MQ=255)

     tAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTc                                                                                      <  2:5017383/93‑1 (MQ=255)

     tAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTc                                                                                      <  2:3603740/93‑1 (MQ=255)

      aaCATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTcc                                                                                     <  2:1907415/93‑1 (MQ=255)

       aCATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATTCTTTTGCTGCTTCCg                                                                                    <  2:2816723/93‑1 (MQ=255)

       aCATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCg                                                                                    <  2:3083915/93‑1 (MQ=255)

       aCATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCg                                                                                    >  1:70913/1‑93 (MQ=255)

        cATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCgg                                                                                   >  1:4500531/1‑93 (MQ=255)

          tCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGt                                                                                 <  1:351002/93‑1 (MQ=255)

           cccGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGtt                                                                                <  2:4297961/93‑1 (MQ=255)

             cGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTAc                                                                              >  1:6632396/1‑93 (MQ=255)

             cGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTAc                                                                              >  1:5625641/1‑93 (MQ=255)

             cGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTAc                                                                              >  1:3027948/1‑93 (MQ=255)

              gCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTAcc                                                                             <  1:4924094/93‑1 (MQ=255)

               ccAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCg                                                                            >  1:1037674/1‑93 (MQ=255)

               ccAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCg                                                                            <  2:4328372/93‑1 (MQ=255)

                cAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGc                                                                           >  2:4649944/1‑93 (MQ=255)

                cAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGc                                                                           <  2:743048/93‑1 (MQ=255)

                cAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGc                                                                           >  2:6321127/1‑93 (MQ=255)

                 aaaGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGcc                                                                          >  2:6334222/1‑93 (MQ=255)

                 aaaGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGcc                                                                          <  2:1405160/93‑1 (MQ=255)

                  aaGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCa                                                                         <  2:2247990/93‑1 (MQ=255)

                   aGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGccca                                                                        >  2:3695131/1‑91 (MQ=255)

                   aGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCaa                                                                        <  1:4823721/93‑1 (MQ=255)

                    ggCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAAt                                                                       >  1:5226451/1‑93 (MQ=255)

                     gCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATg                                                                      <  2:1797939/93‑1 (MQ=255)

                     gCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATg                                                                      <  2:5250806/93‑1 (MQ=255)

                       aTAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGc                                                                    <  1:1482601/93‑1 (MQ=255)

                         aGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCaa                                                                  >  1:3934911/1‑93 (MQ=255)

                         aGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCaa                                                                  >  2:1193574/1‑93 (MQ=255)

                         aGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCAAATGGCaa                                                                  >  2:2551049/1‑93 (MQ=255)

                            acacGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCaaaaa                                                               <  1:1724962/93‑1 (MQ=255)

                             cacGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAAt                                                              >  2:4761701/1‑93 (MQ=255)

                               cGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGt                                                            <  1:6018413/93‑1 (MQ=255)

                                gCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGtt                                                           <  2:2161313/93‑1 (MQ=255)

                                 ccAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTc                                                          <  1:4255528/93‑1 (MQ=255)

                                  cAGACAACTCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCt                                                         >  1:6466921/1‑93 (MQ=255)

                                  cAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCt                                                         >  1:5843928/1‑93 (MQ=255)

                                  cAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCt                                                         >  2:2433251/1‑93 (MQ=255)

                                  cAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCt                                                         >  2:5525100/1‑93 (MQ=255)

                                    gACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTa                                                       >  2:318539/1‑93 (MQ=255)

                                    gACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTa                                                       <  2:1614994/93‑1 (MQ=255)

                                     aCAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTAc                                                      <  2:5855776/93‑1 (MQ=255)

                                        aCGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGa                                                   <  1:1227276/93‑1 (MQ=255)

                                        aCGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGa                                                   <  1:1132555/93‑1 (MQ=255)

                                         cgcgTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGaa                                                  <  2:5579730/93‑1 (MQ=255)

                                          gcgTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAg                                                 >  1:6170810/1‑93 (MQ=255)

                                           cgTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGc                                                <  2:3331261/93‑1 (MQ=255)

                                           cgTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGc                                                <  2:5618296/93‑1 (MQ=255)

                                              aGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCggg                                             <  2:2693683/93‑1 (MQ=255)

                                               gATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGt                                            >  1:3861315/1‑93 (MQ=255)

                                               gATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCCGGt                                            >  1:2775193/1‑93 (MQ=255)

                                               gATCATCAGCAGGGTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCGTACCGAAGCGGGt                                            >  2:3791121/1‑93 (MQ=255)

                                               aatcatcAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGt                                            <  1:3060786/92‑1 (MQ=255)

                                                atcatcAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGtt                                           <  2:3358335/93‑1 (MQ=255)

                                                 tcatcaGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGttt                                          <  2:3531114/93‑1 (MQ=255)

                                                 tcatcaGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGttt                                          <  2:805363/93‑1 (MQ=255)

                                                 tcatcaGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGttt                                          <  1:2087014/93‑1 (MQ=255)

                                                   atcaGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTg                                        >  1:3325534/1‑93 (MQ=255)

                                                   atcaGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTg                                        <  2:4188005/93‑1 (MQ=255)

                                                     cagcagGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGtt                                      >  2:4060962/1‑93 (MQ=255)

                                                      agcagGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTc                                     >  1:3605919/1‑93 (MQ=255)

                                                      agcagGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTc                                     <  2:684948/93‑1 (MQ=255)

                                                      agcagGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTc                                     <  2:2978951/93‑1 (MQ=255)

                                                       gcagGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTcc                                    >  1:5938002/1‑93 (MQ=255)

                                                       gcagGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTcc                                    <  1:6351600/93‑1 (MQ=255)

                                                       gcagGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTcc                                    >  1:1538262/1‑93 (MQ=255)

                                                         agGTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCcac                                  <  1:1354768/93‑1 (MQ=255)

                                                           gTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCacaa                                >  1:2794923/1‑93 (MQ=255)

                                                           gTTCAGGCCGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTCCCacaa                                <  1:5521530/93‑1 (MQ=255)

                                                                   cGACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGa                        <  1:2546604/93‑1 (MQ=255)

                                                                   cGACGTTCACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGCTCTTACCGCAGCGGGTTTTGTTCCACAACACCCAGa                        >  1:1005878/1‑93 (MQ=255)

                                                                    gACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGAc                       >  1:6160091/1‑93 (MQ=255)

                                                                    gACGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGAc                       >  2:2062159/1‑93 (MQ=255)

                                                                      cGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGt                                            >  2:5074774/1‑70 (MQ=255)

                                                                      cGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGt                                            <  1:5074774/70‑1 (MQ=255)

                                                                      cGTTGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGa                     <  1:1786364/93‑1 (MQ=255)

                                                                        ttGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAAc                   <  2:4436646/93‑1 (MQ=255)

                                                                        ttGACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAAc                   >  2:5497761/1‑93 (MQ=255)

                                                                          gACACCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGc                 <  1:455700/93‑1 (MQ=255)

                                                                           acacCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGcc                >  1:1581471/1‑93 (MQ=255)

                                                                           acacCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGcc                >  2:6575119/1‑93 (MQ=255)

                                                                            cacCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCa               >  2:4826140/1‑93 (MQ=255)

                                                                            cacCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCa               <  1:1438866/93‑1 (MQ=255)

                                                                             acCAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCAc              >  2:2318771/1‑93 (MQ=255)

                                                                              ccAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCAcc             >  1:2373892/1‑93 (MQ=255)

                                                                               cAGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCg            >  1:4329761/1‑93 (MQ=255)

                                                                               cAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGACCGCCACCg            >  1:5811357/1‑93 (MQ=255)

                                                                                aGATACTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGc           >  1:4702010/1‑93 (MQ=255)

                                                                                    tcTTTTGCTGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATc       <  1:3875496/92‑1 (MQ=255)

                                                                                      ttttGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCaa     <  1:3363391/93‑1 (MQ=255)

                                                                                         tGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAAcgc  >  2:1010222/1‑93 (MQ=255)

GCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGC  >  NC_000913/2235025‑2235206

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