BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,742	A→C	W86G (TGG→GGG)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,742	0	A	C	99.2%	468.1 / ‑4.1	127	W86G (TGG→GGG)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base C (66/60); total (66/61)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 128 total aligned reads displayed.

AGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGC  >  NC_000913/3379643‑3379838

aGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCa                                                                                                  <  1:2089935/100‑1 (MQ=255)

  gCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCgg                                                                                                <  1:2870051/100‑1 (MQ=255)

     gTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGttt                                                                                             <  2:2641646/100‑1 (MQ=255)

      ttAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGttt                                                                                             >  2:1728740/1‑99 (MQ=255)

       tAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTc                                                                                           >  2:3300175/1‑100 (MQ=255)

       tAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTc                                                                                           >  1:3539284/1‑100 (MQ=255)

        aGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTcc                                                                                          <  2:1395360/100‑1 (MQ=255)

          ctctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcg                                                                                        <  1:1686825/100‑1 (MQ=255)

          ctctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcg                                                                                        >  2:1162600/1‑100 (MQ=255)

          ctctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcg                                                                                        >  2:281591/1‑100 (MQ=255)

            ctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcgcg                                                                                      >  1:2423790/1‑100 (MQ=255)

            ctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcgcg                                                                                      <  1:907514/100‑1 (MQ=255)

            ctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcgcg                                                                                      <  1:156872/100‑1 (MQ=255)

            ctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcgcg                                                                                      <  1:3342990/100‑1 (MQ=255)

            ccgatctcACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCcgcg                                                                                      <  1:2360236/93‑1 (MQ=255)

              tCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTc                                                                                    >  2:2155515/1‑100 (MQ=255)

                 tGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTa                                                                                 >  2:744178/1‑100 (MQ=255)

                   cACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTaca                                                                               >  2:1646850/1‑100 (MQ=255)

                   cACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTaca                                                                               >  2:1722950/1‑100 (MQ=255)

                    aCCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTacac                                                                              >  2:2238769/1‑100 (MQ=255)

                     cccGCAGCATAAAACGGTGAAAGTGTAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACAcc                                                                             >  2:2013895/1‑100 (MQ=255)

                     cccGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGCTTTCCCGCGCCCTACAcc                                                                             >  2:3374394/1‑100 (MQ=255)

                       cGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGc                                                                           >  2:832469/1‑100 (MQ=255)

                       cGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGc                                                                           >  2:823068/1‑100 (MQ=255)

                             aTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCcaca                                                                     >  2:3322005/1‑100 (MQ=255)

                              tAAAACGGTGACAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCACGGTTTTACCGCGTCCTACACCGCCCCACat                                                                    <  1:177553/100‑1 (MQ=255)

                               aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACata                                                                   >  2:2577135/1‑100 (MQ=255)

                               aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACata                                                                   >  1:265726/1‑100 (MQ=255)

                               aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACata                                                                   <  2:159440/100‑1 (MQ=255)

                               aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGCCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACata                                                                   <  1:1589496/100‑1 (MQ=255)

                                 aaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACatata                                                                 <  2:3376045/100‑1 (MQ=255)

                                  aCGGTGCAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACatata                                                                 <  1:3305943/99‑1 (MQ=255)

                                  aCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACatatat                                                                <  2:689069/100‑1 (MQ=255)

                                    ggTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATaa                                                              >  1:221555/1‑100 (MQ=255)

                                    ggTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATaa                                                              <  2:1230083/100‑1 (MQ=255)

                                      tGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGc                                                            >  1:2848784/1‑100 (MQ=255)

                                      tGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGc                                                            >  2:640240/1‑100 (MQ=255)

                                       gAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGcc                                                           <  1:931987/100‑1 (MQ=255)

                                        aaaGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCa                                                          <  1:1674595/100‑1 (MQ=255)

                                         aaGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCAc                                                         >  2:849473/1‑100 (MQ=255)

                                         aaGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCAc                                                         >  2:2126460/1‑100 (MQ=255)

                                          aGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACg                                                        >  1:782038/1‑100 (MQ=255)

                                             gCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCAc                                                     >  1:1906852/1‑100 (MQ=255)

                                               aGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTg                                                   <  1:2577117/100‑1 (MQ=255)

                                                ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTgg                                                  <  1:1646844/100‑1 (MQ=255)

                                                ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTgg                                                  >  1:2787049/1‑100 (MQ=255)

                                                  cgcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCg                                                <  2:741525/100‑1 (MQ=255)

                                                   gcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGt                                               <  1:2391745/100‑1 (MQ=255)

                                                   gcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGt                                               <  1:241144/100‑1 (MQ=255)

                                                    cTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTa                                              >  1:2209976/1‑100 (MQ=255)

                                                     tGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTAt                                             <  1:1350863/100‑1 (MQ=255)

                                                      gTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATg                                            <  1:1162596/100‑1 (MQ=255)

                                                      gTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATg                                            <  2:297713/100‑1 (MQ=255)

                                                        ttCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCt                                          <  1:1777486/100‑1 (MQ=255)

                                                          ccGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCttt                                        <  1:3454283/100‑1 (MQ=255)

                                                           cGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                       <  1:2776038/100‑1 (MQ=255)

                                                            gCTCTCCCGGCAGGCTTTGATGGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                      <  1:264316/100‑1 (MQ=255)

                                                            gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                      <  1:640237/100‑1 (MQ=255)

                                                            gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                      <  1:2779758/100‑1 (MQ=255)

                                                            gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                      >  2:1938368/1‑100 (MQ=255)

                                                               ctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTc                                   >  1:769029/1‑100 (MQ=255)

                                                                 cccGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCtt                                 >  2:1988920/1‑100 (MQ=255)

                                                                   cGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcg                               <  2:1100236/100‑1 (MQ=255)

                                                                   cGGCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcg                               <  2:2885884/100‑1 (MQ=255)

                                                                    ggCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcgc                              >  1:1167922/1‑100 (MQ=255)

                                                                     gCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcgcg                             <  2:3539313/100‑1 (MQ=255)

                                                                     gCAGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcgcg                             <  2:1831150/100‑1 (MQ=255)

                                                                       aGGCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGtt                           >  1:2115075/1‑100 (MQ=255)

                                                                        ggCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGttt                          >  1:3395114/1‑100 (MQ=255)

                                                                        ggCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGACACGCACTGGCGTATGCTTTGTGTCTTCGCGttt                          >  1:1868633/1‑100 (MQ=255)

                                                                         gCTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCAGCGTCCTACACAGCCCCACGTATATAAGCCTCGCCCTGGCCTATGCTTTGTGTCTTCGCGTTTa                         >  1:1569416/1‑100 (MQ=255)

                                                                          cTTTGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAc                        >  2:2804732/1‑100 (MQ=255)

                                                                             tGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                     >  2:1279500/1‑100 (MQ=255)

                                                                             tGATAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                     <  2:2744462/100‑1 (MQ=255)

                                                                             tGATAGAAAAGGTCCATCAGCCCGGTTTTCCAGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                     >  2:2089601/1‑100 (MQ=255)

                                                                               aTAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcac                   <  2:1623746/100‑1 (MQ=255)

                                                                                tAGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcaca                  >  1:3213181/1‑100 (MQ=255)

                                                                                 aGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAg                 <  1:3508652/100‑1 (MQ=255)

                                                                                 aGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAg                 <  1:832465/100‑1 (MQ=255)

                                                                                 aGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAg                 >  2:1192234/1‑100 (MQ=255)

                                                                                 aGAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTGGTGTCTTCGCGTTTACCCCACAg                 >  1:861358/1‑100 (MQ=255)

                                                                                  gAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                >  1:1978853/1‑100 (MQ=255)

                                                                                  gAAAAGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                >  2:374049/1‑100 (MQ=255)

                                                                                  gAAAAGGTCCATCAGCCCGGTTTTCCCGCGGCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                <  1:2105978/100‑1 (MQ=255)

                                                                                     aaGGTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTa             <  2:210751/100‑1 (MQ=255)

                                                                                       ggTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc           >  1:1500250/1‑100 (MQ=255)

                                                                                       ggTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc           <  1:2854693/100‑1 (MQ=255)

                                                                                       ggTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc           >  1:2399675/1‑100 (MQ=255)

                                                                                        gTCCATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCg          <  1:2581418/100‑1 (MQ=255)

                                                                                        gTCCATCAGCCCGGTTTTCCAGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCg          <  1:2089588/100‑1 (MQ=255)

                                                                                          ccATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc        <  2:221556/100‑1 (MQ=255)

                                                                                          ccATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc        <  2:2065087/100‑1 (MQ=255)

                                                                                          ccATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc        >  1:1399188/1‑100 (MQ=255)

                                                                                           cATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAcc       <  2:446971/100‑1 (MQ=255)

                                                                                           cATCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAcc       <  1:1823536/100‑1 (MQ=255)

                                                                                            ntCAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATAGAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAccc      >  1:1938670/2‑100 (MQ=255)

                                                                                              cAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgc    >  2:2501442/1‑100 (MQ=255)

                                                                                              cAGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcg     <  2:2439749/99‑1 (MQ=255)

                                                                                               aGCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcacg   <  2:3395140/100‑1 (MQ=255)

                                                                                                gCCCGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCgcgc  >  2:2442085/1‑100 (MQ=255)

AGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGC  >  NC_000913/3379643‑3379838

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