Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,469,808 A→G 5.6% V112A (GTG→GCG)  algK ← alginate biosynthesis protein AlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,469,8080AG 5.6% 207.5 / 2.3 71V112A (GTG→GCG) algKalginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base G (2/2);  total (33/38)
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.80e-01

CGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAG  >  NC_002947/1469663‑1469953
                                                                                                                                                 |                                                                                                                                                  
cgcgTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGccac                                                                                                                                                   <  5:116096/147‑1 (MQ=255)
 gcgTTGACCTTGGGAAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACcag                                                                                                                                                >  8:20460/1‑149 (MQ=255)
      gACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGATGTCCGGCTTGGCCACCAGCAggc                                                                                                                                           >  6:13675/1‑149 (MQ=255)
       aCCGTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAggcg                                                                                                                                          >  2:196373/1‑149 (MQ=255)
        ccTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGccac                                                                                                                                                   <  3:271903/139‑1 (MQ=255)
        ccTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGccac                                                                                                                                                   >  4:271903/1‑139 (MQ=255)
         cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGccac                                                                                                                                                   >  6:18080/1‑138 (MQ=255)
         cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGccac                                                                                                                                                   <  5:18080/138‑1 (MQ=255)
         cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAggcggc                                                                                                                                        >  1:296362/1‑149 (MQ=255)
         cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAggcggc                                                                                                                                        >  4:204189/1‑149 (MQ=255)
         cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAggc                                                                                                                                           >  8:231163/1‑146 (MQ=255)
         cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAggc                                                                                                                                           <  7:231163/146‑1 (MQ=255)
         cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAcgcgcc                                                                                                                                        >  7:3456/1‑147 (MQ=255)
               gAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCGCCTGCTGGTGGCCAAgcc                                                                                                                                  <  8:33427/149‑3 (MQ=255)
               gAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGgc                                                                                                                                  >  1:162924/1‑149 (MQ=255)
                aaaCTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGgcg                                                                                                                                 <  3:204189/149‑1 (MQ=255)
                 aaCTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTTGCCACCAGCAGGCGGCCAAGgcgc                                                                                                                                <  4:238136/149‑1 (MQ=255)
                   ctctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGgcgcgc                                                                                                                             <  2:151209/149‑1 (MQ=255)
                              gCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTg                                                                                                                  <  2:162924/149‑1 (MQ=255)
                               cTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTgg                                                                                                                 >  1:47832/1‑149 (MQ=255)
                                 gAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGc                                                                                                                >  1:81447/1‑148 (MQ=255)
                                  aGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGa                                                                                                              >  5:311181/1‑149 (MQ=255)
                                    gTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCACGCGCCCAAGCGCGCCTTGGCCGCCTgctggtg                                                                                                            <  8:3456/149‑7 (MQ=38)
                                     tAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGCCCAAGCGCGCCTTGGCCGCCTgctggtg                                                                                                            <  6:248445/148‑7 (MQ=255)
                                       gagCATTGCCAGCGGAATCAGGTTATTGCTTAGGCCCTGTTTCGCGGCCTGCTTCAGCAGGGATTCGCCTTACTAACGTTCGGCCTGGCTGCTGTCCTGCTTGGCCACCAGCAGGCGGACAAGTCGCGCTTGGGAGCGAGCCGACGtcg                                                                                                         >  6:240893/1‑147 (MQ=255)
                                       gagCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTgg                                                                                                         >  5:73394/1‑149 (MQ=255)
                                        agCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTggc                                                                                                        >  1:12627/1‑149 (MQ=255)
                                         gCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTggcg                                                                                                       >  2:175208/1‑149 (MQ=255)
                                             tGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCg                                                                                                   >  4:57936/1‑149 (MQ=255)
                                              gCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCg                                                                                                   <  2:47832/148‑1 (MQ=255)
                                              gCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGGTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCg                                                                                                   <  5:222201/148‑1 (MQ=255)
                                               ccAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCg                                                                                                   >  2:164518/1‑147 (MQ=255)
                                                   cGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAgcgc                                                                                             <  6:2880/149‑1 (MQ=255)
                                                     gAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTggcggc                                                                                                     <  2:59800/139‑1 (MQ=255)
                                                     gAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTggcggc                                                                                                     >  1:59800/1‑139 (MQ=255)
                                                     gAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGa                                                                                           <  1:196373/149‑1 (MQ=255)
                                                          aGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGt                                                                                      <  6:73394/149‑1 (MQ=255)
                                                          aGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGt                                                                                      <  6:311181/149‑1 (MQ=255)
                                                             gTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGc                                                                                   <  2:81447/149‑1 (MQ=255)
                                                                ttGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                >  5:193542/1‑149 (MQ=255)
                                                                       ggCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTcc                                                                         >  4:62013/1‑149 (MQ=255)
                                                                        gCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCt                                                                        <  1:164518/149‑1 (MQ=255)
                                                                        gCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCATCACGTTCGGCCTGGGTGCTGTCAGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCt                                                                        >  8:202030/1‑149 (MQ=255)
                                                                         cccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCtt                                                                       >  1:22709/1‑149 (MQ=255)
                                                                         cccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCtt                                                                       >  5:85588/1‑149 (MQ=255)
                                                                          ccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCtt                                                                       >  8:232338/1‑148 (MQ=255)
                                                                          ccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCtt                                                                       >  3:312020/1‑148 (MQ=255)
                                                                          ccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTg                                                                      <  3:57936/149‑1 (MQ=255)
                                                                              tttCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGAttt                                                                  <  2:12627/149‑1 (MQ=255)
                                                                                tCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGc                                                                >  5:185191/1‑149 (MQ=255)
                                                                                   cgGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAg                                                             <  5:13675/149‑1 (MQ=255)
                                                                                   cgGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAg                                                             <  3:62013/149‑1 (MQ=255)
                                                                                    gtccTGCTTGAGCTGGGGTTCGGCTTCCTCTCGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGGAGGCGGCCAAGGCGCGCTTGTGCGCGTGGAGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAgg                                                            <  5:240893/147‑1 (MQ=255)
                                                                                    ggCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAgg                                                            <  7:232338/149‑1 (MQ=255)
                                                                                          cTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgc                                                      <  4:312020/149‑1 (MQ=255)
                                                                                           ttGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgcg                                                     <  3:229712/149‑1 (MQ=255)
                                                                                              aGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCggg                                                   <  6:185191/148‑1 (MQ=255)
                                                                                               gCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTg                                                 <  8:143065/149‑1 (MQ=255)
                                                                                               gCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTg                                                 >  7:143065/1‑149 (MQ=255)
                                                                                                  gggTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGcc                                              <  6:193542/149‑1 (MQ=255)
                                                                                                     tttCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCTTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCgg                                                    <  1:40392/140‑1 (MQ=255)
                                                                                                     tttCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCgg                                                    >  2:40392/1‑140 (MQ=255)
                                                                                                          gCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCAc                                            >  6:99838/1‑143 (MQ=255)
                                                                                                          gCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCAc                                            <  5:99838/143‑1 (MQ=255)
                                                                                                               cTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCTGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTcc                                  <  7:20460/148‑1 (MQ=255)
                                                                                                                  aCGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGc                                >  6:320968/1‑147 (MQ=255)
                                                                                                                  aCGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCa                              <  1:175208/149‑1 (MQ=255)
                                                                                                                    gTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAgg                            >  7:92812/1‑149 (MQ=255)
                                                                                                                     ttCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                           <  6:85588/149‑1 (MQ=255)
                                                                                                                             ggggGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTgc                   <  7:3126/149‑1 (MQ=255)
                                                                                                                              ggTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTgc                   <  2:296362/148‑1 (MQ=255)
                                                                                                                                     tCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGAGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCg             <  5:320968/147‑1 (MQ=255)
                                                                                                                                         ccTTGGCCGCCTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAg       >  1:125958/2‑149 (MQ=255)
                                                                                                                                           ttGGCCGCCTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGa                                       <  5:85300/115‑1 (MQ=255)
                                                                                                                                           ttGGCCGCCTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGa                                       >  6:85300/1‑115 (MQ=255)
                                                                                                                                              gCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAg  <  2:22709/149‑1 (MQ=255)
                                                                                                                                                caccaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCt                <  8:224674/133‑1 (MQ=255)
                                                                                                                                                caccaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCt                >  7:224674/1‑133 (MQ=255)
                                                                                                                                                caccaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGc            <  2:241661/137‑1 (MQ=255)
                                                                                                                                                caccaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGc            >  1:241661/1‑137 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAG  >  NC_002947/1469663‑1469953

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: