BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,926,882	C→A	19.1%	G653V (GGC→GTC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,926,882	0	C	A	19.1%	406.9 / 76.2	166	G653V (GGC→GTC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base C (71/63); new base A (14/18); total (85/81)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 169 total aligned reads displayed.

TCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGC  >  NC_000913/1926783‑1926980

     aGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATgg                                                                                               <  1:3485903/100‑1 (MQ=255)

      ggCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGc                                                                                              >  1:3369370/1‑100 (MQ=255)

      ggCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGc                                                                                              <  1:3796625/100‑1 (MQ=255)

      cgatctgCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGc                                                                                              <  1:3204518/94‑1 (MQ=255)

       gCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGcc                                                                                             >  1:731532/1‑100 (MQ=255)

        cGTAGCGGGTAATGTTCCCTGCGCCAGCGCTACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                                                                                            <  1:4290351/100‑1 (MQ=255)

         gTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                                                                                            >  1:1892313/1‑99 (MQ=255)

         gTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTg                                                                                           <  2:4331863/100‑1 (MQ=255)

            gCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                                        <  2:3818118/100‑1 (MQ=255)

            gCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                                        <  1:3310426/100‑1 (MQ=255)

             cGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGTCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                                                       <  1:2443840/100‑1 (MQ=255)

             cGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                                                       <  2:3109106/100‑1 (MQ=255)

              gggTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTcc                                                                                      >  1:2564387/1‑100 (MQ=255)

               ggTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                                                                                     >  1:2475208/1‑100 (MQ=255)

               ggTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                                                                                     <  1:3924796/100‑1 (MQ=255)

               ggTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATggactcaggcca                                                                                     <  2:2661679/100‑11 (MQ=255)

                 gtaTGTTCCTCTCGGCATCGCGACCATAATAGCATATTCCATCGTTTTTCCTTCTTACTATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAt                                                                                    <  2:3679014/97‑1 (MQ=255)

                   aTGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                 >  2:3165026/1‑100 (MQ=255)

                    tGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg                                                                                >  2:414461/1‑100 (MQ=255)

                    tGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg                                                                                <  2:828023/100‑1 (MQ=255)

                     gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                               >  1:98876/1‑100 (MQ=255)

                     gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                               <  2:3614100/100‑1 (MQ=255)

                     gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                               <  1:791442/100‑1 (MQ=255)

                     gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                               >  1:2332432/1‑100 (MQ=255)

                     gTTCCCTGCGCCAGCGCGACCAGAAAAGCATAATCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCTCCATCGCCTGAATCCATGTCg                                                                                >  1:1895784/1‑99 (MQ=255)

                      ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCggn                                                                              <  2:1541762/100‑2 (MQ=38)

                      ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGt                                                                              <  2:1388183/100‑1 (MQ=38)

                       tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                             >  1:1852818/1‑100 (MQ=255)

                        cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaa                                                                            <  2:3279350/100‑2 (MQ=37)

                         ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                           <  2:461241/100‑1 (MQ=255)

                         ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaa                                                                           <  1:1014626/100‑3 (MQ=37)

                         ccTCCGTCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCTTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                             <  1:1710916/98‑1 (MQ=255)

                           tGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatc                                                                         <  2:1125792/100‑5 (MQ=25)

                             cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACagag                                                                        >  1:2518828/1‑99 (MQ=255)

                             cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc                                                                       >  2:3563965/1‑100 (MQ=255)

                              gcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAGGTCGGTACAGAGCa                                                                      >  2:1906250/1‑100 (MQ=255)

                                cAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctggtc                                                                    <  1:1858161/100‑10 (MQ=21)

                                cAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctggt                                                                     <  1:3298216/99‑9 (MQ=21)

                                  gcgcGACCAGAAAAGCATATTCCATCGCTACGCCTTAGTACGAGTTAAAGCAACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaa                                                                  >  1:2382811/1‑100 (MQ=255)

                                    gtgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                <  2:4397843/98‑1 (MQ=255)

                                    gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                <  2:4282756/100‑1 (MQ=255)

                                    gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                <  2:3155551/100‑1 (MQ=255)

                                     cgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACGGCCATGACCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                               >  1:3090915/1‑100 (MQ=255)

                                      gACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGa                                                              >  2:2265545/1‑100 (MQ=255)

                                         cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                           <  2:1620539/100‑1 (MQ=255)

                                          aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                          >  1:2303774/1‑100 (MQ=255)

                                          aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                          <  2:1242734/100‑1 (MQ=255)

                                               aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGtacggc                                                     >  2:1728919/1‑95 (MQ=255)

                                               aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCg                                                     >  2:1189190/1‑100 (MQ=255)

                                                gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCgg                                                    <  1:74069/100‑1 (MQ=255)

                                                 cATATTCCATCGCTACGCCTTCGTACGATTTAAGACGACCAGATTTACCGCCATGGCCTGAGTCCATGCCGGTAAAGAACAATAAAAGCCGCTCATCGGt                                                   >  1:4391296/1‑100 (MQ=255)

                                                   tatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAAAGCAATAAAAGATGGTCATCGGttt                                                 >  2:2894410/1‑100 (MQ=255)

                                                    atTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACCGAGCAATAAAAGATGGTAATCGGtttt                                                >  1:4188980/1‑100 (MQ=255)

                                                     ttacATCGCTACGCGTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                               <  1:3549581/97‑1 (MQ=255)

                                                     ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                               >  2:375489/1‑100 (MQ=255)

                                                      tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCGATAAAAGATGGTCATCGGTTTTc                                               <  1:541652/99‑1 (MQ=255)

                                                      tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                              <  1:1001494/100‑1 (MQ=255)

                                                        cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                            <  1:3162338/100‑1 (MQ=255)

                                                        cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                            >  1:2858849/1‑100 (MQ=255)

                                                        cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                            >  1:2856067/1‑100 (MQ=255)

                                                        cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                            >  1:905662/1‑100 (MQ=255)

                                                        cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCACGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                            <  1:4036702/100‑1 (MQ=255)

                                                         aTCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCt                                           <  1:1943791/100‑1 (MQ=255)

                                                          tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTc                                          <  1:2921580/100‑1 (MQ=255)

                                                             cTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                       >  2:2250451/1‑100 (MQ=255)

                                                                    ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTa                                >  2:4225196/1‑100 (MQ=255)

                                                                     tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAg                               <  2:691261/100‑1 (MQ=255)

                                                                      cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGc                              >  1:4198643/1‑100 (MQ=255)

                                                                        tACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGa                            >  1:3069804/1‑100 (MQ=255)

                                                                         ncGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAACGATGGTCATCGGCTTTCAGCTCGCGCAATTTAGCGAc                           >  1:1226103/2‑100 (MQ=255)

                                                                          cGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAcc                          <  2:3978833/100‑1 (MQ=255)

                                                                           gATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAccc                         >  1:164138/1‑100 (MQ=255)

                                                                            aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCa                        >  1:1418886/1‑100 (MQ=255)

                                                                             tttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt                       >  2:4224635/1‑100 (MQ=255)

                                                                               tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAttt                     >  1:4367691/1‑100 (MQ=255)

                                                                               tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAttt                     >  1:4262556/1‑100 (MQ=255)

                                                                               tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAttt                     >  2:2134325/1‑100 (MQ=255)

                                                                                  aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                  <  2:1475253/100‑1 (MQ=255)

                                                                                  aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                  >  2:4077247/1‑100 (MQ=255)

                                                                                  aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                  <  1:230596/100‑1 (MQ=255)

                                                                                  aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAGCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                  <  1:2362509/100‑1 (MQ=255)

                                                                                    cGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCg                >  1:1309254/1‑100 (MQ=255)

                                                                                     gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg               >  2:3530943/1‑100 (MQ=255)

                                                                                     gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg               <  2:3076171/100‑1 (MQ=255)

                                                                                         aGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc           >  2:3796116/1‑100 (MQ=255)

                                                                                          gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc          <  2:2015586/100‑1 (MQ=255)

                                                                                          gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc          <  2:3061583/100‑1 (MQ=255)

                                                                                          gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc           >  1:1349070/1‑99 (MQ=255)

                                                                                            tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCaa        >  2:1054586/1‑100 (MQ=255)

                                                                                            tgtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCAGGTTcac          >  2:3018956/3‑96 (MQ=25)

                                                                                               aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat     >  2:1670430/1‑100 (MQ=38)

                                                                                               aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat     >  2:1695112/1‑100 (MQ=38)

                                                                                               aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat     >  2:545684/1‑100 (MQ=38)

                                                                                               aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat     >  1:2173298/1‑100 (MQ=38)

                                                                                               aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGATTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat     >  2:3385231/1‑100 (MQ=38)

                                                                                                ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt    <  1:1650476/100‑1 (MQ=255)

                                                                                                ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt    >  1:3968915/1‑100 (MQ=38)

                                                                                                ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt    >  2:4345037/1‑100 (MQ=38)

                                                                                                 cGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg   <  1:1853513/100‑1 (MQ=255)

                                                                                                  gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc  <  2:3449574/100‑1 (MQ=255)

TCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGC  >  NC_000913/1926783‑1926980

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