BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,744	G→T	44.3%	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	44.3%	104.6 / 229.8	186	T85N (ACC→AAC)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base G (53/49); new base T (35/47); total (90/96)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

GCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCA  >  NC_000913/3379645‑3379840

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

 cAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGt                                                                                                 <  1:4273810/100‑1 (MQ=255)

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

     tAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTc                                                                                             <  2:1410638/100‑1 (MQ=255)

      aGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTcc                                                                                            <  2:349092/100‑1 (MQ=255)

      aGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTcc                                                                                            >  1:1044256/1‑100 (MQ=255)

       gCTCTTCATGCAGCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTccc                                                                                           >  2:2800799/1‑100 (MQ=255)

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

       gCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTccc                                                                                           <  1:388450/100‑1 (MQ=255)

        ctctTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCcg                                                                                          >  2:3432245/1‑100 (MQ=255)

            tCATGCACCCGCAGCATCAACCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTc                                                                                      <  2:3404570/100‑1 (MQ=255)

            tCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTc                                                                                      >  2:1020168/1‑100 (MQ=255)

            tCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTc                                                                                      >  2:1460966/1‑100 (MQ=255)

             cATGCATCCGCAGCATAAAACGGTGACAGTGCAGGCGCTGTTTCCGCTCTCCCGGCATGCCTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTcc                                                                                     <  1:577469/100‑1 (MQ=255)

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

               tGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAACAGGTCCATCAGCCAGTTTTTCCCGCGTCCTa                                                                                   >  1:150388/1‑100 (MQ=255)

               tGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCt                                                                                    >  1:4023688/1‑99 (MQ=255)

               tGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTa                                                                                   <  1:158110/100‑1 (MQ=255)

                 cACCCGCAGCATAAAACGATGAAAGTGCAAACGCTGTTTCCGCTCTCCCGGAAGACTTTGATAGAAAAGATACATCAACCAGATTTTACCGCGTCCTaaa                                                                                 >  1:905414/1‑98 (MQ=255)

                  aCCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTacac                                                                                <  1:2048779/100‑1 (MQ=255)

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

                    ccGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCg                                                                              >  1:2086849/1‑100 (MQ=255)

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

                      gcagcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGcc                                                                            <  1:4116672/100‑1 (MQ=255)

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

                        agcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGcccc                                                                          >  1:1110037/1‑100 (MQ=255)

                         gcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCca                                                                         <  2:604791/100‑1 (MQ=255)

                          caTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCcac                                                                        <  2:1851287/100‑1 (MQ=255)

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

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

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

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

                             aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACata                                                                     >  2:3896650/1‑100 (MQ=255)

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

                              aaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACatat                                                                    <  1:1020165/100‑1 (MQ=255)

                               aaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCATCCAGGTTTTCCCGCGTCCTACACCGCCCCACatata                                                                   <  1:2558202/100‑1 (MQ=255)

                                 cGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACAtatata                                                                 >  2:4014101/1‑100 (MQ=255)

                                 cGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACAtatata                                                                 <  2:1464692/100‑1 (MQ=255)

                                 cGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACAtatata                                                                 <  1:1831106/100‑1 (MQ=255)

                                  ggTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATaa                                                                >  2:3575672/1‑100 (MQ=255)

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

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

                                    ttaaaGTGCCGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGCTAGCAAAGGTCCTTCCGCCAGTTTTTCCCGCCTCCTACACCGCCCCACatctatccg                                                               >  1:225210/3‑96 (MQ=255)

                                    tgcaaGTGCAGGCGCTGTTTCCTCTCTCCCGGCAGTCTTTTATAGAAAAGGTCCATCATCCAGGTGTTCCCGCGTCCTACACCGCCCCACATATATAAGc                                                              <  1:3602152/97‑1 (MQ=255)

                                     gAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCATACAGTTTTTCCCGCGACCTACACCGCCCCACATATATAAGcc                                                             <  1:2225714/100‑1 (MQ=255)

                                      aaaGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCa                                                            <  2:2086853/100‑1 (MQ=255)

                                       aaGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCAc                                                           >  1:1256318/1‑100 (MQ=255)

                                        aGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACg                                                          >  1:122259/1‑100 (MQ=255)

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

                                              ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTgg                                                    >  2:3274285/1‑100 (MQ=255)

                                              ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTgg                                                    >  2:553289/1‑100 (MQ=255)

                                                cgcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGcn                                                  <  2:3328693/100‑2 (MQ=255)

                                                cgcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCg                                                  <  2:3822951/100‑1 (MQ=255)

                                                cgcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCg                                                  <  1:652165/100‑1 (MQ=255)

                                                 gcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGt                                                 >  2:1312834/1‑100 (MQ=255)

                                                  cTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTa                                                <  2:3857712/100‑1 (MQ=255)

                                                   tGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTAt                                               >  1:1951762/1‑100 (MQ=255)

                                                   tGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTAt                                               >  1:2119575/1‑100 (MQ=255)

                                                   tGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTAt                                               <  1:2731457/100‑1 (MQ=255)

                                                    gTTTCCGCTCTCCCGGCAGGCTTTGATAGATAAAGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCGCTGGCGTATg                                              >  1:2817387/1‑100 (MQ=255)

                                                        ccGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGCTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCt                                            >  1:3097475/1‑98 (MQ=255)

                                                         cGCTCTCCCGGCATGCTTTGATAGAAAAGTTCCATCATCCAGTTATTACCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                         <  1:3484603/100‑1 (MQ=255)

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

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

                                                         cGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                         >  1:916015/1‑100 (MQ=255)

                                                           ctctCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                       <  1:2656971/100‑1 (MQ=255)

                                                           ctctCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                       <  2:768444/100‑1 (MQ=255)

                                                             ctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTACTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTgtgt                                      >  1:1346535/1‑99 (MQ=255)

                                                             cacccGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTc                                     >  2:3240692/3‑100 (MQ=255)

                                                               cccGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCtt                                   >  1:329837/1‑100 (MQ=255)

                                                                ccGGGAGGCGTTGATAGAGCAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACATCGCCCCACATATAGAAGCCACGCACTGGCGTATGCTTtgtg                                       <  1:3932433/95‑1 (MQ=255)

                                                                 cGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcg                                 <  2:2570300/100‑1 (MQ=255)

                                                                  ggCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcgc                                <  2:2123282/100‑1 (MQ=255)

                                                                   gCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcgcg                               <  2:1604112/100‑1 (MQ=255)

                                                                        cTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAc                          >  2:28334/1‑100 (MQ=255)

                                                                         tttGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcc                         <  1:1922893/100‑1 (MQ=255)

                                                                          ttGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAccc                        >  2:2702329/1‑100 (MQ=255)

                                                                           tGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                       >  2:2087425/1‑100 (MQ=255)

                                                                            gATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCca                      >  2:2978606/1‑100 (MQ=255)

                                                                             attgAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGTCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                       <  1:1141460/95‑1 (MQ=255)

                                                                             aTAGAAAAGCTCCATCAGCCAGGTTTTCCCTCTCTCTACACCGCCCCACATATATACGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcac                     <  1:963654/100‑1 (MQ=255)

                                                                              tataaaaGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcaca                    <  1:1022325/97‑1 (MQ=255)

                                                                                gAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                  <  1:1727678/100‑1 (MQ=255)

                                                                                gAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                  >  1:785076/1‑100 (MQ=255)

                                                                                gAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                  >  1:3482296/1‑100 (MQ=255)

                                                                                  aaaGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCtt                <  1:2938570/100‑1 (MQ=255)

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

                                                                                  aaaGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCtt                <  2:1196741/100‑1 (MQ=255)

                                                                                    aGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCCTTGTGTCTTCGCGTTTACCCCACAGCTTAc              <  1:3269882/100‑1 (MQ=255)

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

                                                                                      gTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCg            <  1:2161285/100‑1 (MQ=255)

                                                                                       tCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAGGCCACGCACTGGCATATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGa           >  1:2375891/1‑100 (MQ=255)

                                                                                        ccATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc          >  2:3943847/1‑100 (MQ=255)

                                                                                         cATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAcc         >  2:2160163/1‑100 (MQ=255)

                                                                                         cATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAcc         >  2:2317945/1‑100 (MQ=255)

                                                                                          aTCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGATTTGTCTCTTCGCGTTTACCCCACAGCTTACCGAccc        >  1:2874279/1‑100 (MQ=255)

                                                                                            cAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgc      >  2:962661/1‑100 (MQ=255)

                                                                                             aGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgcg     >  2:1075167/1‑100 (MQ=255)

                                                                                               ccAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGcc   <  2:1996835/100‑1 (MQ=255)

                                                                                               ccAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGcc   <  1:3417788/100‑1 (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: `‑`