BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,049,989	A→C	16.0%	intergenic (+237/+144)	polA → / ← yihA	5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease/cell division GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,989	0	A	C	16.0%	388.7 / 71.0	156	intergenic (+237/+144)	polA/yihA	5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease/cell division GTP‑binding protein
Reads supporting (aligned to +/- strand): ref base A (71/60); new base C (16/9); total (87/69)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

Only 100 of 158 total aligned reads displayed.

GTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTAC  >  NC_000913/4049894‑4050088

gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                                 >  1:1121731/1‑100 (MQ=255)

gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                                 <  2:397766/100‑1 (MQ=255)

gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                                 >  1:221835/1‑100 (MQ=255)

   gcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                              >  2:3384123/1‑100 (MQ=255)

    cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                             >  2:362501/1‑100 (MQ=255)

    cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                             >  2:2313113/1‑100 (MQ=255)

      tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                           <  1:3229150/100‑1 (MQ=255)

      tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTGCACGTAATCGGATTTGGCTGAATATTGTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                           <  1:1443248/100‑1 (MQ=255)

      naGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc                                                                                           >  1:339661/2‑100 (MQ=255)

       aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCg                                                                                          <  1:17792/100‑1 (MQ=255)

        gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCg                                                                                          <  1:3093174/99‑1 (MQ=255)

        gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGt                                                                                         >  1:1694450/1‑100 (MQ=255)

         ggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCggc                                                                                        <  2:742718/100‑3 (MQ=255)

         ggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCggc                                                                                        <  2:1907963/100‑3 (MQ=255)

             cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                    <  2:2879919/100‑1 (MQ=255)

              agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCggctaaaa                                                                                   <  2:3094352/100‑8 (MQ=38)

              agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTa                                                                                   <  1:1060204/100‑1 (MQ=255)

               gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                                  >  1:1810619/1‑100 (MQ=255)

               gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                                  <  1:1551568/100‑1 (MQ=255)

                 ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAtt                                                                                 >  2:2405886/1‑99 (MQ=255)

                 ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                                <  2:2712830/100‑1 (MQ=255)

                 ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                                >  2:1390409/1‑100 (MQ=255)

                  gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                               >  1:2327317/1‑100 (MQ=255)

                   tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                              <  2:1857210/100‑1 (MQ=255)

                   tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                              >  1:224053/1‑100 (MQ=255)

                   tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                              >  2:488066/1‑100 (MQ=255)

                    aaGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc                                                                             <  2:353293/100‑1 (MQ=255)

                     aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCg                                                                            <  2:774327/100‑1 (MQ=255)

                       aTGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaa                                                                          >  1:701186/1‑100 (MQ=255)

                       aTGTTCTATCTTTCAGAACTTTTACTTCACGTAATAGGATTTGGCTGAATATTTTAGTCGCCCCAGTCAGTAATGACTGGGGCGTCTTTTATTAGGc                                                                             >  1:2844235/1‑97 (MQ=255)

                        tGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTGTATTGGGCgaaa                                                                         >  1:1707012/1‑100 (MQ=255)

                           tCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaa                                                                      >  2:881168/1‑100 (MQ=255)

                             tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                    >  2:1133056/1‑100 (MQ=255)

                              aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaag                                                                   <  2:850942/100‑1 (MQ=255)

                              aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaag                                                                   <  1:881165/100‑1 (MQ=255)

                                cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                                 <  1:2639467/100‑1 (MQ=255)

                                cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCACCAGTCAGTAATGAATGGGAAGTTTTTTCTTGGGAGcaagacaaga                                                                  >  2:281169/1‑99 (MQ=255)

                                 tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaaga                                                                  >  2:1970601/1‑98 (MQ=255)

                                 tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                                >  1:433424/1‑100 (MQ=255)

                                 tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                                <  2:1351250/100‑1 (MQ=255)

                                   tCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCg                                                              <  2:1083565/100‑1 (MQ=255)

                                      gACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGAGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                                           >  2:1385737/1‑100 (MQ=255)

                                         cTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGc                                                        <  1:2566549/100‑1 (MQ=255)

                                          ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGcc                                                       <  1:2125999/100‑1 (MQ=255)

                                               cTTCGCGTAATCGGATTTGGCTGCATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCTAAAGAAAAGATCCGTAATGCCTGATg                                                  <  1:1043768/100‑1 (MQ=255)

                                               cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                                  >  2:2416180/1‑100 (MQ=255)

                                               cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                                  <  2:1638385/100‑1 (MQ=255)

                                               cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                                  >  1:503225/1‑100 (MQ=255)

                                                ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgc                                                 <  2:3176508/100‑1 (MQ=255)

                                                 tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                <  2:1121734/100‑1 (MQ=255)

                                                 tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                >  2:1630996/1‑100 (MQ=255)

                                                 tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                >  1:1493420/1‑100 (MQ=255)

                                                   aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                              >  2:2519393/1‑100 (MQ=255)

                                                   aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                              >  2:2804506/1‑100 (MQ=255)

                                                      tAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                           >  1:1932508/1‑100 (MQ=255)

                                                        aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGtt                                         <  2:2896269/100‑1 (MQ=255)

                                                         tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTTTTTTTTTTTAGAAAAGATCCGTAATGCCTGATGCGCTATGttt                                        >  1:1384194/1‑100 (MQ=255)

                                                          cGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTa                                       <  2:1143804/100‑1 (MQ=255)

                                                          cGGATTTAGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTa                                       >  2:1776242/1‑100 (MQ=255)

                                                           ggATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTTGTTTTTTTAAAAGATATGTAATGCCTGATATGCTATTTTTAt                                      >  1:273561/1‑100 (MQ=255)

                                                           ggATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTAt                                      >  2:1321746/1‑100 (MQ=255)

                                                            gATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                     <  1:526385/100‑1 (MQ=255)

                                                             aTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAGTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                    <  2:1114647/100‑1 (MQ=255)

                                                                 ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                                <  2:1932521/100‑1 (MQ=255)

                                                                  gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCa                               <  2:1621565/100‑1 (MQ=255)

                                                                      aaTATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACgg                           >  2:2817767/1‑100 (MQ=255)

                                                                         ataaaaaacgccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAg                        >  2:3015215/9‑100 (MQ=37)

                                                                          ttttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGa                       >  1:1948320/1‑100 (MQ=255)

                                                                               gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGt                  <  1:2723628/100‑1 (MQ=255)

                                                                                ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa                 <  2:1001958/100‑1 (MQ=255)

                                                                                 cgccCCAGTCATTACTGACTGGGGCGTTTTTTTTTTTTTTTTTTAAAAGAACAGTAATACCTGAAAAGCTATGTTTATCAGGCCAACACTAGAACTg                   >  1:1491747/1‑97 (MQ=2)

                                                                                 cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTTTTTTTAATTTAATATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa                >  1:1789783/1‑100 (MQ=255)

                                                                                 cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa                >  2:2125031/1‑100 (MQ=255)

                                                                                 cgccCCAGTCAGTAACGACTGGGGCGTTTTTTATTTGTCGATAGCAATGTCCCGCGACGCCTGATTCCGTACGATTATCAgcc                                 >  1:3337537/1‑81 (MQ=255)

                                                                                  gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGTGCTATGTTTATCAGGCCAACGGTAGAATTGTaa                >  1:2126023/1‑99 (MQ=255)

                                                                                  gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt               >  2:2847084/1‑100 (MQ=255)

                                                                                  gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt               >  2:3043798/1‑100 (MQ=255)

                                                                                  gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt               >  2:1701110/1‑100 (MQ=255)

                                                                                  gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt               >  2:1882994/1‑100 (MQ=255)

                                                                                  gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt               >  2:852135/1‑100 (MQ=255)

                                                                                  gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGTGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt               >  1:1839928/1‑100 (MQ=255)

                                                                                  gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt               >  2:488330/1‑100 (MQ=255)

                                                                                   ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc              >  2:2994899/1‑100 (MQ=255)

                                                                                   ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc              >  2:2697828/1‑100 (MQ=255)

                                                                                   ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc              <  2:2923493/100‑1 (MQ=255)

                                                                                   ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc              >  2:710607/1‑100 (MQ=255)

                                                                                   ccccAGTCAGTAATGACTGGGGCGTTTTTTATTAGGCGAAAGAAAAGATCAGTAATGACTGATGCGCTATGTTTATCAGTACAACGCTAGACTTGTAATc              >  1:1166712/1‑100 (MQ=255)

                                                                                    cccAGTCAGTAATGACTGGGGCGTTTTTTTTTTTTTTTTTTAAAAGAAAAGTAATGCCTGATGCGCTATGTTTATCAGGCTAACGGTAGAATTGttttct             >  1:1295365/1‑95 (MQ=255)

                                                                                     ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa            >  1:2874770/1‑100 (MQ=255)

                                                                                     ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa            >  1:1950814/1‑100 (MQ=255)

                                                                                     ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa            <  2:2440389/100‑1 (MQ=255)

                                                                                      cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAt           <  2:1948333/100‑1 (MQ=255)

                                                                                      cagtcagtAACGATTTGGGCTTTTTTTATTGGGCAAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAATGGTAGAATTGTAATCTAt           <  1:3228454/100‑1 (MQ=255)

                                                                                       agtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt          >  2:267984/1‑100 (MQ=255)

                                                                                       agtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt          >  2:2647587/1‑100 (MQ=255)

                                                                                         tcagtAATGACTGGCGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTAGCAGGCCAACGGTAGAATTGTAATCTAt           <  1:1835500/97‑1 (MQ=255)

                                                                                          cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa       <  2:320786/100‑1 (MQ=255)

                                                                                           agtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt      <  2:255459/100‑1 (MQ=255)

                                                                                              aaTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTa   <  1:992945/100‑1 (MQ=255)

                                                                                               aTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTAc  >  1:1569681/1‑100 (MQ=255)

GTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTAC  >  NC_000913/4049894‑4050088

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