Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,469,670 A→G 6.7% V158A (GTC→GCC)  algK ← alginate biosynthesis protein AlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,469,6700AG 6.7% 174.3 / 2.8 60V158A (GTC→GCC) algKalginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

GGTGGTGTTCAGCGCGGCCTTGCAGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAG  >  NC_002947/1469522‑1469812
                                                                                                                                                    |                                                                                                                                              
ggtggtGTTCAGCGCGGCCTTGCAGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGa                                                                                                                                                >  5:213669/1‑149 (MQ=255)
      gtTCAGCGCGGCCTTGCAGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTgg                                                                                                                                          <  4:42726/149‑1 (MQ=255)
          aGCGCGGCCTTGCAGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGa                                                                                                                                       >  1:99462/1‑148 (MQ=255)
           gcgcGGCCTTGCAGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCAGCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGaaa                                                                                                                                     <  2:148009/149‑1 (MQ=255)
               ggCCTTGCAGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAActct                                                                                                                                 <  2:251076/149‑1 (MQ=255)
                       aGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAg                                                                                                                         >  5:6710/1‑149 (MQ=255)
                              ttCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCt                                                                                                                       <  3:84648/144‑1 (MQ=255)
                               tCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAgg                                                                                                                   >  6:19865/1‑147 (MQ=255)
                                ccACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGa                                                                                                                                                >  2:37179/1‑117 (MQ=255)
                                ccACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGa                                                                                                                                                <  1:37179/117‑1 (MQ=255)
                                 cACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAACACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTaga                                                                                                               <  3:102747/149‑1 (MQ=255)
                                 cACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAACACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTa                                                                                                                 >  4:102747/1‑147 (MQ=255)
                                           aGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAg                                                                                                     <  4:158883/149‑1 (MQ=255)
                                             gTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCg                                                                                                   >  1:248832/1‑149 (MQ=255)
                                                cTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGaa                                                                                                >  7:156253/1‑149 (MQ=255)
                                                cTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCACGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGCTCTGCTGCTGCGCGTTGACCTTCGGGAAACTCTGCGGCTAGCTGAGGTAGAGCATTGcacacgga                                                                                                 >  8:52483/1‑141 (MQ=255)
                                                 tGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAAt                                                                                               >  2:169474/1‑149 (MQ=255)
                                                   gTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCAtt                                                                                                          >  7:112052/1‑136 (MQ=255)
                                                   gTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCAtt                                                                                                          <  8:112052/136‑1 (MQ=255)
                                                     cATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAgg                                                                                           >  6:221512/1‑149 (MQ=255)
                                                                tGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTg                                                                                >  8:92998/1‑149 (MQ=255)
                                                                  ggTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTgg                                                                               <  5:258194/148‑1 (MQ=255)
                                                                     gCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCt                                                                           >  3:16986/1‑149 (MQ=255)
                                                                           cagcagCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTcg                                                                     >  6:104624/1‑149 (MQ=255)
                                                                            agcagcACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTcgc                                                                    >  8:83009/1‑149 (MQ=255)
                                                                                gcaCCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGGCCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCCGCGGAATCAGGGTATTGCTTTGGCCCGGTTTCGCGGcc                                                                >  2:120259/1‑149 (MQ=255)
                                                                                gcaCCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGGCCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGcc                                                                >  2:124385/1‑149 (MQ=255)
                                                                                gcaCCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGcc                                                                >  3:92283/1‑149 (MQ=255)
                                                                                    cTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCt                                                            >  4:142391/1‑149 (MQ=255)
                                                                                     tGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCtt                                                           >  4:249629/1‑149 (MQ=255)
                                                                                       ggCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGa                                                         <  6:213669/149‑1 (MQ=255)
                                                                                       ggCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCACGGCCTGCTTGa                                                         <  8:156253/149‑1 (MQ=255)
                                                                                        gCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGGTGGCCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAg                                                        <  1:124385/149‑1 (MQ=255)
                                                                                        gCCAGGCCCGCTTCCGGGGTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGGCCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGGATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAg                                                        <  1:120259/149‑1 (MQ=255)
                                                                                          cAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCa                                                      <  2:99462/149‑1 (MQ=255)
                                                                                             gCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAgg                                                    >  6:43743/1‑148 (MQ=255)
                                                                                             gCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAgg                                                    <  5:43743/148‑1 (MQ=255)
                                                                                               ccGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGtt                                                 >  8:241046/1‑149 (MQ=255)
                                                                                               ccGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGtt                                                 <  5:19865/149‑1 (MQ=255)
                                                                                                 gCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTc                                               <  7:241046/149‑1 (MQ=255)
                                                                                                  cTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCg                                              <  3:142391/149‑1 (MQ=255)
                                                                                                   ttCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCgg                                             >  5:37684/1‑149 (MQ=255)
                                                                                                       gggTTGCCGGCGGCGTGCCACTGGTCGATCTGCTGCTGGGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTc                                         <  7:92998/149‑1 (MQ=255)
                                                                                                         gTTGGCGTCGGCGCGCCTGTGGTCGATCTGCTGCTTCGCGTTTAGCTAGGGGAAACTCTGCCGGTAGCGGAGGTAGAGGATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCATGGTTTCGGCTTCCt                                       <  6:205479/149‑1 (MQ=255)
                                                                                                            gCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCa                                     <  5:221512/148‑1 (MQ=255)
                                                                                                                   gcgcgcCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGc                             <  1:169474/149‑1 (MQ=255)
                                                                                                                    cgcgcCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGcc                            >  4:127042/1‑149 (MQ=255)
                                                                                                                     gcgcCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCt                           >  1:241032/1‑149 (MQ=255)
                                                                                                                       gcCACTGGTCGTTCTGCTGCTGCGCGTGGACCTTGGGGATACTCTGCGGGGAGCTGAGGTAGAGCATTGCCAGCGGTATCAGGGTATTGCGTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTgg                         <  7:52483/149‑1 (MQ=255)
                                                                                                                       gcCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTgg                         <  4:92283/149‑1 (MQ=255)
                                                                                                                       gcCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTgg                         <  6:6710/149‑1 (MQ=255)
                                                                                                                          aCTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTg                      <  5:104624/149‑1 (MQ=255)
                                                                                                                              gTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCAc                                    <  2:17148/131‑1 (MQ=255)
                                                                                                                              gTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCAc                                    >  1:17148/1‑131 (MQ=255)
                                                                                                                                      gctgctgcGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTg          >  5:247577/1‑149 (MQ=255)
                                                                                                                                       ctgctgcGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCgg              <  4:16986/144‑1 (MQ=255)
                                                                                                                                       ctgctgcGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGc             <  3:127042/145‑1 (MQ=255)
                                                                                                                                       ctgctgcGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCtt           <  6:37684/147‑1 (MQ=255)
                                                                                                                                       ctgctgcGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTgg         <  7:83009/149‑1 (MQ=255)
                                                                                                                                       ctgctgcGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTg          <  2:248832/148‑1 (MQ=255)
                                                                                                                                          ctgcGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGcc       <  6:247577/148‑1 (MQ=255)
                                                                                                                                                gTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACcag  >  1:216583/1‑147 (MQ=255)
                                                                                                                                                gTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACcag  >  8:188414/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGTGGTGTTCAGCGCGGCCTTGCAGATCTTTTCCACTTCACCCAGGTGCTGGTCATAGGTGCCCTGGGTGCGGTACAGCAGCACCTGGGCCAGGCCCGCTTCCGGGTTGCCGGCGGCGCGCCACTGGTCGATCTGCTGCTGCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAG  >  NC_002947/1469522‑1469812

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