Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,356,318 G→C 8.0% G264G (GGG→GGC modC → molybdate ABC transporter ATP‑binding protein ModC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,356,3180GC 8.0% 277.2 / 12.2 88G264G (GGG→GGCmodCmolybdate ABC transporter ATP‑binding protein ModC
Reads supporting (aligned to +/- strand):  ref base G (43/38);  new base C (4/3);  total (47/41)
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.82e-01

ATCGCCAGCGGCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTGCGCGAGAGCTGCCC  >  NC_002947/4356169‑4356467
                                                                                                                                                     |                                                                                                                                                     
aTCGCCAGCGGCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCggg                                                                                                                                                       <  2:1308802/150‑1 (MQ=255)
 tCGCCAGCGGCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCCGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGCTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGc                                                                                                                                                      >  1:1234545/1‑150 (MQ=255)
  cGCCAGCGGCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGcc                                                                                                                                                     <  1:203224/150‑1 (MQ=255)
   gCCAGCGGCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCg                                                                                                                                                    >  2:858297/1‑150 (MQ=255)
      aGCGGCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCgctg                                                                                                                                                 >  1:1586172/1‑150 (MQ=255)
          gCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCtgctgc                                                                                                                                             <  2:1586172/150‑1 (MQ=255)
           cccGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGACCGGCAGCAGCGGGCCGCTGCTgcg                                                                                                                                            <  1:933330/150‑1 (MQ=255)
                  ggCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGCCGAAGCCGGCGTCGTGTTCGAAGGTCGGGTGGTTGGGCACGACCCGCATTAGGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCAcc                                                                                                                                     >  2:274287/1‑150 (MQ=255)
                    cGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCCGCAGCAGCGGCCCGCtgctgccgggca                                                                                                                                       <  1:1279105/146‑7 (MQ=255)
                    cGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCCGCAGCAGCGGCCCGCtgctgccgggca                                                                                                                                       >  2:1279105/1‑140 (MQ=255)
                     gAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGTTTACGACGCCGTCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCAc                                                                                                                                  <  2:1101635/150‑1 (MQ=255)
                         cccTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACgcgg                                                                                                                              <  2:686780/150‑1 (MQ=255)
                          ccTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACgcggc                                                                                                                             >  1:2012057/1‑150 (MQ=255)
                           cTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGAATCACCCACgcggcg                                                                                                                            >  2:1683431/1‑150 (MQ=255)
                              gCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAg                                                                                                                         >  2:577161/1‑150 (MQ=255)
                                 cacactGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTg                                                                                                                      >  2:1169009/6‑150 (MQ=255)
                                  gCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGAAGGGCGGTTGCGCACGCCGCTTATTACGGACTGCTGGACCTGTCACTGCCCTGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGa                                                                                                                     >  1:1741559/1‑150 (MQ=255)
                                  gCCTTGACCTGTCGCTGGCCCAGGGAGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGACTGCCCGGCAGCCGCCGGCCGCTGATGAGAAACACCCACGCGGCGCAGGTGa                                                                                                                     >  2:1697563/1‑150 (MQ=255)
                                   ccTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAg                                                                                                                         <  1:565132/145‑1 (MQ=255)
                                   ccTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAg                                                                                                                         >  2:565132/1‑145 (MQ=255)
                                     ttGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATgg                                                                                                                  >  2:864858/1‑150 (MQ=255)
                                      tGACCTGTCGCTGGCCCACGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATggg                                                                                                                 >  1:573171/1‑150 (MQ=255)
                                          cTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGgcagc                                                                                                             >  2:2367575/1‑150 (MQ=255)
                                           tGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGgcagca                                                                                                            >  2:1144062/1‑150 (MQ=255)
                                               gCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCt                                                                                                        <  2:2012057/150‑1 (MQ=255)
                                               gCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCCCATCACCTGCGCCGCGTGGGTGATGCgcagcagcgg                                                                                                        <  1:858297/150‑5 (MQ=37)
                                                  ggCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTgcg                                                                                                     >  1:721327/1‑150 (MQ=255)
                                                   gCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCtgctgc                                                                                                                                             <  2:1030684/109‑1 (MQ=255)
                                                   gCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCtgctgc                                                                                                                                             >  1:1030684/1‑109 (MQ=255)
                                                   gCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTgcgc                                                                                                    >  1:1265695/1‑150 (MQ=255)
                                                     ccAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGt                                                                                                  >  1:2409039/1‑150 (MQ=255)
                                                     ccAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGt                                                                                                  >  2:1126443/1‑150 (MQ=255)
                                                          gCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAgg                                                                                             <  1:577161/150‑1 (MQ=255)
                                                          aacacgacgCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAgg                                                                                             >  1:16752/4‑150 (MQ=255)
                                                               gacgCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAg                                                                                        <  2:299335/150‑1 (MQ=255)
                                                                    cGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGgcgc                                                                                                                           <  1:1101272/110‑1 (MQ=255)
                                                                    cGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGgcgc                                                                                                                           >  2:1101272/1‑110 (MQ=255)
                                                                     ggCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCgg                                                                                  <  1:2367575/150‑1 (MQ=255)
                                                                     ggCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCgg                                                                                  <  2:160382/150‑1 (MQ=255)
                                                                      gCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGCCCGCtgctgcc                                                                                                                                            <  1:2311398/91‑2 (MQ=255)
                                                                      gCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGCCCGCtgctgcc                                                                                                                                            >  2:2311398/1‑90 (MQ=255)
                                                                        gtcgtGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGAt                                                                               <  1:1012995/150‑1 (MQ=255)
                                                                            tgtTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCa                                                                           <  1:1697563/150‑1 (MQ=255)
                                                                              ttagAAGTTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGc                                                                         >  2:1901773/4‑150 (MQ=255)
                                                                              ttCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCTCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGc                                                                         >  1:1260077/1‑150 (MQ=255)
                                                                              ttCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGc                                                                         >  1:1317367/1‑150 (MQ=255)
                                                                                  aaGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTgg                                                                     <  2:2409039/150‑1 (MQ=255)
                                                                                        gggTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTgg                                                                     >  1:424569/1‑144 (MQ=255)
                                                                                        gggTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTgg                                                                     <  2:424569/144‑1 (MQ=255)
                                                                                             gtTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGgcggcg                                                          >  2:497493/1‑150 (MQ=255)
                                                                                             gtTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTCGCACTGgcggcg                                                          >  2:1513330/1‑150 (MQ=255)
                                                                                                 ggCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATa                                                      <  1:237030/150‑1 (MQ=255)
                                                                                                   cttgACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAg                                                                                              <  1:1861974/105‑1 (MQ=255)
                                                                                                   cttgACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAg                                                                                              >  2:1861974/4‑108 (MQ=255)
                                                                                                   cACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAgc                                                    <  1:1487210/150‑1 (MQ=255)
                                                                                                    aCGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAgcg                                                   <  1:1901773/150‑1 (MQ=255)
                                                                                                         ccGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAg                                                                                                                         >  1:2087490/1‑75 (MQ=255)
                                                                                                         ccGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAg                                                                                                                         <  2:2087490/75‑1 (MQ=255)
                                                                                                          cGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGa                                             >  2:1636430/1‑150 (MQ=255)
                                                                                                          cGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGa                                             <  2:16752/150‑1 (MQ=255)
                                                                                                          cGCATGACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGGGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGa                                             <  1:439199/150‑1 (MQ=255)
                                                                                                           gCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCAcgcg                                                                                                                               >  1:1868958/1‑67 (MQ=255)
                                                                                                           gCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCAcgcg                                                                                                                               <  2:1868958/67‑1 (MQ=255)
                                                                                                             aTTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCa                                          <  1:1636430/150‑1 (MQ=255)
                                                                                                              ttACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCAt                                         >  2:128276/1‑150 (MQ=255)
                                                                                                                   gCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATa                                                      >  1:1019423/1‑132 (MQ=255)
                                                                                                                   gCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATa                                                      <  2:1019423/132‑1 (MQ=255)
                                                                                                                    ccTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTggc                                                              >  2:1246567/1‑123 (MQ=255)
                                                                                                                    ccTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTggc                                                              <  1:1246567/123‑1 (MQ=255)
                                                                                                                    ccTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATc                                        <  2:1982166/145‑1 (MQ=255)
                                                                                                                    ccTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATc                                        >  1:1982166/1‑145 (MQ=255)
                                                                                                                     ctgctgGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAAc                                  <  1:497493/150‑1 (MQ=255)
                                                                                                                      tgctgGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAAcc                                 >  2:53922/1‑150 (MQ=255)
                                                                                                                        ctgGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATc                                        >  2:39361/1‑141 (MQ=255)
                                                                                                                        ctgGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATc                                        <  1:39361/141‑1 (MQ=255)
                                                                                                                          ggACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAg                                            <  1:1862784/135‑1 (MQ=255)
                                                                                                                          ggACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAg                                            >  2:1862784/1‑135 (MQ=255)
                                                                                                                             ccTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGcc                          >  2:1615742/1‑150 (MQ=255)
                                                                                                                               tGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCgg                        >  2:1794540/1‑150 (MQ=255)
                                                                                                                                gcgcCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGt                       <  1:1126443/150‑1 (MQ=255)
                                                                                                                                  gcCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGc                     <  2:1265695/150‑1 (MQ=255)
                                                                                                                                        ccGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTgc               >  1:282853/1‑150 (MQ=255)
                                                                                                                                          ggCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTgcgc             <  2:721327/150‑1 (MQ=255)
                                                                                                                                           gcagcaGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTACCGGTGCGTTTGcgcg            >  2:852566/1‑150 (MQ=255)
                                                                                                                                           gcagcaGCGGCCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTggcggc                                                           <  1:2234704/103‑1 (MQ=255)
                                                                                                                                           gcagcaGCGGCCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTggcggc                                                           >  2:2234704/1‑103 (MQ=255)
                                                                                                                                           gcagcaGCGGCCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTGcgcg            >  1:268277/1‑150 (MQ=255)
                                                                                                                                            cagcagCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTGCGCga           <  2:1260077/150‑1 (MQ=255)
                                                                                                                                              gcagcGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTGCGCgaga         <  1:274287/150‑1 (MQ=255)
                                                                                                                                                     gCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTGCGCGAGAGCTGccc  <  1:1794540/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGCCAGCGGCCCGATCGGCGAAACCCTGGCTCGCCTTGACCTGTCGCTGGCCCAGGGCGACGACGCCGGCGTCGTGTTCGAAGGTAGGGTGGTTGGGCACGACCCGCATTACGGCCTGCTGGACCTGCGCCTGCCCGGCAGCAGCGGGCCGCTGCTGCGCATCACCCACGCGGCGCAGGTGATGGGCAGCACCCTGCGCGTCAAGGTGCAGGCCCGGGATGTCAGCCTGGCACTGGCGGCGGATAGCGCATCGAGCATCCTCAACCGCCTGCCGGTGCGTGTGCGCGAGAGCTGCCC  >  NC_002947/4356169‑4356467

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