Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,602,521 A→G 7.4% Q123R (CAG→CGG) ‡ PP_5469 → SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,602,5210AG 7.4% 188.2 / 4.5 67Q123R (CAG→CGG) ‡PP_5469SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base G (3/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

CAACCTGCCGCTGGACCAGGCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTCGCGCAAGTGGAAGCAGA  >  NC_002947/1602375‑1602667
                                                                                                                                                  |                                                                                                                                                  
cAACCTGCCGCTGGACCAGGCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAgg                                                                                                                                                  <  4:295035/149‑1 (MQ=255)
cAACCTGCCGCTGGACCAGGCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAgg                                                                                                                                                  <  3:309320/149‑1 (MQ=255)
     tGCCGCTGGACCAGGCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGCCACcg                                                                                                                                             >  7:98794/1‑149 (MQ=255)
       ccGCTGCACCAGGCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGAAGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAACGGCAAACCGATGCTCCAGCGCCACGGTGCCTGGGCACcgt                                                                                                                                            <  4:34628/148‑2 (MQ=255)
        ccccggACCACGCCCACGTGCTGGTGCCCGAACCGCCCCCGGACAGCTTCGCCCATCCCCACCTGTTCGATGCTGACGCCGAAGTCCAGCCGACCTTCAGCAAGAAGGGCAAACCGATGCGCCACCGCCACGCTGCCGAGGCACcgcgc                                                                                                                                          <  8:98794/145‑1 (MQ=255)
          cTGGACCAGGCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGt                                                                                                                                        <  2:94867/149‑1 (MQ=255)
                  ggCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTAGCCGATg                                                                                                                                >  1:122473/1‑149 (MQ=255)
                    cccAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATgcg                                                                                                                              >  8:234286/1‑149 (MQ=255)
                     ccAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATgcgg                                                                                                                             >  4:157539/1‑149 (MQ=255)
                                  gCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGa                                                                                                                >  4:329295/1‑149 (MQ=255)
                                   ccGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGa                                                                                                                <  5:328230/148‑1 (MQ=255)
                                           ctgcCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCgag                                                                                                         >  8:313440/1‑147 (MQ=255)
                                            tgcCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCgag                                                                                                         >  7:84640/1‑146 (MQ=255)
                                              ccGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGagaag                                                                                                      >  7:40638/1‑147 (MQ=255)
                                               cGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAATAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGaagaag                                                                                                   <  7:161556/149‑1 (MQ=255)
                                                    aGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTa                                                                                              >  3:316032/1‑149 (MQ=255)
                                                          cGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGa                                                                                        <  3:157539/149‑1 (MQ=255)
                                                          cGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTAGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGa                                                                                        <  2:122473/149‑1 (MQ=255)
                                                           gCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAg                                                                                       >  1:63006/1‑149 (MQ=255)
                                                                   cACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGATCTGTcgcg                                                                               <  7:234286/149‑1 (MQ=255)
                                                                    aCCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCTTTATCTGGAGCTGTCGCgg                                                                              >  8:74450/1‑149 (MQ=255)
                                                                            gATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCt                                                                      >  6:242843/1‑149 (MQ=255)
                                                                                 tGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTgcgc                                                                  >  1:278833/1‑148 (MQ=255)
                                                                                 tGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTgcgc                                                                  >  5:281599/1‑148 (MQ=255)
                                                                                   aCGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGa                                                                <  5:242843/148‑1 (MQ=255)
                                                                                   aCGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGAc                                                               <  4:114368/149‑1 (MQ=255)
                                                                                    cGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGAcc                                                              <  7:313440/149‑1 (MQ=255)
                                                                                       cGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGc                                                                                                  >  7:1338/1‑110 (MQ=255)
                                                                                       cGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGc                                                                                                  <  8:1338/110‑1 (MQ=255)
                                                                                        gAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAgaagaa                                                                                                    <  4:113681/107‑1 (MQ=255)
                                                                                        gAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAgaagaa                                                                                                    >  3:113681/1‑107 (MQ=255)
                                                                                          aGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGgtg                                                        <  6:281599/149‑1 (MQ=255)
                                                                                          aGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGgt                                                         <  2:278833/148‑1 (MQ=255)
                                                                                                gCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGAcc                                                   >  8:304520/1‑148 (MQ=255)
                                                                                                gCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCg                                                  >  3:35296/1‑149 (MQ=255)
                                                                                                 cTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGa                                                 <  3:261438/149‑1 (MQ=255)
                                                                                                    aCCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGc                                                                             <  5:226312/118‑1 (MQ=255)
                                                                                                    aCCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGc                                                                             >  6:226312/1‑118 (MQ=255)
                                                                                                     ccTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGca                                             <  8:40638/149‑1 (MQ=255)
                                                                                                        tCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGcac                                            >  7:107802/1‑147 (MQ=255)
                                                                                                         cAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGcaca                                           <  5:318848/147‑1 (MQ=255)
                                                                                                            cAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGgcg                                      >  2:154261/1‑149 (MQ=255)
                                                                                                             aagaagGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTCTCTGCGCGACCTCGGTGTGACCGACGCACAAGgcg                                      >  6:99343/1‑148 (MQ=255)
                                                                                                                   ggCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGAt                               >  4:282590/1‑149 (MQ=255)
                                                                                                                   ggCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGAt                               >  4:64443/1‑149 (MQ=255)
                                                                                                                    gCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCAGAAGGCGCACTGATc                              <  7:304520/149‑1 (MQ=255)
                                                                                                                    gCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATc                              <  4:35296/149‑1 (MQ=255)
                                                                                                                    gCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATc                              <  8:84640/149‑1 (MQ=255)
                                                                                                                      aaaCCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTc                                                                                  >  8:22724/1‑95 (MQ=255)
                                                                                                                      aaaCCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTc                                                                                  <  7:22724/95‑1 (MQ=255)
                                                                                                                      aaaCCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATccc                            <  6:180417/149‑1 (MQ=255)
                                                                                                                       aaCCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCg                           >  2:126376/1‑149 (MQ=255)
                                                                                                                         ccGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTc                         <  1:154261/149‑1 (MQ=255)
                                                                                                                         ccGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTc                         <  8:107802/149‑1 (MQ=255)
                                                                                                                                tCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACaa                                          >  6:148580/1‑125 (MQ=255)
                                                                                                                                tCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACaa                                          <  5:148580/125‑1 (MQ=255)
                                                                                                                                tCCAGCGCCACGGTGCCCAGCCACCGCGCGTGGCCGATGCGGGCCGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGAGTGACCGCCGCACAAGGCGCACTCATCCCGTCCATGTcg                  >  8:260540/1‑149 (MQ=255)
                                                                                                                                 ccAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGTGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTcgc                 >  1:113854/1‑149 (MQ=255)
                                                                                                                                       gcgcGGTGCCTGGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTCGCGCAAGt           >  6:83551/4‑149 (MQ=255)
                                                                                                                                       ccACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTCGCGCAAGt           <  4:316032/149‑1 (MQ=255)
                                                                                                                                        cACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTCGCGCAAGTg          <  7:74450/149‑1 (MQ=255)
                                                                                                                                         gcGGTGCCTGGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGt                    <  6:14161/137‑1 (MQ=255)
                                                                                                                                         gcGGTGCCTGGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGt                    >  5:14161/2‑138 (MQ=255)
                                                                                                                                         gcGGTGCCTGGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTCGCGCAAGTg          >  6:232256/2‑148 (MQ=255)
                                                                                                                                               cccAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCAt                      <  6:20791/130‑1 (MQ=255)
                                                                                                                                               cccAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCAt                      >  5:20791/1‑130 (MQ=255)
                                                                                                                                                 cAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTCGCGCAAGTGGAAGCAGa  >  4:33518/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAACCTGCCGCTGGACCAGGCCCAGGTGCTGGTTGCCGAACTGCTGCCGGACAGCTTCCGCAATGCCCACCTGTTCGATGCTGACGGCGAAGTGCAGCTGACCTTCAGCAAGAAGGGCAAACCGATGCTCCAGCGCCACGGTGCCCAGGCACCGCGCGTGGCCGATGCGGGCAGCGGGCATGACCGCGAGAAGAAGCGTTATCTGGAGCTGTCGCGGCCGTTTCTGCGCGACCTCGGTGTGACCGACGCACAAGGCGCACTGATCCCGTCCATGTCGCGCAAGTGGAAGCAGA  >  NC_002947/1602375‑1602667

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