Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,468,805 T→C 5.4% Q446Q (CAA→CAG)  algK ← alginate biosynthesis protein AlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,468,8050TC 5.4% 220.4 / 2.1 74Q446Q (CAA→CAG) algKalginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base C (2/2);  total (37/37)
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.06e-01

CCGATGCCGGCTTTCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTGGCCCTCGCT  >  NC_002947/1468657‑1468951
                                                                                                                                                      |                                                                                                                                                  
ccGATGCCGGCTTTCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGt                                                                                                                                                    >  5:110414/1‑149 (MQ=255)
      ccGGCTTTCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTgg                                                                                                                                              <  4:160317/149‑1 (MQ=255)
         gCTTTCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGcc                                                                                                                                           >  4:266692/1‑149 (MQ=255)
          cTTTCACGAAGGGATTGAGCGTCATAGCGAGTCCGTCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAACTGCTggac                                                                                                                                            >  1:237472/1‑146 (MQ=255)
          cTTTCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAACTGCTGGAcc                                                                                                                                           <  2:237472/148‑1 (MQ=255)
           tttCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTg                                                                                                                                         >  4:91845/1‑149 (MQ=255)
             tCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAACTGCTGGACCAgg                                                                                                                                        <  3:334245/148‑1 (MQ=255)
              cACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGt                                                                                                                                       <  7:256394/148‑1 (MQ=255)
               aCGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGtt                                                                                                                                      <  8:73004/148‑1 (MQ=255)
                cGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTgc                                                                                                                                    <  6:110414/149‑1 (MQ=255)
                 gAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTTGCGTCGTGTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCATGCGCCAGCATGCCGCGTGCCCGGTTTTCCTGTTCCTGCAGTTGCTGGGCCTGGTTgc                                                                                                                                    >  1:268803/1‑149 (MQ=255)
                 gAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTgcg                                                                                                                                   >  6:328914/1‑149 (MQ=255)
                  aaGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTgcgc                                                                                                                                  <  6:222563/149‑1 (MQ=255)
                      gATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGAt                                                                                                                              <  4:349656/149‑1 (MQ=255)
                        ttGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGa                                                                                                                               >  6:169490/1‑146 (MQ=255)
                        ttGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGa                                                                                                                               <  5:169490/146‑1 (MQ=255)
                          gAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAg                                                                                                                           <  3:266692/148‑1 (MQ=255)
                           atcGTCATAGCGAGGCC‑TCACCTTCTACTTC‑TTAC‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGgt                                                                                                                         <  7:355825/146‑1 (MQ=255)
                             cGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGgtgt                                                                                                                       <  3:27848/149‑1 (MQ=255)
                                 aTAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGc                                                                                                                                                 <  3:88963/119‑1 (MQ=255)
                                 aTAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGc                                                                                                                                                 >  4:88963/1‑119 (MQ=255)
                                  tAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGtt                                                                                                                  >  5:367010/1‑149 (MQ=255)
                                         tCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTc                                                                                                           <  2:268803/149‑1 (MQ=255)
                                               aCCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAg                                                                                                        <  3:353149/147‑1 (MQ=255)
                                                 cGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAgctgc                                                                                                    >  6:49518/1‑149 (MQ=255)
                                                  gTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAgctgct                                                                                                   >  5:229255/1‑149 (MQ=255)
                                                  gTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAgctgct                                                                                                   <  8:225042/149‑1 (MQ=255)
                                                   tCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCt                                                                                                                                          >  1:203645/1‑109 (MQ=255)
                                                   tCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCt                                                                                                                                          <  2:203645/109‑1 (MQ=255)
                                                   tCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGctgctg                                                                                                  >  7:348359/1‑149 (MQ=255)
                                                       cttcttTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTgcag                                                                                               <  5:49518/148‑1 (MQ=255)
                                                       cttcttTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGACGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGATGGTCGAGCTGCTgcag                                                                                               >  7:301851/1‑148 (MQ=255)
                                                        ttcttTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAgc                                                                                                       <  8:284949/139‑1 (MQ=255)
                                                        ttcttTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAgc                                                                                                       >  7:284949/1‑139 (MQ=255)
                                                        ttcttTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTgcag                                                                                               >  1:55938/1‑147 (MQ=255)
                                                        ttcttTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTgcag                                                                                               >  6:184713/1‑147 (MQ=255)
                                                             ttt‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGtt                                                                                          <  1:155685/147‑1 (MQ=255)
                                                                  gtcgtcTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCt                                                                                     >  8:203041/1‑148 (MQ=255)
                                                                  gtcgtcTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCt                                                                                     <  7:86444/148‑1 (MQ=255)
                                                                           cAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTgcggcg                                                                           >  6:228565/1‑149 (MQ=255)
                                                                            aGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGt                                                                          >  6:225967/1‑149 (MQ=255)
                                                                                 ttCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGt                                                                     >  4:160028/1‑149 (MQ=255)
                                                                                   cGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGt                                                                          >  8:58500/1‑142 (MQ=255)
                                                                                   cGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGt                                                                          <  7:58500/142‑1 (MQ=255)
                                                                                   cGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTg                                                                   >  7:338647/1‑149 (MQ=255)
                                                                                      gggCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGcc                                                                >  5:11793/1‑149 (MQ=255)
                                                                                      gggCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGcc                                                                >  7:139278/1‑149 (MQ=255)
                                                                                       ggCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCt                                                               <  3:160028/149‑1 (MQ=255)
                                                                                         gcAGGGCGCGAATAGCACCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTgg                                                                                                             >  1:269853/2‑101 (MQ=255)
                                                                                         gcAGGGCGCGAATAGCACCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTgg                                                                                                             <  2:269853/100‑1 (MQ=255)
                                                                                         gcAGGGCGCGAATAGCACCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGc                                                                 <  1:134460/144‑1 (MQ=255)
                                                                                         gcAGGGCGCGAATAGCACCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGc                                                                 >  2:134460/2‑145 (MQ=255)
                                                                                             ggTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGa                                                           >  3:265654/1‑147 (MQ=255)
                                                                                                  tATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTCCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTgg                                                    >  8:63077/1‑149 (MQ=255)
                                                                                                    ttCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTggg                                                   <  6:11793/148‑1 (MQ=255)
                                                                                                    ttCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTgg                                                    <  3:91845/147‑1 (MQ=255)
                                                                                                    ttCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGc                                                  <  2:55938/149‑1 (MQ=255)
                                                                                                     tCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTg                                                                   <  7:301230/131‑1 (MQ=255)
                                                                                                     tCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTg                                                                   >  8:301230/1‑131 (MQ=255)
                                                                                                        cgcCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAac                                              <  5:328914/149‑1 (MQ=255)
                                                                                                          cccTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAacac                                            <  6:229255/149‑1 (MQ=255)
                                                                                                            cTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACAccc                                          <  8:338647/149‑1 (MQ=255)
                                                                                                               cgcCAGCCTGCAGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTAGCGCTCATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCAGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAgg                                       >  2:334008/1‑149 (MQ=255)
                                                                                                                   aGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGtt                                   <  5:228565/149‑1 (MQ=255)
                                                                                                                           gcgTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGtt                           <  5:184713/149‑1 (MQ=255)
                                                                                                                                cACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTgcggg                      >  5:83134/1‑149 (MQ=255)
                                                                                                                                         ccaggcccAGCAACTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCg             >  5:20895/7‑149 (MQ=255)
                                                                                                                                           tGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGGCAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTg           >  6:320573/1‑149 (MQ=255)
                                                                                                                                           tGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTg           >  5:19761/1‑149 (MQ=255)
                                                                                                                                             gTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTgg          <  5:337535/148‑1 (MQ=255)
                                                                                                                                             gTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTg           <  8:139278/147‑1 (MQ=255)
                                                                                                                                              tCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTGGcc        <  4:265654/149‑1 (MQ=255)
                                                                                                                                                  gcagTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTGGCCCTCg    >  5:232383/1‑149 (MQ=255)
                                                                                                                                                  gcagTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTGGCCCTCg    >  1:86369/1‑149 (MQ=255)
                                                                                                                                                     gTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTGGCCCTCGCt  >  1:321584/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCGATGCCGGCTTTCACGAAGGGATTGAGCGTCATAGCGAGTCC‑TCACCGTCTACTTCTTTTT‑CGTCGTCTTGCAGGGCTTCGAGGGCCAGGGTGCTATTCGCGCCCTGCGCCAGGCTGCCGCGTGCACGCTTTTCCTGGTCCAGCAGTTGCTGGGCCTGGTTGCGCTGATCAGGTGTCAGTTGCTGGTCGAGCTGCTGCAGCAGTTCGGCTGATTGCGGCGTCGGGTTGGCCTGCGACAGCTGGGCGAACACCCAGGCGTTGCCCGGTTGCGGGCGGATGCCGTGGCCCTCGCT  >  NC_002947/1468657‑1468951

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