BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,750,373	A→G	16.4%	A112A (GCT→GCC)	grpE ←	heat shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,750,373	0	A	G	16.4%	369.8 / 62.9	165	A112A (GCT→GCC)	grpE	heat shock protein
Reads supporting (aligned to +/- strand): ref base A (69/68); new base G (10/17); total (80/85)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 172 total aligned reads displayed.

AGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCAATA  >  NC_000913/2750274‑2750472

aGTTTCGGCGATCACTTCAACGCTAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCATACATATCCGGGTTa                                                                                                     <  2:2131671/100‑1 (MQ=255)

aGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTa                                                                                                     <  1:2654394/100‑1 (MQ=255)

 gTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAg                                                                                                    <  2:1938000/100‑1 (MQ=255)

       gCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTAt                                                                                              <  2:1083955/100‑1 (MQ=255)

        cGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTGGCTTTATc                                                                                             >  2:2606167/1‑100 (MQ=255)

        cGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATc                                                                                             >  1:3210392/1‑100 (MQ=255)

        cGATCACTTCAACGCCAAACATACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATc                                                                                             >  2:2337974/1‑100 (MQ=255)

         gATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCa                                                                                            >  2:1785867/1‑100 (MQ=255)

          aTCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAg                                                                                           <  2:981384/100‑1 (MQ=255)

            cACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGcc                                                                                         <  1:2764417/100‑1 (MQ=255)

            cACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGcc                                                                                         <  1:3096805/100‑1 (MQ=255)

             aCTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCa                                                                                        <  1:2337989/100‑1 (MQ=255)

             aCTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCa                                                                                        >  2:453634/1‑100 (MQ=255)

             aCTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCa                                                                                        <  1:315053/100‑1 (MQ=255)

             aCTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCa                                                                                        >  1:840805/1‑100 (MQ=255)

              catCAACGCCAAACCTACGCACCCCAACCACCCTCGACTTCAGCGTCAGCTCAATGCCTTCAACCCCCCCCCACATATCCGCATCCGCCTTATCACCCAc                                                                                       >  2:747153/3‑100 (MQ=255)

              cTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCAc                                                                                       >  1:865345/1‑100 (MQ=255)

              cTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCAc                                                                                       >  2:1012015/1‑100 (MQ=255)

              cTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCAc                                                                                       >  1:3291448/1‑100 (MQ=255)

                 cAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACt                                                                                      >  2:3404080/1‑98 (MQ=255)

                 cAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTc                                                                                    >  2:120603/1‑100 (MQ=255)

                     gCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAgc                                                                                <  1:2054378/100‑1 (MQ=255)

                      ccAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAgcg                                                                               >  2:2986884/1‑100 (MQ=255)

                      ccAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAgcg                                                                               <  1:3647746/100‑1 (MQ=255)

                      caaaaTTTACGGACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATACGGGTTGGCTTTATCAGCCACTTcc                                                                                   <  2:1681975/94‑1 (MQ=255)

                       cAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAgcgc                                                                              >  2:3715266/1‑100 (MQ=255)

                        aaaCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCa                                                                             <  1:921178/100‑1 (MQ=255)

                        aaaCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCa                                                                             >  1:1975376/1‑100 (MQ=255)

                            ttACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGAt                                                                         <  2:251414/100‑1 (MQ=255)

                            ttACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGAt                                                                         <  1:60541/100‑1 (MQ=255)

                               cGCACAACGTCGAGCATCGGGTTCAGCGCCGGCTCAATGCCTTCAACCAGCGCAGAGATACCCGGGGCAGCCTTGCTACCCACGGCCGGCGCACGCTcc                                                                       >  2:2186617/1‑99 (MQ=255)

                               cGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCa                                                                      <  1:722049/100‑1 (MQ=255)

                                 cacaACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAgg                                                                    >  1:1796747/1‑100 (MQ=255)

                                 cacaACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAgg                                                                    >  2:2095583/1‑100 (MQ=255)

                                 cacaACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAgg                                                                    >  1:2180677/1‑100 (MQ=255)

                                 cacaACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAgg                                                                    <  1:2831406/100‑1 (MQ=255)

                                  acaacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGc                                                                   <  2:1386881/100‑1 (MQ=255)

                                  acaacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGc                                                                   >  1:3354161/1‑100 (MQ=255)

                                   caacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCt                                                                  >  1:947309/1‑100 (MQ=255)

                                   caacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCt                                                                  >  2:3644300/1‑100 (MQ=255)

                                        tCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaa                                                             <  2:149268/100‑1 (MQ=255)

                                         ccAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATCTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaat                                                            <  2:1981304/100‑1 (MQ=255)

                                            gCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGACCCAGGCTATCAATCAc                                                         <  1:1065038/100‑1 (MQ=255)

                                             cATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCAcc                                                        >  1:3344305/1‑100 (MQ=255)

                                             cATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCAcc                                                        <  1:1132098/100‑1 (MQ=255)

                                              aTCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTCCCAGCGCACGATCCAGGCTATCAATCACCg                                                       >  2:1130390/1‑100 (MQ=255)

                                               tCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCgg                                                      >  2:569689/1‑100 (MQ=255)

                                                cGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgc                                                     <  1:1533859/100‑1 (MQ=255)

                                                cGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgc                                                     <  2:2800076/100‑1 (MQ=255)

                                                cGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgc                                                     <  1:2104318/100‑1 (MQ=255)

                                                 gACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgca                                                    <  2:2381871/100‑1 (MQ=255)

                                                 gACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgca                                                    >  2:243038/1‑100 (MQ=255)

                                                 gACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgca                                                    >  1:1606754/1‑100 (MQ=255)

                                                 gACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgca                                                    <  2:3577540/100‑1 (MQ=255)

                                                 gACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGCGTTAGCTTTATCAGCCACTTACAGCGCACGATCCAGGCTATCAATCACCGgc                                                     >  1:1305591/1‑99 (MQ=255)

                                                  aCTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgcag                                                   <  1:2522807/100‑1 (MQ=255)

                                                   cTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgcagc                                                  >  2:96021/1‑100 (MQ=255)

                                                       aGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAAtt                                              >  1:2829592/1‑100 (MQ=255)

                                                        gCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTc                                             <  2:2881705/100‑1 (MQ=255)

                                                           tCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGtt                                          <  2:20330/100‑1 (MQ=255)

                                                           tCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGtt                                          <  1:1435031/100‑1 (MQ=255)

                                                            cAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtg                                         <  2:1225730/100‑1 (MQ=255)

                                                            cAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtg                                         <  1:2067633/100‑1 (MQ=255)

                                                             aGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtga                                        <  1:2885643/100‑1 (MQ=255)

                                                              gCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtgat                                       >  2:2011242/1‑100 (MQ=255)

                                                               cTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtgatg                                      <  2:1296563/100‑1 (MQ=255)

                                                                tCAATGCCTTCAACCATCGCAGACATATCCGGGTGGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtg                                         <  2:2068622/96‑1 (MQ=255)

                                                                 cAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGaa                                    >  1:1719940/1‑100 (MQ=255)

                                                                  aaTGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAAt                                   >  1:2248244/1‑100 (MQ=255)

                                                                    tGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAAttt                                 <  1:1137838/100‑1 (MQ=255)

                                                                      ccTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTtct                               <  2:1118787/100‑1 (MQ=255)

                                                                      ccTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTtct                               >  1:1008943/1‑100 (MQ=255)

                                                                       cTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTtctc                              >  1:2355009/1‑100 (MQ=255)

                                                                       cTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTtctc                              >  1:3150038/1‑100 (MQ=255)

                                                                       cTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTtctc                              >  1:3748616/1‑100 (MQ=255)

                                                                          cAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAg                           <  2:1137596/100‑1 (MQ=255)

                                                                             ccATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAgcgc                        <  2:191083/100‑1 (MQ=255)

                                                                                 cGCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGaa                     >  2:613274/1‑99 (MQ=255)

                                                                                  gCAGACATATCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAAtt                   >  2:834374/1‑100 (MQ=255)

                                                                                    aGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTtg                 >  2:3477571/1‑100 (MQ=255)

                                                                                     gACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTtgt                <  2:3023462/100‑1 (MQ=255)

                                                                                      aCATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTtgtg               <  1:1634878/100‑1 (MQ=255)

                                                                                      aCATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTtgtg               <  2:3143350/100‑1 (MQ=255)

                                                                                       cATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTACAGCGCGAATTTGTgg              >  2:2232626/1‑100 (MQ=255)

                                                                                        atatCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTggg             <  1:560107/100‑1 (MQ=255)

                                                                                         tatCCGGGTTGGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGc            >  2:1763484/1‑100 (MQ=255)

                                                                                         tatCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCATTTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGc            <  1:2690937/100‑1 (MQ=255)

                                                                                          atCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCt           >  2:2697595/1‑100 (MQ=255)

                                                                                          atCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCt           >  1:1193537/1‑100 (MQ=255)

                                                                                            ccGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCCGGCTATCAATCACCGGCCGCAATTCGTTGCAGAATTTCCCCAGCGCGACTTTGTGGCCt           >  2:3205041/1‑98 (MQ=255)

                                                                                               ggTTAGCTTTATCAGCCACTTCCAGCGCACTTTCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCttttgc      <  2:149529/100‑3 (MQ=255)

                                                                                               ggTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTc      <  2:1731237/100‑1 (MQ=255)

                                                                                               ggTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTc      >  2:1618787/1‑100 (MQ=255)

                                                                                                gTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGTTTTTTCa     <  1:2001405/100‑1 (MQ=255)

                                                                                                gTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCa     <  1:2384977/100‑1 (MQ=255)

                                                                                                gTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCa     >  2:440641/1‑100 (MQ=255)

                                                                                                gTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCa     >  1:71184/1‑100 (MQ=255)

                                                                                                  tAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCAat   <  1:1012024/100‑1 (MQ=255)

                                                                                                   aGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCAata  >  1:1489778/1‑100 (MQ=255)

                                                                                                   aGCTTTAGCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCaaca  <  2:215740/100‑3 (MQ=255)

AGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAATTTGTGGGCTTTTTCAATA  >  NC_000913/2750274‑2750472

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