BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,958,357	G→A	S698S (AGC→AGT)	recC ←	RecC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,958,357	0	G	A	99.1%	1151.2 / ‑5.8	327	S698S (AGC→AGT)	recC	RecC
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (169/155); total (170/157)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

ATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGA  >  NC_000913/2958267‑2958446

 taGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtc                                                                                        >  1:9190996/1‑93 (MQ=255)

 taGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGCCTACGGTCGCCACGCTTCGGTTTCTGACtc                                                                                        >  1:8649708/1‑93 (MQ=255)

 taGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAAAAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCCGTTTCTGACtc                                                                                        >  2:7149450/1‑93 (MQ=255)

   gagTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcat                                                                                      >  2:6301078/1‑93 (MQ=255)

   gagTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGATACGGTCGCCACGCTTCGGTTTCTGACtcat                                                                                      >  2:3665557/1‑93 (MQ=255)

     gTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcatca                                                                                    >  2:4940465/1‑93 (MQ=255)

     gTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcatca                                                                                    <  1:2285292/93‑1 (MQ=255)

     gTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcatca                                                                                    >  1:7931417/1‑93 (MQ=255)

     gTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcatca                                                                                    <  1:5903237/93‑1 (MQ=255)

         ttGCTGCGCGGAAAGTAACGCTTCCACGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTc                                                                                <  1:633238/93‑1 (MQ=255)

            ctgcGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCaaa                                                                             <  2:9904454/93‑1 (MQ=255)

            ctgcGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCaaa                                                                             <  2:8370627/93‑1 (MQ=255)

                gcgGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGccc                                                                         >  1:6955189/1‑93 (MQ=255)

                  ggAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCaa                                                                       >  1:5835407/1‑93 (MQ=255)

                  ggAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCaa                                                                       <  2:5366101/93‑1 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      <  1:3808942/93‑1 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      <  1:6451638/93‑1 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      <  1:9313601/93‑1 (MQ=255)

                     aaTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATgg                                                                    >  2:1816834/1‑93 (MQ=255)

                     aaTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATgg                                                                    <  1:1420868/93‑1 (MQ=255)

                      aTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGgc                                                                   >  1:4804270/1‑93 (MQ=255)

                      aTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGgc                                                                   >  1:2492356/1‑93 (MQ=255)

                      aTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGgc                                                                   >  2:9397971/1‑93 (MQ=255)

                       ttAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGgcg                                                                  >  2:2678187/1‑93 (MQ=255)

                       ttAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGgcg                                                                  >  1:3220016/1‑93 (MQ=255)

                          aCGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCaa                                                               <  2:7795854/93‑1 (MQ=255)

                            gCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGc                                                             >  1:6644704/1‑93 (MQ=255)

                             cTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCt                                                            >  1:7232928/1‑93 (MQ=255)

                             cTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCt                                                            <  1:7375938/93‑1 (MQ=255)

                              ttCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTg                                                           <  2:9677604/93‑1 (MQ=255)

                              ttCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTg                                                           >  1:4632203/1‑93 (MQ=255)

                               tCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGa                                                          <  1:9841537/93‑1 (MQ=255)

                                ccAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGAc                                                         >  1:4137620/1‑93 (MQ=255)

                                 cAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                        <  2:6488627/93‑1 (MQ=255)

                                 cAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                        >  1:4043639/1‑93 (MQ=255)

                                  aGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACgg                                                       >  1:2982845/1‑92 (MQ=255)

                                  aGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGt                                                       <  2:4278601/93‑1 (MQ=255)

                                  aGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGt                                                       <  2:4823414/93‑1 (MQ=255)

                                   ggAACAGATGGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTg                                                      <  1:1941056/93‑1 (MQ=255)

                                   ggAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTg                                                      <  2:4366826/93‑1 (MQ=255)

                                   ggAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTg                                                      >  2:7957865/1‑93 (MQ=255)

                                    gAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTgg                                                     >  1:5295651/1‑93 (MQ=255)

                                      aCAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGAt                                                   >  2:1895874/1‑93 (MQ=255)

                                       cAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATa                                                  <  2:6037623/93‑1 (MQ=255)

                                         gATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATaaa                                                <  2:6940931/93‑1 (MQ=255)

                                         gATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATaaa                                                <  2:808675/93‑1 (MQ=255)

                                         gATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATaaa                                                <  2:8775413/93‑1 (MQ=255)

                                         gATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATaaa                                                <  2:5194756/93‑1 (MQ=255)

                                          aTAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCTGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGAGGTGGATAAAc                                               <  1:8068208/93‑1 (MQ=255)

                                          aTAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTTTGACTCATCAGGTCACAGCCCAATGGCGCAAGATGACGTGGATAAAc                                               >  2:4137392/1‑93 (MQ=255)

                                           tAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACg                                              >  2:5255468/1‑93 (MQ=255)

                                           tAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACg                                              >  1:5464491/1‑93 (MQ=255)

                                           tAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACg                                              >  1:2718092/1‑93 (MQ=255)

                                            aGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGc                                             >  2:8707203/1‑93 (MQ=255)

                                            aGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGc                                             <  2:7988436/93‑1 (MQ=255)

                                             gCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGcc                                            >  1:4255978/1‑93 (MQ=255)

                                             gCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGcc                                            <  1:4403578/93‑1 (MQ=255)

                                              cGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcg                                           <  2:9218695/93‑1 (MQ=255)

                                               ggTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgt                                          >  1:4876800/1‑93 (MQ=255)

                                                gtcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtc                                         <  2:4502004/93‑1 (MQ=255)

                                                gtcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAACCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtc                                         <  2:3131501/93‑1 (MQ=255)

                                                   gtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGtt                                      <  2:1204193/93‑1 (MQ=255)

                                                   gtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGtt                                      >  2:2053379/1‑93 (MQ=255)

                                                    tcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTc                                     >  1:4854759/1‑93 (MQ=255)

                                                    tcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTc                                     <  1:8215979/93‑1 (MQ=255)

                                                     cATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCa                                    <  1:8800585/93‑1 (MQ=255)

                                                     cATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCa                                    >  2:6521859/1‑93 (MQ=255)

                                                       tCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCCtt                                  >  2:2316422/1‑93 (MQ=255)

                                                       tCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCAtt                                  >  1:7701257/1‑93 (MQ=255)

                                                       tCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCAtt                                  <  1:8561938/93‑1 (MQ=255)

                                                       tCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCAtt                                  >  2:9019525/1‑93 (MQ=255)

                                                        cgcgACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTc                                 >  2:3873430/1‑93 (MQ=255)

                                                         gcgACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTcc                                <  1:9428706/93‑1 (MQ=255)

                                                           gACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCca                              <  1:7099926/93‑1 (MQ=255)

                                                              ggCTACGGTCGCTACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagca                           <  1:9387364/93‑1 (MQ=255)

                                                              ggCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagca                           <  1:7981219/93‑1 (MQ=255)

                                                                cTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAgg                         <  2:254049/93‑1 (MQ=255)

                                                                 tACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGc                        <  1:7544887/93‑1 (MQ=255)

                                                                 tACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGc                        <  2:8446716/93‑1 (MQ=255)

                                                                   cGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCaa                      <  2:5921704/93‑1 (MQ=255)

                                                                      tCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAAcc                   >  2:3027071/1‑93 (MQ=255)

                                                                       cGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCa                  <  1:7305381/93‑1 (MQ=255)

                                                                        gCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCAc                 <  1:3562746/93‑1 (MQ=255)

                                                                         ccACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACt                >  2:8293026/1‑93 (MQ=255)

                                                                           aCGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACttt              >  2:3406623/1‑93 (MQ=255)

                                                                            cGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTg             <  2:1303652/93‑1 (MQ=255)

                                                                            cGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTg             <  2:7318623/93‑1 (MQ=255)

                                                                            cGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTg             >  1:2063287/1‑93 (MQ=255)

                                                                            cGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTg             >  2:4224670/1‑93 (MQ=255)

                                                                              cTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGaa           <  2:5693947/93‑1 (MQ=255)

                                                                              cTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGaa           <  1:2581607/93‑1 (MQ=255)

                                                                               ttCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAAc          <  1:1165495/93‑1 (MQ=255)

                                                                               ttCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAAc          <  1:5633460/93‑1 (MQ=255)

                                                                               ttCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAAc          >  1:1013617/1‑93 (MQ=255)

                                                                                tCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACg         >  2:1734148/1‑93 (MQ=255)

                                                                                tCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACg         >  2:2217251/1‑93 (MQ=255)

                                                                                  ggTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGa       <  1:3552542/93‑1 (MQ=255)

                                                                                   gTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGaa      <  2:7484841/93‑1 (MQ=255)

                                                                                   gTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGaa      <  2:3611814/93‑1 (MQ=255)

                                                                                      tCTGACTCGTCACGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTg   <  1:5134700/93‑1 (MQ=255)

ATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGA  >  NC_000913/2958267‑2958446

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