BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,530,925	T→C	16.6%	Q103Q (CAA→CAG)	zipA ←	FtsZ stabilizer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,530,925	0	T	C	16.6%	396.4 / 70.1	182	Q103Q (CAA→CAG)	zipA	FtsZ stabilizer
Reads supporting (aligned to +/- strand): ref base T (69/80); new base C (15/15); total (86/96)
Fisher's exact test for biased strand distribution p-value = 8.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 186 total aligned reads displayed.

CTGTTCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGATTCA  >  NC_000913/2530826‑2531022

cTGTTCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGt                                                                                                   >  2:3637266/1‑100 (MQ=255)

  gTTCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGtt                                                                                                  >  1:2479987/1‑99 (MQ=255)

   ttCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGa                                                                                                >  2:896504/1‑100 (MQ=255)

    tCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGAc                                                                                               <  2:1580496/100‑1 (MQ=255)

     cAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGAcg                                                                                              >  2:1779357/1‑100 (MQ=255)

     cAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGAcg                                                                                              <  2:182800/100‑1 (MQ=255)

     cAGGCTGAGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGAcg                                                                                              >  1:2174631/1‑100 (MQ=255)

       ggCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTAGCTGCGCTTCAGGCGGCTGCTGGACCGGGTGAcgcg                                                                                            >  1:1855496/1‑100 (MQ=255)

       ggCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGAcgcg                                                                                            >  1:1928328/1‑100 (MQ=255)

       ggCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGAcgcg                                                                                            <  1:2111226/100‑1 (MQ=255)

         cTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGGCCGGTTGACGCGGc                                                                                          >  2:491171/1‑100 (MQ=255)

             ggCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCt                                                                                         >  2:3386850/1‑97 (MQ=255)

              gCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTgcgc                                                                                     >  1:1829042/1‑100 (MQ=255)

               cTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCa                                                                                    <  1:3002764/100‑1 (MQ=255)

               cTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCa                                                                                    <  1:1107534/100‑1 (MQ=255)

                   tAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACg                                                                                >  1:3126797/1‑100 (MQ=255)

                   tAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACg                                                                                <  1:1908867/100‑1 (MQ=255)

                     ggcaggcTGGTTCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCATACg                                                                                <  2:170758/98‑1 (MQ=255)

                     ggcaggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGc                                                                               >  2:230431/1‑99 (MQ=255)

                     ggcaggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCCGACGCGGCTGCGCAGACGCa                                                                              >  2:936249/1‑100 (MQ=255)

                       caggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTCCCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCa                                                                              >  2:3859085/1‑98 (MQ=255)

                       caggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCAt                                                                             >  2:824398/1‑99 (MQ=255)

                          gctgctGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAgg                                                                         >  2:2984399/1‑100 (MQ=255)

                          gctgctGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAgg                                                                         <  2:3713833/100‑1 (MQ=255)

                                gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGTCTGACGCGACTGCGCAGACGCATAAGGCGGTTg                                                                   >  1:2257415/1‑100 (MQ=255)

                                gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGGTg                                                                   >  2:1086745/1‑100 (MQ=255)

                                gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCGTCAGGCGGCTGCTGGACCGGTAGACGCGGCTGCGCGGACGCATAAGGCGGGTg                                                                   >  2:2789309/1‑100 (MQ=255)

                                 cACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTgg                                                                  <  1:4041541/100‑1 (MQ=255)

                                  aCCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTAACGCGGCTGCGCAGACGCATAAGGCGGTTgg                                                                  >  2:2904224/1‑99 (MQ=255)

                                   ccGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTa                                                                >  2:2547057/1‑100 (MQ=255)

                                     ggCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACt                                                              >  2:1144155/1‑100 (MQ=255)

                                       cTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTgg                                                            >  2:1131788/1‑100 (MQ=255)

                                       cTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTgg                                                            <  1:3637290/100‑1 (MQ=255)

                                       cTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTgg                                                            <  2:3016797/100‑1 (MQ=255)

                                         gcgcTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGgtg                                                          >  1:3426719/1‑100 (MQ=255)

                                          cgcTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGgtgt                                                         <  1:1900652/100‑1 (MQ=255)

                                           gcTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGtt                                                        <  2:369947/100‑1 (MQ=255)

                                             cgatctcGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgc                                                      <  1:943670/94‑1 (MQ=255)

                                               gATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcgg                                                    <  1:1549610/100‑1 (MQ=255)

                                               gATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcgg                                                    <  2:1350211/100‑1 (MQ=255)

                                               gATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcgg                                                    <  1:1203428/100‑1 (MQ=255)

                                                aTGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggc                                                   <  1:1502844/100‑1 (MQ=255)

                                                aTGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggc                                                   >  1:2138812/1‑100 (MQ=255)

                                                 tGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggcg                                                  >  1:283770/1‑100 (MQ=255)

                                                  gCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCga                                                 >  1:2566025/1‑100 (MQ=255)

                                                   cGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgac                                                <  1:1713248/100‑1 (MQ=255)

                                                   cGGAGCATGTCGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgac                                                >  1:4000338/1‑100 (MQ=255)

                                                    ggAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacg                                               <  1:2606358/100‑1 (MQ=255)

                                                      aGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacgga                                             <  2:1041172/100‑1 (MQ=255)

                                                          tGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGag                                         >  2:3433560/1‑100 (MQ=255)

                                                           gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                        >  1:2925913/1‑100 (MQ=255)

                                                           gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                        >  1:3621143/1‑100 (MQ=255)

                                                           gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                        >  1:1493303/1‑100 (MQ=255)

                                                           gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                        <  2:247953/100‑1 (MQ=255)

                                                           gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGGTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                        <  1:4247/100‑1 (MQ=255)

                                                               cggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGcc                                    <  2:2731134/100‑1 (MQ=255)

                                                               cggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGcc                                    >  2:3254943/1‑100 (MQ=255)

                                                               cggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGcc                                    >  1:2396587/1‑100 (MQ=255)

                                                               cggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGcc                                    <  1:746357/100‑1 (MQ=255)

                                                                ggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCT                                   <  2:4000312‑M1/100‑1 (MQ=255)

                                                                ggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCT                                   >  1:1950312‑M1/1‑100 (MQ=255)

                                                                 gcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTC                                  >  1:1651670‑M1/1‑100 (MQ=255)

                                                                     tACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGC                              <  2:3091874‑M1/100‑1 (MQ=255)

                                                                     tACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGC                              <  1:3401974‑M1/100‑1 (MQ=255)

                                                                     tACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGC                              <  1:724753‑M1/100‑1 (MQ=255)

                                                                       ccTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTc                            >  1:816914/1‑100 (MQ=255)

                                                                        cTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTcc                           <  2:3771980/100‑1 (MQ=255)

                                                                        cTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTcc                           <  1:2774252/100‑1 (MQ=255)

                                                                        cTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTcc                           <  1:2547071/100‑1 (MQ=255)

                                                                            gcTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCccctgag                       >  2:2105237‑M1/1‑93 (MQ=255)

                                                                            gcTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAg                       <  1:1512722/100‑1 (MQ=255)

                                                                             cTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGc                      >  1:3073699/1‑100 (MQ=255)

                                                                             cTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGc                      <  1:536826/100‑1 (MQ=255)

                                                                             cTTCAGGCGGCTGCTAGACTGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGc                      <  1:4068338/100‑1 (MQ=255)

                                                                                 aGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTa                  <  2:1456381/100‑1 (MQ=255)

                                                                                  ggcggcTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAg                 <  2:2232467/100‑1 (MQ=255)

                                                                                  ggcggcTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAg                 >  2:4010765/1‑100 (MQ=255)

                                                                                  ggcggcTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAg                 >  2:479548/1‑100 (MQ=255)

                                                                                  ggcggcTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAg                 <  1:3035684/100‑1 (MQ=255)

                                                                                  ggcggcTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGGTGGGACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGtt                   >  2:4052623/1‑98 (MQ=255)

                                                                                  ggcggcTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAg                 >  1:2251758/1‑100 (MQ=255)

                                                                                      gctgctGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCgg             >  1:1779981/1‑100 (MQ=255)

                                                                                          ctgGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGgca         >  1:3242459‑M2/1‑97 (MQ=255)

                                                                                           tgGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcat        <  2:2366155‑M1/100‑23 (MQ=255)

                                                                                           tgGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcat        <  1:1286842‑M1/100‑23 (MQ=255)

                                                                                            ggACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatg       <  1:3066580‑M1/100‑24 (MQ=255)

                                                                                            ggACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGgcatg       <  2:1204382‑M2/100‑6 (MQ=255)

                                                                                             gACCGTGTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTGGCGGCGACGGACGCGCAGCCTCATGCTCCTGAGCGTTAGCCGGGgcatga      >  2:709273‑M2/1‑94 (MQ=255)

                                                                                             gACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatga      <  2:58642‑M1/100‑25 (MQ=255)

                                                                                              aCCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatgag     <  1:1363558‑M1/100‑26 (MQ=255)

                                                                                              aCCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatgag     <  2:553754‑M1/100‑26 (MQ=255)

                                                                                              aCCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatgag     <  1:238314‑M1/100‑26 (MQ=255)

                                                                                               ccGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatga      <  2:2042123‑M1/98‑25 (MQ=255)

                                                                                                cGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatgaggc   <  1:502471‑M1/100‑28 (MQ=255)

                                                                                                cGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatgaggc   <  1:1029624‑M1/100‑28 (MQ=255)

                                                                                                cGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGATTc   >  2:2518199/1‑100 (MQ=255)

                                                                                                cGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGATTc   <  1:2409764/100‑1 (MQ=255)

                                                                                                 ggTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCcccggctaacgctcaggagcatgaggct  <  1:3615159‑M1/100‑29 (MQ=255)

                                                                                                 ggTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGATTCa  >  1:3266857/1‑100 (MQ=255)

                                                                                                 ggTTGACGCGGCTGCGCAGACGCATAAGGCGGGAGGTACTGGTGTTGCGGCGACGGACGAGCATCCTCATGCTCCTGAGCGTTAGCCGGGGCATGATTCa  >  2:1867856/1‑100 (MQ=255)

CTGTTCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGATTCA  >  NC_000913/2530826‑2531022

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: `‑`