Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,600,653 C→T 7.6% S177S (TCG→TCA bglX ← beta‑D‑glucoside glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,600,6530CT 7.6% 151.4 / 2.3 53S177S (TCG→TCAbglXbeta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base T (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

TTTTCGGCAGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTTGGCGACGGCCT  >  NC_002947/1600505‑1600801
                                                                                                                                                     |                                                                                                                                                    
ttttCGGCAGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATtctc                                                                                                                                                      <  4:83740/149‑1 (MQ=255)
ttttCGGCAGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATtctc                                                                                                                                                      <  5:198234/149‑1 (MQ=255)
    cGGCAGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGaa                                                                                                                                                   <  4:130480/148‑1 (MQ=255)
     ggCAGGCTCATATCGACCGTGGTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTCCTGCCCAGGACAGAGCGCCCCCTTACCTGGCCGATTCTCGAAg                                                                                                                                                  >  2:276258/1‑148 (MQ=255)
      gCAGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTc                                                                                                                                                <  5:119226/149‑1 (MQ=255)
      gCAGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTc                                                                                                                                                <  6:232879/149‑1 (MQ=255)
        aGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCag                                                                                                                                              <  3:231396/149‑1 (MQ=255)
         ggCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCagg                                                                                                                                             <  2:210745/149‑1 (MQ=255)
          gCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGaa                                                                                                                                                   <  5:196249/142‑1 (MQ=255)
          gCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGaa                                                                                                                                                   >  6:196249/1‑142 (MQ=255)
          gCTCATATCGACCGTGTTTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGCACGAGCGCACCATTACCTGGCGGATTCTCGAAGTCagg                                                                                                                                             >  1:287961/1‑147 (MQ=255)
           cTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCaggta                                                                                                                                           <  7:196480/149‑1 (MQ=255)
            tCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCaggtag                                                                                                                                          >  2:19090/1‑149 (MQ=255)
            tCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCaggt                                                                                                                                            <  3:286138/147‑1 (MQ=255)
            tCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCaggt                                                                                                                                            >  4:286138/1‑147 (MQ=255)
              atatCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCaggtaggt                                                                                                                                        <  5:277720/149‑1 (MQ=255)
               tatCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGgtg                                                                                                                                       >  1:239326/1‑149 (MQ=255)
                atCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGAAGATGGCAATGATGCTGTCGGGGTTGGCCGGGCTTCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTAGAAGTCAGGTAGgtgt                                                                                                                                      >  8:145180/1‑149 (MQ=255)
                   gACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCtt                                                                                                                                   >  2:207978/1‑149 (MQ=255)
                         gttgtAGGCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGa                                                                                                                             <  3:33748/149‑1 (MQ=255)
                          ttgtAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGaa                                                                                                                            <  5:213299/149‑1 (MQ=255)
                          ttgtAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGa                                                                                                                             >  6:277259/1‑148 (MQ=255)
                             tAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTTGAAGTCAGGTAGGTGTCTTCGCCGAAAcc                                                                                                                         >  2:23660/1‑149 (MQ=255)
                             tAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTTGAAGTCAGGTAGGTGTCTTCGCCGAAAcc                                                                                                                         >  6:383338/1‑149 (MQ=255)
                                      ccgccTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTGGGCCGGGCTGCTGCCCGGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTtgt                                                                                                                <  6:264958/149‑1 (MQ=255)
                                               aCGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGccc                                                                                                         >  1:175743/1‑147 (MQ=255)
                                               aCGGCGCCATACAAGGCGAAATGCTTGACGATGGACATGATGCTGTCGGGGTTGGCCGGGCTGCTGCACTGGAACGAGCGCACCATTACCTGGCCCATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGccc                                                                                                         >  2:126934/1‑147 (MQ=255)
                                                     ccATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCgg                                                                                                 >  6:85750/1‑149 (MQ=255)
                                                     ccATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCgg                                                                                                 <  1:9450/149‑1 (MQ=255)
                                                            aGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACgg                                                                                          >  1:20542/1‑149 (MQ=255)
                                                            aGGCGAAATGCTTGACGATGGCCATGATGCGGGCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACgg                                                                                          <  5:126923/149‑1 (MQ=255)
                                                             ggCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACggg                                                                                         >  6:267696/1‑149 (MQ=255)
                                                                          aCGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATccacca                                                                             <  5:267696/148‑1 (MQ=255)
                                                                                    tgatgCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAAc                                                                  <  5:85750/149‑1 (MQ=255)
                                                                                    tgatgCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAAc                                                                  <  2:123717/149‑1 (MQ=255)
                                                                                              ggggTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCa                                                        <  3:369001/149‑1 (MQ=255)
                                                                                                ggTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCTCCATTGGCGCGAACGTCATGTCCAgg                                                      <  1:198931/149‑1 (MQ=255)
                                                                                                       cGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCgg                                               <  5:277259/149‑1 (MQ=255)
                                                                                                       cGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCgg                                               <  3:323909/149‑1 (MQ=255)
                                                                                                           ctgctgCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTTGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGa                                           <  1:23660/149‑1 (MQ=255)
                                                                                                           ctgctgCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTTGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGa                                           <  5:383338/149‑1 (MQ=255)
                                                                                                           ctgctgCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGa                                           <  1:276258/149‑1 (MQ=255)
                                                                                                              ctgcCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAg                                          <  2:20542/147‑1 (MQ=255)
                                                                                                               tgcCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGa                                           <  6:169707/145‑1 (MQ=255)
                                                                                                               tgcCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGa                                           >  5:169707/1‑145 (MQ=255)
                                                                                                                 cccTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTc                                     >  6:170389/1‑149 (MQ=255)
                                                                                                                 cccTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTc                                     >  2:298483/1‑149 (MQ=255)
                                                                                                                   cTGGAACGAGCGCACCATTACCTGGCCGATGTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTc                                     >  1:223700/1‑148 (MQ=255)
                                                                                                                    tGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGAt                                  >  4:207384/1‑149 (MQ=255)
                                                                                                                          gAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCgcggc                            >  6:293520/1‑149 (MQ=255)
                                                                                                                                aCCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCGTCGATCGCGGAGGTGCg                      >  8:346321/1‑149 (MQ=255)
                                                                                                                                        cTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGACACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCTCTATCCACAATTGGCGCGAAAGTCATGTCCAGGGCGTCGGCCGAGGCTACGATCGCGGCGCTGCGGCCGACCt              >  7:164112/1‑149 (MQ=255)
                                                                                                                                           gCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTTg            <  1:207978/148‑1 (MQ=255)
                                                                                                                                           gCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTTg            <  2:175743/148‑1 (MQ=255)
                                                                                                                                              gATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTgggaca        >  5:318916/1‑146 (MQ=255)
                                                                                                                                                    tcGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCtt                                      >  4:354010/1‑114 (MQ=255)
                                                                                                                                                    tcGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCtt                                      <  3:354010/114‑1 (MQ=255)
                                                                                                                                                    tcGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTTg            <  3:353996/140‑1 (MQ=255)
                                                                                                                                                    tcGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTTg            >  4:353996/1‑140 (MQ=255)
                                                                                                                                                     cGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTTGGCGACGGCCt  >  4:175806/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTTCGGCAGGCTCATATCGACCGTGTTGTAGTCGCGCCCGCCTTCCACGGCGCCATACAAGGCGAAATGCTTGACGATGGCCATGATGCTGTCGGGGTTGGCCGGGCTGCTGCCCTGGAACGAGCGCACCATTACCTGGCCGATTCTCGAAGTCAGGTAGGTGTCTTCGCCGAAACCTTCGCTTGTGCGGCCCCAGCGCGGGTCACGGGCGATATCCACCATTGGCGCGAACGTCATGTCCAGGGCGTCGGCCGAGGCTTCGATCGCGGCGGTGCGGCCGACCTTGGCGACGGCCT  >  NC_002947/1600505‑1600801

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: