BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,539	T→C	intergenic (‑250/+21)	b3460 ← / ← rpoH	LivJ/sigma(32) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,539	0	T	C	99.5%	1168.0 / ‑5.5	374	intergenic (‑250/+21)	b3460/rpoH	LivJ/sigma(32) factor
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (221/151); total (222/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 387 total aligned reads displayed.

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

aCAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtc                                                                                         <  2:4981205/93‑1 (MQ=255)

aCAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtc                                                                                         >  1:234119/1‑93 (MQ=255)

 cAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctct                                                                                        >  1:362115/1‑93 (MQ=255)

  aaTTTGCTTGGAGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       <  2:6592698/93‑1 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       <  1:5505804/93‑1 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       >  2:2202190/1‑93 (MQ=255)

   aTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGTAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGTctctct                                                                                      <  2:3715147/93‑1 (MQ=255)

   aTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctctct                                                                                      >  1:278765/1‑93 (MQ=255)

     cTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGc                                                                                    >  2:5730961/1‑93 (MQ=255)

      tGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCt                                                                                   <  1:5360458/93‑1 (MQ=255)

         gTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTaa                                                                                >  1:7247577/1‑93 (MQ=255)

         gTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTaa                                                                                >  1:2977329/1‑93 (MQ=255)

         gTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTaa                                                                                <  1:1612586/93‑1 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               >  2:2301811/1‑93 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               >  1:6382243/1‑93 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               >  2:4740062/1‑93 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               >  2:7058224/1‑93 (MQ=255)

           ggggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAtg                                                                              >  2:1211299/1‑92 (MQ=255)

           ggggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATa                                                                              <  2:1655720/93‑1 (MQ=255)

            gggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAg                                                                             <  2:1472503/93‑1 (MQ=255)

             ggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAACCCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGc                                                                            <  2:1553026/93‑1 (MQ=255)

              gATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCg                                                                           >  2:6915718/1‑93 (MQ=255)

              gATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCg                                                                           >  1:5189024/1‑93 (MQ=255)

               aTTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCgg                                                                          >  2:6449745/1‑93 (MQ=255)

                ttGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGa                                                                         <  2:3056719/93‑1 (MQ=255)

                 tGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaa                                                                        <  2:1515293/93‑1 (MQ=255)

                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      <  1:454758/93‑1 (MQ=255)

                    cGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAtt                                                                     >  1:1172779/1‑93 (MQ=255)

                     gTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTa                                                                    <  1:875050/93‑1 (MQ=255)

                     gTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTa                                                                    <  1:1838560/93‑1 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   <  1:484291/93‑1 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   >  2:2025062/1‑93 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   >  2:1886673/1‑93 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   >  1:3794679/1‑93 (MQ=255)

                        ttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGc                                                                 <  2:6237147/93‑1 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            >  2:1056940/1‑93 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            >  2:7485652/1‑93 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            <  1:5565123/93‑1 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           >  2:7123778/1‑93 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           >  2:6330911/1‑93 (MQ=255)

                               aGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATg                                                          <  1:4490744/93‑1 (MQ=255)

                               aGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATg                                                          <  2:6017561/93‑1 (MQ=255)

                                ggggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATgg                                                         >  2:83714/1‑93 (MQ=255)

                                ggggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATgg                                                         <  1:4239316/93‑1 (MQ=255)

                                  ggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcg                                                       >  1:6108170/1‑92 (MQ=255)

                                  ggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgca                                                       <  1:7436549/93‑1 (MQ=255)

                                    aaaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagc                                                     <  1:1023433/93‑1 (MQ=255)

                                    aaaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagc                                                     <  1:334464/93‑1 (MQ=255)

                                        tGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGc                                                 >  2:6163902/1‑93 (MQ=255)

                                         gATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCa                                                >  2:4205410/1‑93 (MQ=255)

                                          attattACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCaa                                               >  1:6441371/1‑93 (MQ=255)

                                           ttattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAt                                              >  2:7218120/1‑93 (MQ=255)

                                           ttattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAt                                              >  1:7188419/1‑93 (MQ=255)

                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           >  1:5885238/1‑93 (MQ=255)

                                               taCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttttt                                          >  1:2205869/1‑93 (MQ=255)

                                               taCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttttt                                          >  1:2606984/1‑93 (MQ=255)

                                                aCAGAGGCCCAAAAAACAAAAACCCCGGCCTCTGATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttttt                                         <  1:3950013/93‑1 (MQ=255)

                                                 catagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATTGTATCACGCATTTTTTTc                                        >  1:7023773/4‑93 (MQ=255)

                                                      gcgcAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcgc                                   <  1:5673754/90‑1 (MQ=255)

                                                      gCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcgc                                   >  2:7743358/1‑93 (MQ=255)

                                                       cccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCCATTTTTTCATcgcg                                  >  2:6063589/1‑93 (MQ=255)

                                                         cAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGtt                                >  1:4733762/1‑93 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              >  2:2950251/1‑93 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCCCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              >  2:7615223/1‑93 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             >  2:5479158/1‑93 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             <  1:2520798/93‑1 (MQ=255)

                                                             aaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            >  2:5314670/1‑93 (MQ=255)

                                                             aaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            >  2:4566675/1‑93 (MQ=255)

                                                             aaacaaaaacCCCGGACCCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            >  1:6182281/1‑93 (MQ=255)

                                                              aacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtttt                           >  1:3293797/1‑93 (MQ=255)

                                                              aacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtttt                           >  1:6753501/1‑93 (MQ=255)

                                                               acaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTc                          <  1:4192225/93‑1 (MQ=255)

                                                                caaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTcc                         <  2:2479303/93‑1 (MQ=255)

                                                                 aaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCa                        >  2:4680333/1‑93 (MQ=255)

                                                                   aaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGc                      <  1:4095349/93‑1 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     <  1:1391517/93‑1 (MQ=255)

                                                                     acCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTg                    <  2:5085261/93‑1 (MQ=255)

                                                                      ccccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTgg                   <  2:7880686/93‑1 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  >  1:6552198/1‑93 (MQ=255)

                                                                        ccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCg                 >  2:516437/1‑93 (MQ=255)

                                                                        ccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCg                 >  2:4502555/1‑93 (MQ=255)

                                                                        ccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCg                 >  1:7466487/1‑93 (MQ=255)

                                                                        ccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCg                 >  1:3137672/1‑93 (MQ=255)

                                                                         cGGACTCTCTTCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGt                >  2:3148728/1‑93 (MQ=255)

                                                                         cGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGt                >  2:1404214/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               >  2:1784999/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               >  1:7734276/1‑93 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              >  1:6875866/1‑93 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              >  1:5618298/1‑93 (MQ=255)

                                                                            aCTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTaca             <  2:3523928/93‑1 (MQ=255)

                                                                            aCTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTaca             <  2:330595/93‑1 (MQ=255)

                                                                               ctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGc          >  1:5467254/1‑93 (MQ=255)

                                                                               ctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGc          >  1:509293/1‑93 (MQ=255)

                                                                                 cATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTc        <  1:1368540/93‑1 (MQ=255)

                                                                                 cATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTc        <  1:730009/93‑1 (MQ=255)

                                                                                  aTCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCa       >  1:1764280/1‑93 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACCCGCTCAg      >  2:3918446/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:2056387/1‑93 (MQ=255)

                                                                                       gggCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGa  >  1:7932141/1‑93 (MQ=255)

                                                                                       gggCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGa  >  1:7157336/1‑93 (MQ=255)

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

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