BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,744	G→T	82.1%	T85N (ACC→AAC)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,744	0	G	T	82.1%	362.0 / 82.7	168	T85N (ACC→AAC)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base G (15/15); new base T (72/65); total (88/80)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

Only 100 of 170 total aligned reads displayed.

GCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCA  >  NC_000913/3379645‑3379840

gCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAgg                                                                                                  <  1:699080/100‑1 (MQ=255)

  aGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGtt                                                                                                <  1:2141418/100‑1 (MQ=255)

   gTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGtttt                                                                                               >  2:2311081/1‑100 (MQ=255)

     tAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTc                                                                                             <  2:3280707/100‑1 (MQ=255)

      aGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCACTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTc                                                                                             <  2:2268408/99‑1 (MQ=255)

       gCTCTTCATGCACCCGCAGCATAAACCGGTGAAAGGGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTccc                                                                                           >  2:3984980/1‑100 (MQ=255)

       gCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTccc                                                                                           >  2:1524737/1‑100 (MQ=255)

        ctctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTtttt                                                                                              >  2:94797/1‑96 (MQ=255)

        ctctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTccc                                                                                           >  1:617425/1‑99 (MQ=255)

          ctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCcgcg                                                                                        <  2:1437616/100‑1 (MQ=255)

            tCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTc                                                                                      >  1:3722255/1‑100 (MQ=255)

             cATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTcc                                                                                     >  1:4049726/1‑100 (MQ=255)

             cATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTcc                                                                                     <  1:1718257/100‑1 (MQ=255)

              aTGCACCCGCCGCATAAAACGGCGAACGTGCAGGCGCTCTTTCCGCTCTCCCGGCAGGCTTTGATCGCAAAGGTCCATCAGCCAGCTCTTCCcgccccct                                                                                    >  2:62416/1‑95 (MQ=255)

               tGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCt                                                                                    >  2:3094441/1‑99 (MQ=255)

               tGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTa                                                                                   >  2:2619349/1‑100 (MQ=255)

                gCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTac                                                                                  <  1:2717304/100‑1 (MQ=255)

                   cccGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACAcc                                                                               >  2:2533544/1‑100 (MQ=255)

                     cGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGc                                                                             >  1:816815/1‑100 (MQ=255)

                      gcagcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGcc                                                                            >  1:637939/1‑100 (MQ=255)

                      gcagcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCGTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTacac                                                                                >  2:2001042/1‑96 (MQ=255)

                        agcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGcccc                                                                          <  1:432265/100‑1 (MQ=255)

                          caTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCcac                                                                        >  1:2065360/1‑100 (MQ=255)

                           aTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCcaca                                                                       >  1:2420445/1‑100 (MQ=255)

                           aTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCcaca                                                                       <  2:2133802/100‑1 (MQ=255)

                           aTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCcaca                                                                       >  2:368937/1‑100 (MQ=255)

                           aTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCcaca                                                                       <  1:3806008/100‑1 (MQ=255)

                            tAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAAGCTTTGATGGAAAAGGTCCATAAGCCAGGTTTTCCAGCGTCCTACAACGCCCcaca                                                                       >  1:398381/1‑99 (MQ=255)

                            tAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCGCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTCCCCGCGTCCTACACCGCCCCACat                                                                      <  2:15984/100‑1 (MQ=255)

                            tAAAACGGTGAAAGTGCAGGCGATGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACat                                                                      >  2:2412728/1‑100 (MQ=255)

                             aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACata                                                                     <  2:2993423/100‑1 (MQ=255)

                              aaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACAGCCCCACatat                                                                    >  1:1139745/1‑100 (MQ=255)

                              aaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACatat                                                                    >  1:561154/1‑100 (MQ=255)

                              aaaCGGTGAAAGTGCAGGCGCTGTTTCCGCCCTCCCGGCAGGCTTTGATAGAAAAGGACCATCAGCCCGTTATTCCCGCGTCCTACACCACCCCACatat                                                                    >  2:3457203/1‑100 (MQ=255)

                               aaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACatata                                                                   <  1:2711207/100‑1 (MQ=255)

                                aCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACatatat                                                                  <  2:3329396/100‑1 (MQ=255)

                                 cGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACAtatata                                                                 <  1:1780540/100‑1 (MQ=255)

                                 cGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACAtatata                                                                 >  2:3458710/1‑100 (MQ=255)

                                   gTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAg                                                               <  2:2101487/100‑1 (MQ=255)

                                   gTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAg                                                               <  1:1604258/100‑1 (MQ=255)

                                   gTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCACTTTTTCCCGCGTCCTACACCGCCCCACATATATAAg                                                               <  1:2942711/100‑1 (MQ=255)

                                   gTGAAAGTGCAGGCGCTGTTTCCGCTCGCCCGGCAGGCTTTGGTAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAg                                                               <  2:3303059/100‑1 (MQ=255)

                                       aaGTGCATGCTCTGTTTCCGCCCTCCCGGCAGGCTTTTATATAAAAGGTCCATCATCCAGTTTTTACCGCGTCCTACACCACCCCACATATATAAGCCAc                                                           <  2:794119/100‑1 (MQ=255)

                                         gTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGc                                                         <  2:724304/100‑1 (MQ=255)

                                         gTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGc                                                         >  1:3090992/1‑100 (MQ=255)

                                         gTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGc                                                         >  2:1796430/1‑100 (MQ=255)

                                              ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTgg                                                    >  2:875340/1‑100 (MQ=255)

                                              ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTgg                                                    <  1:3277939/100‑1 (MQ=255)

                                              ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTgg                                                    <  1:1180874/100‑1 (MQ=255)

                                                cgcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCg                                                  <  1:1552691/100‑1 (MQ=255)

                                                 gcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTCCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGt                                                 <  2:1606100/100‑1 (MQ=255)

                                                  cTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGt                                                 >  1:775687/1‑99 (MQ=255)

                                                  cTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTa                                                <  2:3982516/100‑1 (MQ=255)

                                                    gTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATg                                              <  2:134806/100‑1 (MQ=255)

                                                    gTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATg                                              <  2:3167240/100‑1 (MQ=255)

                                                      ttCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCt                                            >  2:590552/1‑100 (MQ=255)

                                                      ttCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACAGCCCCACATATATAAGCCACGAACTGGCGTATGCt                                            >  1:2478371/1‑100 (MQ=255)

                                                      ttCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTCGCGTATGCt                                            >  2:2833538/1‑100 (MQ=255)

                                                       tCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCt                                            >  2:3595317/1‑99 (MQ=255)

                                                        ccGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCttt                                          >  2:3930290/1‑100 (MQ=255)

                                                        ccGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCttt                                          <  1:222918/100‑1 (MQ=255)

                                                        ccGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCAACATATATAAGCCACGCACTGGCGTATGCttt                                          >  2:526800/1‑100 (MQ=255)

                                                         cgcgCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                         <  2:1152666/100‑1 (MQ=255)

                                                         cGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                         <  2:1240318/100‑1 (MQ=255)

                                                         cGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                         <  2:3346201/100‑1 (MQ=255)

                                                         cGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                         >  1:3465137/1‑100 (MQ=255)

                                                          gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                        <  1:849052/100‑1 (MQ=255)

                                                           ctctCCCGTCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                       <  2:2561248/100‑1 (MQ=255)

                                                            gctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                       <  2:3596792/98‑1 (MQ=255)

                                                             ctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTc                                     >  1:1657479/1‑100 (MQ=255)

                                                              tcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCt                                    <  2:2597868/100‑1 (MQ=255)

                                                                   gCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTGCACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcgcg                               >  1:3931407/1‑100 (MQ=255)

                                                                    cAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGt                              <  2:1042241/100‑1 (MQ=255)

                                                                      ggCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGttt                            <  2:2186612/100‑1 (MQ=255)

                                                                        cTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAc                          >  1:2304159/1‑100 (MQ=255)

                                                                         tttGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcc                         >  1:911945/1‑100 (MQ=255)

                                                                         tttGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcc                         <  1:3735115/100‑1 (MQ=255)

                                                                           tGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                       <  2:1765707/100‑1 (MQ=255)

                                                                           tGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                       >  2:1730221/1‑100 (MQ=255)

                                                                            gATAGTAAAGGTCCATCAGACAGTGTTCCCCGCGTCCTACACCACCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCca                      <  2:3589673/100‑1 (MQ=255)

                                                                            gATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGtagatcgga                      >  1:1983303/1‑92 (MQ=255)

                                                                              tAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcaca                    <  1:107637/100‑1 (MQ=255)

                                                                               aGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATCTATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAg                   >  2:3378263/1‑100 (MQ=255)

                                                                               aGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAg                   >  1:3314936/1‑100 (MQ=255)

                                                                                gAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                  >  1:3021159/1‑100 (MQ=255)

                                                                                gAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCACGTTTACCCCACAGc                  >  2:570135/1‑100 (MQ=255)

                                                                                  aaaGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCtt                >  1:1709680/1‑100 (MQ=255)

                                                                                     ggTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc             >  1:2515282/1‑100 (MQ=255)

                                                                                     ggTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc             <  1:3930320/100‑1 (MQ=255)

                                                                                     ggTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc             <  1:3248619/100‑1 (MQ=255)

                                                                                      gTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCg            >  1:2238698/1‑100 (MQ=255)

                                                                                        ccATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc          <  1:2055000/100‑1 (MQ=255)

                                                                                        ccATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc          >  1:2851283/1‑100 (MQ=255)

                                                                                        caaTCAGCCAGTTTTCCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc          <  2:166007/98‑1 (MQ=255)

                                                                                           tCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCACGTTTACCCCACAGCTTACCGACCcg       >  2:2651676/1‑100 (MQ=255)

                                                                                             aGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgcg     <  1:3563524/100‑1 (MQ=255)

                                                                                             aGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgcg     <  1:3734320/100‑1 (MQ=255)

                                                                                             aGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgcg     <  1:2981070/100‑1 (MQ=255)

                                                                                              gCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCgcgc    >  1:546612/1‑100 (MQ=255)

                                                                                                cAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCa  >  1:2419689/1‑100 (MQ=255)

GCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCA  >  NC_000913/3379645‑3379840

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