BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,623,776	G→A	E245K (GAG→AAG)	ppk →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,623,776	0	G	A	99.0%	741.6 / ‑3.5	206	E245K (GAG→AAG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (92/112); total (93/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 209 total aligned reads displayed.

AGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTT  >  NC_000913/2623679‑2623870

aGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAg                                                                                              <  1:236425/100‑1 (MQ=255)

      cttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATg                                                                                        <  1:4838652/100‑1 (MQ=255)

      cttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATg                                                                                        <  1:771859/100‑1 (MQ=255)

        ttGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTc                                                                                      >  1:4404672/1‑100 (MQ=255)

          gATTATGACGCGCTGAATGCCTATTCAGTGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCtt                                                                                    >  1:5499369/1‑100 (MQ=255)

           aTTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTc                                                                                   >  1:5694373/1‑100 (MQ=255)

            ttATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTcc                                                                                  <  1:2984151/100‑1 (MQ=255)

             tATGACGCGCTCAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGtt                                                                                            <  2:2472468/89‑1 (MQ=255)

              aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTg                                                                                           <  2:3787479/89‑1 (MQ=255)

              aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTg                                                                                           <  2:5352309/89‑1 (MQ=255)

              aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAACCCAGCCTGATGAAGTTGATGTCTTCCAg                                                                                <  1:5033693/100‑1 (MQ=255)

               tGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGa                                                                                          >  2:5553947/1‑89 (MQ=255)

                gACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGAt                                                                                         >  2:3942887/1‑89 (MQ=255)

                 aCGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATg                                                                                        <  2:485928/89‑1 (MQ=255)

                 aCGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATg                                                                                        <  2:4718090/89‑1 (MQ=255)

                 aCGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATg                                                                                        <  2:5115630/89‑1 (MQ=255)

                  cgcgCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGtctc                                                                            <  1:2320777/100‑1 (MQ=255)

                  cgcgCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGtctc                                                                            <  1:1244064/100‑1 (MQ=255)

                       tGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAagca                                                                       >  1:284409/1‑100 (MQ=255)

                        gAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAagcag                                                                      >  1:5084446/1‑100 (MQ=255)

                           tGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGt                                                                   <  1:2447032/100‑1 (MQ=255)

                           tGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGt                                                                   >  1:1890630/1‑100 (MQ=255)

                            gCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGtct                                                                             <  2:416687/89‑1 (MQ=255)

                             ccTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGttt                                                                 <  1:117165/100‑1 (MQ=255)

                              cTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCa                                                                           <  2:5449282/89‑1 (MQ=255)

                              cTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTa                                                                <  1:3823497/100‑1 (MQ=255)

                               tATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTaa                                                               >  1:362419/1‑100 (MQ=255)

                                aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGGCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAag                                                                         <  2:2909804/89‑1 (MQ=255)

                                aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAAc                                                              >  1:1314189/1‑100 (MQ=255)

                                 ttCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAAct                                                             <  1:1610848/100‑1 (MQ=255)

                                   cAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAagcag                                                                      >  2:1061858/1‑89 (MQ=255)

                                   cAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAagcag                                                                      >  2:1575467/1‑89 (MQ=255)

                                   cAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgc                                                           <  1:1233444/100‑1 (MQ=255)

                                   cAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgc                                                           >  1:5590222/1‑100 (MQ=255)

                                   cAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgc                                                           >  1:4540678/1‑100 (MQ=255)

                                    aaTGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgct                                                          >  1:935927/1‑100 (MQ=255)

                                     aTGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCg                                                                    >  2:58329/1‑89 (MQ=255)

                                     aTGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgctg                                                         >  1:2214878/1‑100 (MQ=255)

                                      tGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGa                                                        >  1:2052577/1‑100 (MQ=255)

                                       gAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGtt                                                                  >  2:5043921/1‑89 (MQ=255)

                                       gAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAg                                                       >  1:395305/1‑100 (MQ=255)

                                       gAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAg                                                       <  1:3445594/100‑1 (MQ=255)

                                        aaGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGttt                                                                 <  2:3053570/89‑1 (MQ=255)

                                         aGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGcc                                                     <  1:47718/100‑1 (MQ=255)

                                          gATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTaa                                                               >  2:2928283/1‑89 (MQ=255)

                                          gATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCg                                                    <  1:2662258/100‑1 (MQ=255)

                                           aTGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAAc                                                              <  2:1651102/89‑1 (MQ=255)

                                           aTGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCgg                                                   <  1:3865795/100‑1 (MQ=255)

                                            tGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAAct                                                             >  2:5200389/1‑89 (MQ=255)

                                              aCCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgc                                                           <  2:1261860/89‑1 (MQ=255)

                                               cccGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                                               >  1:108552/1‑100 (MQ=255)

                                                      tGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGtt                                        >  1:3315622/1‑100 (MQ=255)

                                                        ccGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTTACTGCTGAGCCGGTGCGTTTTGTTTa                                      <  1:190841/100‑1 (MQ=255)

                                                        ccGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTg                                                 >  2:342974/1‑89 (MQ=255)

                                                         cGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGc                                                <  2:4933802/89‑1 (MQ=255)

                                                            aTACGATTTAGTGCATGAGGTGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAg                                  <  1:3271803/100‑1 (MQ=255)

                                                             tACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgc                                 >  1:2748312/1‑100 (MQ=255)

                                                             tACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgc                                 <  1:3220281/100‑1 (MQ=255)

                                                               cGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTg                                          >  2:1131448/1‑89 (MQ=255)

                                                               cGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgcgc                               <  1:133972/100‑1 (MQ=255)

                                                               cGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgcgc                               >  1:4579965/1‑100 (MQ=255)

                                                                gATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGt                                         <  2:2813046/89‑1 (MQ=255)

                                                                gATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGcgcg                              <  1:5553947/100‑1 (MQ=255)

                                                                 aTTTAGTGCATGAGCTGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGtt                                        >  2:1580146/1‑89 (MQ=255)

                                                                 aTTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGtt                                        <  2:666474/89‑1 (MQ=255)

                                                                  tttAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGttt                                       <  2:4639083/89‑1 (MQ=255)

                                                                  tttAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGttt                                       <  2:4451127/89‑1 (MQ=255)

                                                                  tttAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGttt                                       <  2:4085738/89‑1 (MQ=255)

                                                                  tttAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGttt                                       <  2:62935/89‑1 (MQ=255)

                                                                   ttAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTa                                      <  2:2620276/89‑1 (MQ=255)

                                                                   ttAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTa                                      <  2:3472111/89‑1 (MQ=255)

                                                                       tGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAg                                  <  2:5370295/89‑1 (MQ=255)

                                                                       tGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGcc                       <  1:4066724/100‑1 (MQ=255)

                                                                        gCATGAGATGGAAGCCCGCCTGATGAAGTGGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGCTCTGCTTCTCAgc                                 >  2:1943801/1‑89 (MQ=255)

                                                                        gCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGccc                      >  1:3124237/1‑100 (MQ=255)

                                                                        gCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGccc                      <  1:3531509/100‑1 (MQ=255)

                                                                        gCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGccc                      >  1:324167/1‑100 (MQ=255)

                                                                         cATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCa                     <  1:4643351/100‑1 (MQ=255)

                                                                         cATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCa                     >  1:4181279/1‑100 (MQ=255)

                                                                          aTGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgcgc                               >  2:3095355/1‑89 (MQ=255)

                                                                          aTGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCaa                    <  1:5155792/100‑1 (MQ=255)

                                                                           tGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGcgcg                              <  2:4516528/89‑1 (MQ=255)

                                                                             agaTGGGAGGCAGCCGGTTGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTAAGCCGGTGCGATTTGTTTATGAGCGCGag                            >  2:2786992/1‑88 (MQ=255)

                                                                             agaTGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATgc                 >  1:2512325/1‑100 (MQ=255)

                                                                              gaTGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTTCGTTTTGTTTATCAGCGCGata                           >  2:715762/1‑89 (MQ=255)

                                                                                tGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCt              <  1:370605/100‑1 (MQ=255)

                                                                                 ggAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGc                        >  2:523846/1‑89 (MQ=255)

                                                                                 ggAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGc                        >  2:2136769/1‑89 (MQ=255)

                                                                                 ggAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTg             >  1:5283498/1‑100 (MQ=255)

                                                                                  gAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGcc                       <  2:4314413/89‑1 (MQ=255)

                                                                                  gAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGcc                       >  2:3276265/1‑89 (MQ=255)

                                                                                  gAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTgg            >  1:1554633/1‑100 (MQ=255)

                                                                                  gAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTgg            >  1:180743/1‑100 (MQ=255)

                                                                                     gccagccTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTg         >  1:924666/1‑100 (MQ=255)

                                                                                       cagccTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATg                  <  2:1953880/89‑1 (MQ=255)

                                                                                        agccTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAg      >  1:3188410/1‑100 (MQ=255)

                                                                                        agccTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAg      <  1:4453968/100‑1 (MQ=255)

                                                                                         gccTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAgt     >  1:702266/1‑100 (MQ=255)

                                                                                            tgatgaAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTg             >  2:227941/1‑89 (MQ=255)

                                                                                            tgatgaAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAAAGCGCTg             >  2:2732741/1‑89 (MQ=255)

AGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTT  >  NC_000913/2623679‑2623870

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