Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,460,565 G→A 6.4% F300F (TTC→TTT algJ ← alginate O‑acetylase AlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,460,5650GA 6.4% 187.3 / 2.3 62F300F (TTC→TTTalgJalginate O‑acetylase AlgJ
Reads supporting (aligned to +/- strand):  ref base G (25/33);  new base A (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCGGGGCCGGCAGCAGGT  >  NC_002947/1460419‑1460713
                                                                                                                                                   |                                                                                                                                                      
taggggAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAACGCATCGCATTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAgg                                                                                                                                                        >  1:94497/3‑147 (MQ=255)
       aTTCCCACACCACCACTTGTGGTGCGGCGTTCTTTAAGGCATCGCTTTGCTGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCTTAGTTGGCGACGTCGCTGGGCTGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCa                                                                                                                                                <  5:165671/148‑1 (MQ=255)
           ccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGc                                                                                                                                            <  1:194791/148‑1 (MQ=255)
           ccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCg                                                                                                                                           >  1:128236/1‑149 (MQ=255)
            cacaCCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCAAGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCgg                                                                                                                                          >  2:237981/1‑149 (MQ=255)
                ccaccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTg                                                                                                                                      <  5:78067/149‑1 (MQ=255)
                    cacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCTGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGgctc                                                                                                                                  <  4:326017/149‑3 (MQ=255)
                    cacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCAGCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGgcgc                                                                                                                                  <  3:363908/149‑1 (MQ=255)
                       ttGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGt                                                                                                                               <  6:228579/149‑1 (MQ=255)
                         gtggtgCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTa                                                                                                                             <  2:128236/148‑1 (MQ=255)
                              gcggcgTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGt                                                                                                                       >  1:76138/1‑149 (MQ=255)
                              gaggcgTTCTTGACGGATTCGCTTTGCATGTACTTGATCATCCGGATCAGCGGGCTATTTCCGTGTTTGGGGTAGTTGGCGACGTCGCTGCGCAGTGCCTGCTGCAGCGCTTCAATGAATTTCCAGTGCTGGTTGGCGCTGTAGCTgg                                                                                                                        <  2:94497/146‑1 (MQ=255)
                               cggcgTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTACTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGt                                                                                                                       <  3:238288/148‑1 (MQ=255)
                                ggcgTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGTCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGc                                                                                                                     <  8:273955/149‑1 (MQ=255)
                                      cTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCCGCGCCAAGGAAGTTCCAGTGCGGGTTgg                                                                                                                                     <  1:29258/129‑1 (MQ=255)
                                      cTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTgg                                                                                                                                     >  2:29258/1‑128 (MQ=255)
                                         gAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCa                                                                                                               <  5:101558/146‑1 (MQ=255)
                                         gAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAg                                                                                                              <  2:296619/147‑1 (MQ=255)
                                             gCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTCCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAAc                                                                                                        >  5:385376/1‑149 (MQ=255)
                                             gCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCAGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAAc                                                                                                        >  2:359197/1‑149 (MQ=255)
                                               aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCg                                                                                                      <  7:210477/149‑1 (MQ=255)
                                               aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCg                                                                                                      >  3:211505/1‑149 (MQ=255)
                                               aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGAAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCg                                                                                                      >  4:231223/1‑149 (MQ=255)
                                               aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGAAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCg                                                                                                      >  4:226334/1‑149 (MQ=255)
                                                     ttGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCt                                                                                                >  1:187816/1‑149 (MQ=255)
                                                     ttGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCt                                                                                                >  2:324562/1‑149 (MQ=255)
                                                        cAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGAAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGc                                                                                              <  3:231223/148‑1 (MQ=255)
                                                        cAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGAAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGc                                                                                              <  3:226334/148‑1 (MQ=255)
                                                           gTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGAt                                                                                                  <  4:30487/141‑1 (MQ=255)
                                                           gTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGAt                                                                                                  >  3:30487/1‑141 (MQ=255)
                                                              cTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCg                                                                                        <  7:227055/148‑1 (MQ=255)
                                                              cTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCg                                                                                        <  8:285380/148‑1 (MQ=255)
                                                                 gAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCg                                                                                     >  1:182058/1‑148 (MQ=255)
                                                                      tCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGaacaac                                                                              <  2:76138/149‑1 (MQ=255)
                                                                             cagcGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGc                                                                       >  2:99863/1‑149 (MQ=255)
                                                                                cGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACACCGCGTCGCCgg                                                                    >  8:39910/1‑149 (MQ=255)
                                                                                    ccATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGTGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTc                                                                >  7:359990/1‑149 (MQ=255)
                                                                                    ccATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGt                                                                 >  2:490/1‑148 (MQ=255)
                                                                                     cATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCa                                                               <  7:128998/149‑1 (MQ=255)
                                                                                          ccGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTc                                                          >  5:367209/1‑149 (MQ=255)
                                                                                             tCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGc                                                       <  5:35619/149‑1 (MQ=255)
                                                                                               ttCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAAcc                                                                                                       >  6:147449/1‑100 (MQ=255)
                                                                                               ttCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAAcc                                                                                                       <  5:147449/100‑1 (MQ=255)
                                                                                                 cGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCgg                                                                                                     >  2:10335/1‑100 (MQ=255)
                                                                                                 cGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCgg                                                                                                     <  1:10335/100‑1 (MQ=255)
                                                                                                 cGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATc                                                    <  1:359197/148‑1 (MQ=255)
                                                                                                  ggCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGa                                                  >  3:188741/1‑149 (MQ=255)
                                                                                                         ttGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGAccggccg                                           <  2:182058/149‑1 (MQ=255)
                                                                                                            gCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGggt                                        <  2:187816/149‑1 (MQ=255)
                                                                                                             cGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCg                                                         <  4:210262/131‑1 (MQ=255)
                                                                                                             cGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCg                                                         >  3:210262/1‑131 (MQ=255)
                                                                                                              gACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGggtgg                                      <  4:211505/149‑1 (MQ=255)
                                                                                                                      tGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTc                              >  3:98134/1‑149 (MQ=255)
                                                                                                                       gcgcAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCt                             <  1:99863/149‑1 (MQ=255)
                                                                                                                       gcgcAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCt                             <  1:490/149‑1 (MQ=255)
                                                                                                                        cgcAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCt                             <  4:188741/148‑1 (MQ=255)
                                                                                                                           aaggCCTGCTGCATCGCGCCCAGGAAGTTCCACTGCCGGGTTGCCCTTTATCCTGTTCCAACCACTTCAAACGGGATATGCCTGTCGGCCAACACACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCa                          <  1:296803/147‑1 (MQ=255)
                                                                                                                             ggCCTGCTGCAGTGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGt                       <  8:359990/149‑1 (MQ=255)
                                                                                                                                 tgctgcAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTc                   >  4:196557/1‑149 (MQ=255)
                                                                                                                                     gcAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCggg               >  1:103340/1‑149 (MQ=255)
                                                                                                                                      cAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCgggg              >  6:305823/1‑149 (MQ=255)
                                                                                                                                        gcgcgcCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTGAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGAccggccgg                                          >  1:63677/1‑120 (MQ=255)
                                                                                                                                        gcgcgcCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGAccggccgg                                          <  2:63677/119‑1 (MQ=255)
                                                                                                                                                  ggAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCGGGGCCGGCAGCAgg   <  3:196557/149‑1 (MQ=255)
                                                                                                                                                   gAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCGGGGCCGGCAGCAGGt  <  1:237981/149‑1 (MQ=255)
                                                                                                                                                   gAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCGGGGCCGGCAGCAGGt  >  2:204462/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCGGGGCCGGCAGCAGGT  >  NC_002947/1460419‑1460713

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: