Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 683,823 T→C 8.2% S10G (AGC→GGC)  PP_0585 ← MerR family metal‑sensitive transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947683,8230TC 8.2% 269.7 / 11.0 98S10G (AGC→GGC) PP_0585MerR family metal‑sensitive transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (51/38);  new base C (3/5);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCA  >  NC_002947/683675‑683968
                                                                                                                                                    |                                                                                                                                                 
gAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCt                                                                                                                                                   >  3:162318/1‑149 (MQ=255)
 aaaaGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTg                                                                                                                                                  >  8:110497/1‑149 (MQ=255)
      cccAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTgcggcg                                                                                                                                             <  5:10061/149‑1 (MQ=255)
           gTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGgcgggcgg                                                                                                                                        <  6:46447/149‑1 (MQ=255)
           gTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGgcgggcgg                                                                                                                                        >  3:321555/1‑149 (MQ=255)
           gTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTCGTGCTGAGCCCGCTGCGgcgggcgg                                                                                                                                        >  6:277783/1‑149 (MQ=255)
            tCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCGCagcgggctca                                                                                                                                       <  8:50431/149‑4 (MQ=255)
             cgcgGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCGCagcgggctcag                                                                                                                                      <  4:219674/149‑5 (MQ=255)
              gcgGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCGCagcgggctcagc                                                                                                                                     <  4:87213/149‑6 (MQ=255)
              gcgGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCGCagcgggctcag                                                                                                                                      <  6:433468/148‑5 (MQ=255)
               cgGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTg                                                                                                                                    <  4:321555/149‑1 (MQ=255)
                gggAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTgg                                                                                                                                   >  7:341872/1‑149 (MQ=255)
                    gcgcTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGa                                                                                                                               <  7:196853/149‑1 (MQ=255)
                      gcTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATg                                                                                                                             >  7:53486/1‑149 (MQ=255)
                       cTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCGCagcgggctcagc                                                                                                                                     <  1:395948/140‑6 (MQ=255)
                       cTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCGCagcgggctcagc                                                                                                                                     >  2:395948/1‑135 (MQ=255)
                          gatgaAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCa                                                                                                                         >  3:28376/1‑149 (MQ=255)
                             gaAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCt                                                                                                                                     <  1:220864/134‑1 (MQ=255)
                             gaAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCt                                                                                                                                     >  2:220864/1‑134 (MQ=255)
                               aGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTg                                                                                                                                                  >  6:59658/1‑119 (MQ=255)
                               aGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTg                                                                                                                                                  <  5:59658/119‑1 (MQ=255)
                                 gCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGcc                                                                                                                  >  3:400061/1‑149 (MQ=255)
                                  ccAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGcct                                                                                                                 >  3:43372/1‑149 (MQ=255)
                                   cAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGcctc                                                                                                                >  6:294041/1‑149 (MQ=255)
                                    aGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGcctcc                                                                                                               >  3:177946/1‑149 (MQ=255)
                                      gTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGcctcct                                                                                                              <  8:53486/148‑1 (MQ=255)
                                      gTTGTGCAAGTCTTGCGCATTGTACAGGAGGTAGCCTCTGTCGCTGCGAGTGGCGGTTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGTTGCTGAGACAGCTGCGGCGGGCGGACTGGCCGATGTTCATGGCGCCTCCTg                                                                                                             >  8:295264/1‑149 (MQ=255)
                                       ttGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTgg                                                                                                            <  4:315606/149‑1 (MQ=255)
                                          tgCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGTATCATCTTGGTGCTGAGCCCTCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCtgtgtg                                                                                                         <  7:295264/149‑1 (MQ=255)
                                            cAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGtgct                                                                                                       >  2:227534/1‑149 (MQ=255)
                                             aaGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGtgctg                                                                                                      >  5:237366/1‑149 (MQ=255)
                                             aaGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGtgctg                                                                                                      >  1:303381/1‑149 (MQ=255)
                                                  ttGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGTGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggt                                                                                                 >  6:72770/1‑149 (MQ=255)
                                                  ttGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCGTGGCGCCTCCTGGGTGCTGCtggt                                                                                                 >  3:287559/1‑149 (MQ=255)
                                                  ttGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTgg                                                                                                            <  3:385614/138‑1 (MQ=255)
                                                  ttGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTgg                                                                                                            >  4:385614/1‑138 (MQ=255)
                                                    gcgcCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGCGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggt                                                                                                 <  5:235426/147‑1 (MQ=255)
                                                     cgcCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggtgg                                                                                               <  2:417202/148‑1 (MQ=255)
                                                       ccTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGtt                                                                                            >  6:121848/1‑149 (MQ=255)
                                                        cTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGttt                                                                                           >  4:5731/1‑149 (MQ=255)
                                                          ggTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTc                                                                                          >  5:25202/1‑148 (MQ=255)
                                                           gTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGTTGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTcc                                                                                         <  7:110497/148‑1 (MQ=255)
                                                             aCAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCa                                                                                        >  6:395923/1‑147 (MQ=255)
                                                             aCAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCa                                                                                        <  5:395923/147‑1 (MQ=255)
                                                                ggCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGttt                                                                                   >  4:138208/1‑149 (MQ=255)
                                                                 gCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCAc                                                                                       >  3:8728/1‑144 (MQ=255)
                                                                 gCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCAc                                                                                       <  4:8728/144‑1 (MQ=255)
                                                                 gCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGttt                                                                                   <  5:294041/148‑1 (MQ=255)
                                                                   ggTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTaa                                                                                <  5:121848/149‑1 (MQ=255)
                                                                    gTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAg                                                                               >  5:376124/1‑149 (MQ=255)
                                                                    gTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGAGCTGCTGGTGGGTTTCCACGTTTTAAg                                                                               >  6:207669/1‑149 (MQ=255)
                                                                     tAGCCGCTGGCGCGGCGAGTGGCGGGTGTGAGCAGGCCGATGGATTCGTAGTAACCGATCATCTTGGTGCTGACCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGCTTCCACGTTTTAAGc                                                                              <  4:43372/149‑1 (MQ=255)
                                                                        ccGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGCGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATCGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCca                                                                            >  8:77202/1‑148 (MQ=255)
                                                                         cGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCGTGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCcaca                                                                          <  4:287559/149‑1 (MQ=255)
                                                                         cGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCca                                                                            >  6:18071/1‑147 (MQ=255)
                                                                          gCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAg                                                                         >  3:311924/1‑149 (MQ=255)
                                                                            tGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTg                                                                       >  5:353610/1‑149 (MQ=255)
                                                                              tCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGACCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGGTTTAAGCCACAGTGCa                                                                     >  4:258954/1‑149 (MQ=255)
                                                                                    cGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCt                                                                <  5:126119/148‑1 (MQ=255)
                                                                                    cGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCAGCTTGGTGCGGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCt                                                                <  8:341872/148‑1 (MQ=255)
                                                                                     gAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTga                                                              <  6:244743/149‑1 (MQ=255)
                                                                                      aGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTgac                                                             >  7:105185/1‑149 (MQ=255)
                                                                                          gCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGacgacg                                                         >  3:128395/1‑149 (MQ=255)
                                                                                             ggTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCt                                                       >  1:229528/1‑148 (MQ=255)
                                                                                               tttGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTga                                                              >  2:79422/1‑139 (MQ=255)
                                                                                               tttGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTga                                                              <  1:79422/139‑1 (MQ=255)
                                                                                               tttGAGCAGGCCGATGGATTCGTAGTAAAGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGa                                                     >  1:27166/1‑148 (MQ=255)
                                                                                                         ccGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAgg                                          <  6:25202/149‑1 (MQ=255)
                                                                                                          cGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAggg                                         >  3:402949/1‑149 (MQ=255)
                                                                                                               gATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCAt                                     <  4:400061/148‑1 (MQ=255)
                                                                                                                aTTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATgg                                   >  1:158809/1‑149 (MQ=255)
                                                                                                                   cGTAGTAACGGATCATCTTGGTGCTGAGCGCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGAGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGAGGCTGACACTGGACAGGGCCATggcg                                 <  5:207669/148‑1 (MQ=255)
                                                                                                                    gtagtaACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATggcgg                                <  6:237366/148‑1 (MQ=255)
                                                                                                                    gtagtaACGCATCAGCTTGGGGCTGAGCCCGCGGCGGCGGGCGGCCTGGCCGATGTTCATGCCGCCTCCTGGGTGCTGCTGGTGGGTTTCCAGGTTTTAAGCCACAGATCATTGCTGACGACGCTGACACTGGACAGGGCCATggcggc                               <  5:277783/149‑1 (MQ=255)
                                                                                                                       gtaACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTCTTAAGCCACAGTGCATTGGTGACGACGCTGACACTGGACAGGGCCATGGCGGCTcc                            >  2:281038/1‑149 (MQ=255)
                                                                                                                        taACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTcc                            <  4:28376/148‑1 (MQ=255)
                                                                                                                        taACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCt                           <  4:177946/149‑1 (MQ=255)
                                                                                                                        taACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCt                           >  8:297450/1‑149 (MQ=255)
                                                                                                                            ggATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGc                         >  5:423773/1‑147 (MQ=255)
                                                                                                                              atcatcTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAg                      >  1:420444/1‑148 (MQ=255)
                                                                                                                                 atcTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTATGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTAcc                  >  3:205188/1‑149 (MQ=255)
                                                                                                                                 atcTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTAcc                  >  2:414651/1‑149 (MQ=255)
                                                                                                                                    ttGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCATCCTGGGTGCTGCTGGTGGGATTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCACCTGCCAGTACCggg               >  5:442590/1‑149 (MQ=255)
                                                                                                                                      ggTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCggg               <  4:162318/147‑1 (MQ=255)
                                                                                                                                      ggTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGt              <  6:376124/148‑1 (MQ=255)
                                                                                                                                       gTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGtt             <  2:229528/148‑1 (MQ=255)
                                                                                                                                         gCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAg          <  3:306310/149‑1 (MQ=255)
                                                                                                                                         gCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAg          >  5:23474/1‑149 (MQ=255)
                                                                                                                                         aggccgCCCGCCGCAGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACACGGCCATGGCGGCTCCTGCCAGTACCGGGGTGAg          >  4:107013/6‑149 (MQ=255)
                                                                                                                                          cTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAg          <  1:414651/148‑1 (MQ=255)
                                                                                                                                          cTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAg          <  3:138208/148‑1 (MQ=255)
                                                                                                                                          cTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGGTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAgg         <  3:258954/149‑1 (MQ=255)
                                                                                                                                            gAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGt        >  6:92532/1‑148 (MQ=255)
                                                                                                                                            gAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGt        >  7:67530/1‑148 (MQ=255)
                                                                                                                                            ccgCCCGCCGCAGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGt        >  2:418872/3‑148 (MQ=255)
                                                                                                                                               cccGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTGCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACCGGGCCATGGCGGCTCCTGCCCGTACCGGGTTGAGGTaagc    >  6:5588/1‑147 (MQ=255)
                                                                                                                                                 cGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCa  <  5:92532/149‑1 (MQ=255)
                                                                                                                                                  gCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCggg               >  2:292141/1‑135 (MQ=255)
                                                                                                                                                  gCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCggg               <  1:292141/135‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCA  >  NC_002947/683675‑683968

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