BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,026,011	T→G	H82Q (CAT→CAG)	guaD →	guanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,011	0	T	G	97.7%	578.3 / 3.1	172	H82Q (CAT→CAG)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand): ref base T (1/3); new base G (67/101); total (68/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

Only 100 of 173 total aligned reads displayed.

TTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGA  >  NC_000913/3025914‑3026106

ttGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACagat                                                                                               >  2:1399795/1‑100 (MQ=255)

 tGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATc                                                                                              >  2:730320/1‑100 (MQ=255)

 tGGCGAATGGGAAAACGCAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATc                                                                                              <  2:1081572/100‑1 (MQ=255)

  ggCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATcc                                                                                             <  2:2699054/100‑1 (MQ=255)

  ggCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATcc                                                                                             >  2:2909702/1‑100 (MQ=255)

   gCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCa                                                                                            <  1:1056471/100‑1 (MQ=255)

    cGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCAt                                                                                           >  1:3835799/1‑100 (MQ=255)

    cGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCAt                                                                                           <  1:1956691/100‑1 (MQ=255)

      aaTGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTa                                                                                         >  2:2026709/1‑100 (MQ=255)

         gggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATcc                                                                                      <  1:3153392/100‑1 (MQ=255)

         gggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATcc                                                                                      <  1:705029/100‑1 (MQ=255)

         gggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATcc                                                                                      <  1:1760095/100‑1 (MQ=255)

          ggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCg                                                                                     >  1:851179/1‑100 (MQ=255)

          ggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCg                                                                                     <  2:2722530/100‑1 (MQ=255)

          ggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCg                                                                                     >  2:28004/1‑100 (MQ=255)

          ggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCg                                                                                     >  2:2151053/1‑100 (MQ=255)

          ggAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCg                                                                                     >  1:1372760/1‑100 (MQ=255)

          ggAAAACGGAAAGCATCAAATTCCTGACACCAGTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCg                                                                                     <  2:3138082/100‑1 (MQ=255)

            aaaaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCa                                                                                   >  2:3473064/1‑100 (MQ=255)

            aaaaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCa                                                                                   <  1:457366/100‑1 (MQ=255)

             aaaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCaa                                                                                  >  1:65746/1‑100 (MQ=255)

               aCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAg                                                                                <  1:3063436/100‑1 (MQ=255)

                 ggAAAGCATCAAATTCCTTACACCATTCCCGTGCGCGACTATCGCCGCAAACTGATAGTACCGGGCTTTGTCGACACACAGATCCATTATCCGCAAAGTg                                                                              <  1:975806/100‑1 (MQ=255)

                    aaGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGaaa                                                                           <  1:3927810/100‑1 (MQ=255)

                     aGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAt                                                                          <  2:1534505/100‑1 (MQ=255)

                     aGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAt                                                                          <  2:2445443/100‑1 (MQ=255)

                      gCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtg                                                                         >  2:2142418/1‑100 (MQ=255)

                      gCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtg                                                                         >  1:3848752/1‑100 (MQ=255)

                       cATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtgg                                                                        >  1:409487/1‑100 (MQ=255)

                        aTCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtggt                                                                       <  1:1749691/100‑1 (MQ=255)

                         tCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtggtg                                                                      <  1:170485/100‑1 (MQ=255)

                         tCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtggtg                                                                      >  1:2566452/1‑100 (MQ=255)

                         tCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTAGATACACAGACCCATTATCCGCAAAGTGAAAtggtg                                                                      >  2:3779968/1‑100 (MQ=255)

                          cAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtggtgg                                                                     >  2:1018738/1‑100 (MQ=255)

                          cAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtggtgg                                                                     <  2:464382/100‑1 (MQ=255)

                            aaTTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAAATGATAGTACCGGGCATTGTCGATACACAGATCCATTATCCGCCAAGTGAAATGGTgggg                                                                   >  2:3307371/1‑100 (MQ=255)

                              ttCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAAtggtgtcagg                                                                 >  1:3958762/1‑95 (MQ=255)

                              ttCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGc                                                                 <  2:3222029/100‑1 (MQ=255)

                               tCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGcc                                                                <  1:2696380/100‑1 (MQ=255)

                               tCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGcc                                                                <  2:2425895/100‑1 (MQ=255)

                               gccTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGAGCCATTATCCGCAAAGTGAAATGGTGGGGGcc                                                                <  2:3946280/99‑1 (MQ=255)

                                 ctaacacCATTCGCTTGAGCGACTATCGCGGCAAACTGTTAGTACCGGGCTTTGTCGATACACAGATCCATTATACGCGAAGTGAAATGGTGGGGGCCTa                                                              <  2:3729431/97‑1 (MQ=255)

                                  tGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTAt                                                             <  1:1928630/100‑1 (MQ=255)

                                   gACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCCTTATCCGCAAAGTGAAATGGAGGGGGCCtgtg                                                            >  1:4078541/1‑100 (MQ=255)

                                    acacCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATgg                                                           <  1:471716/100‑1 (MQ=255)

                                    acacCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATgg                                                           <  2:3774234/100‑1 (MQ=255)

                                        cATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAg                                                       <  2:777102/100‑1 (MQ=255)

                                         aTTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGc                                                      <  1:3228457/100‑1 (MQ=255)

                                          ttCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCa                                                     <  1:1645434/100‑1 (MQ=255)

                                              cgTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTg                                                 <  1:2631471/100‑1 (MQ=255)

                                               gTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGc                                                <  2:546233/100‑1 (MQ=255)

                                               gTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGc                                                <  2:3420469/100‑1 (MQ=255)

                                                 gcgcGACTATCGCGGCAAACTGATAGTATCGGGCTTTGTCGATACACAGATCCAGTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTg                                              <  2:2609142/100‑1 (MQ=255)

                                                 gcgcGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCTCAAAGTGAAATGGTGGGGGCCTATGGTGATCAATTTCTg                                              >  2:2196076/1‑100 (MQ=255)

                                                 gcgcGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTg                                              <  1:881215/100‑1 (MQ=255)

                                                 gcgcGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTg                                              <  1:3494905/100‑1 (MQ=255)

                                                 gcgcGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTg                                              <  1:3000029/100‑1 (MQ=255)

                                                  cgcgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTgg                                             <  2:1611202/100‑1 (MQ=255)

                                                    cgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAg                                           <  2:2632462/100‑1 (MQ=255)

                                                    cgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAg                                           <  2:921142/100‑1 (MQ=255)

                                                    cgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAg                                           <  2:188567/100‑1 (MQ=255)

                                                    cgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAg                                           >  2:1734096/1‑100 (MQ=255)

                                                    cgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAg                                           <  1:2323551/100‑1 (MQ=255)

                                                    cgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAg                                           <  1:2484613/100‑1 (MQ=255)

                                                        tATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGt                                       <  1:1525643/100‑1 (MQ=255)

                                                         aTCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGtt                                      <  1:248118/100‑1 (MQ=255)

                                                          tCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTg                                     <  2:4002145/100‑1 (MQ=255)

                                                          tCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTg                                     <  2:2476674/100‑1 (MQ=255)

                                                           cgcgGCAAACTGATAGTACTGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGa                                    <  2:3577250/100‑1 (MQ=255)

                                                            gcgGCAAACTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGaa                                   <  1:3961319/100‑1 (MQ=255)

                                                                cAAACTGATAGTACCCGGCTTTGTCGATACACAGACCCATTATCCGCAACGCGCAATGCCGGGGGCCTCTGGCGCGCAATTGCTCGAGCCGTTGACTaaa                               >  2:1155035/1‑100 (MQ=255)

                                                                 aaaCTGCTAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAac                              <  2:3801169/100‑1 (MQ=255)

                                                                  aaCTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAaca                             <  1:2100934/100‑1 (MQ=255)

                                                                  aaCTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAaca                             <  2:2251066/100‑1 (MQ=255)

                                                                  aaCTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAaca                             >  1:2708541/1‑100 (MQ=255)

                                                                   aCTGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacac                            <  1:147246/100‑1 (MQ=255)

                                                                     tGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAACGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGGGGTTGAATAAacacac                          >  1:2100395/1‑100 (MQ=255)

                                                                     tGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacacac                          <  2:2372149/100‑1 (MQ=255)

                                                                     tGATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacacac                          >  2:1241804/1‑100 (MQ=255)

                                                                      gATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACAcc                         <  2:3857261/100‑1 (MQ=255)

                                                                      gATAGTACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACAcc                         <  1:937118/100‑1 (MQ=255)

                                                                      gATAGTACCGGGCTTTGTCGATACACAGATCCATAATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACAcc                         >  2:3696764/1‑100 (MQ=255)

                                                                           tACCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTccc                    <  1:3776699/100‑1 (MQ=255)

                                                                            aCCGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCt                   <  2:641888/100‑1 (MQ=255)

                                                                             ccGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTa                  <  1:1716714/100‑1 (MQ=255)

                                                                              cGGGCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTAc                 >  2:3848332/1‑100 (MQ=255)

                                                                               gggCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACt                <  2:1805762/100‑1 (MQ=255)

                                                                               gggCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACt                <  2:3924837/100‑1 (MQ=255)

                                                                                ggCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTg               <  2:2809892/100‑1 (MQ=255)

                                                                                 gCTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGa              <  2:1924306/100‑1 (MQ=255)

                                                                                  cTTTGTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGaa             >  1:312860/1‑100 (MQ=255)

                                                                                      gTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtc         >  2:773925/1‑100 (MQ=255)

                                                                                      gTCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtc         <  1:3297584/100‑1 (MQ=255)

                                                                                       tCGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtcg        <  1:28446/100‑1 (MQ=255)

                                                                                        cGATACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtcgt       <  2:1452956/100‑1 (MQ=255)

                                                                                          aTACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTa     <  2:1657719/100‑1 (MQ=255)

                                                                                           tACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTAt    <  1:909602/100‑1 (MQ=255)

                                                                                           tACACAGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTAt    <  1:542673/100‑1 (MQ=255)

                                                                                             cacaGATCCATTATCCGCAAAGTGAAATGGTGTGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGa  >  1:3882842/1‑100 (MQ=255)

                                                                                             cacaGATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGa  >  2:1346238/1‑100 (MQ=255)

TTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGA  >  NC_000913/3025914‑3026106

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