Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,100,001 T→G 6.6% intergenic (‑69/‑136) PP_1875 ← / → PP_1876 two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,100,0010TG 6.6% 259.0 / 7.2 91intergenic (‑69/‑136)PP_1875/PP_1876two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base G (3/3);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAG  >  NC_002947/2099856‑2100131
                                                                                                                                                 |                                                                                                                                  
ttGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCG                                                                                                                                <  2:216578‑M1/150‑1 (MQ=255)
ttGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATCTCCGTATATCGCCAAGCGGGGGCCGCTTGGCc                                                                                                                                >  2:895600‑M1/1‑149 (MQ=255)
 tgTATGACAGGCGATTGGTGAACTTGATATCAATACGTCATTAGCTGTTTGCCTGTTCCTTGCGCTGCGCCCGATCTGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGgcgc                                                                                                                               <  1:2377258/150‑1 (MQ=255)
 tGTATGACAGGCGATTGGTGAAGTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGC                                                                                                                               <  2:205728‑M1/150‑1 (MQ=255)
  gTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCC                                                                                                                              <  2:886427‑M1/150‑1 (MQ=255)
  gTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCC                                                                                                                              >  2:2297318‑M1/1‑150 (MQ=255)
    aTGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCC                                                                                                                            <  1:2246025‑M1/150‑1 (MQ=255)
         aGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATATCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCAAAGCGGGGGCCGCTTGGCGCCCCTCCGC                                                                                                                       >  1:756978‑M1/1‑150 (MQ=255)
           gCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGG                                                                                                                     >  2:177651‑M1/1‑150 (MQ=255)
             gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                   <  1:801362‑M1/150‑1 (MQ=255)
                   tGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCC                                                                                                             <  2:1999990‑M1/150‑1 (MQ=255)
                        ttGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTcc                                                                                                        >  2:820367/1‑150 (MQ=255)
                        ttGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCAAGCCCGcgga                                                                                                        >  1:1249931‑M1/1‑146 (MQ=255)
                        ttGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcc                                                                                                        >  2:1702974‑M1/1‑146 (MQ=255)
                          gATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCca                                                                                                      <  1:2325555/150‑1 (MQ=255)
                           atatCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCAAGCCCGcggaggg                                                                                                     <  2:1249931‑M1/150‑8 (MQ=255)
                           atatCCACGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCcac                                                                                                     >  1:765002/1‑150 (MQ=255)
                            tatCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGTGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCcaca                                                                                                    >  2:1229312/1‑150 (MQ=255)
                             atCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAg                                                                                                   <  1:790792/150‑1 (MQ=255)
                               ccATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacaggg                                                                                                 >  2:688554‑M1/1‑139 (MQ=255)
                                 aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTAc                                                                                               >  1:739844/1‑150 (MQ=255)
                                   gAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAg                                                                                             >  2:1119137/1‑150 (MQ=255)
                                      tCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGAGGGCTTGCCCATTCCTACAGGGACCGcgt                                                                                          >  1:272732/1‑148 (MQ=255)
                                       cccGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                         <  2:832232/150‑1 (MQ=255)
                                       cccGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                         <  2:739844/150‑1 (MQ=255)
                                        ccGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGC                                                                                                                               >  1:706304‑M1/1‑111 (MQ=255)
                                        ccGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGC                                                                                                                               <  2:706304‑M1/111‑1 (MQ=255)
                                         cGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                       <  2:157162/150‑1 (MQ=255)
                                                gTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGc                                                                                >  2:1858944/1‑150 (MQ=255)
                                                  ttCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggc                                                                              >  2:1834723‑M1/1‑120 (MQ=255)
                                                     ccGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCt                                                                                           >  2:1449601/1‑134 (MQ=255)
                                                     ccGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCt                                                                                           <  1:1449601/134‑1 (MQ=255)
                                                     ccGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                           <  1:1895851/150‑1 (MQ=255)
                                                     ccGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaa                                                                           >  1:865581‑M1/1‑117 (MQ=255)
                                                      cGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                       >  1:16322/1‑137 (MQ=255)
                                                      cGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                       <  2:16322/137‑1 (MQ=255)
                                                         tCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGTGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCaa                                                                       <  1:1229312/150‑1 (MQ=255)
                                                           cTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGc                                                                     <  1:2161876/150‑1 (MQ=255)
                                                           cTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGc                                                                     >  1:778923/1‑150 (MQ=255)
                                                            ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                         <  1:2165622/129‑1 (MQ=255)
                                                            ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                         >  2:2165622/1‑129 (MQ=255)
                                                            ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                       >  2:176827/1‑131 (MQ=255)
                                                            ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                       <  1:176827/131‑1 (MQ=255)
                                                            ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGcc                                                                    >  2:286764/1‑150 (MQ=255)
                                                            ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGcc                                                                    <  2:765002/150‑1 (MQ=255)
                                                              gcgctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCt                                                                  <  1:2297318/150‑1 (MQ=255)
                                                                gctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgt                                                                <  1:820367/150‑1 (MQ=255)
                                                                 ctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccccagggcttgcgttcagcctggggctaagctgtacctgtg                                                               <  2:1432489‑M1/150‑46 (MQ=255)
                                                                    cgcAAGCATAGCCGAATAGCGCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtggga                                                            >  1:1208768‑M1/1‑102 (MQ=255)
                                                                      cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaac                                                          <  1:688554‑M1/150‑51 (MQ=255)
                                                                      cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctatacctgtgggaac                                                          <  2:1984535‑M1/150‑51 (MQ=255)
                                                                      cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGc                                                          >  2:1613071‑M2/1‑149 (MQ=255)
                                                                      cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGG                                                                                                                    >  2:1325375‑M1/1‑92 (MQ=255)
                                                                      cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGG                                                                                                                    <  1:1325375‑M1/92‑1 (MQ=255)
                                                                       aaGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacg                                                         <  1:1834723‑M1/150‑52 (MQ=255)
                                                                       aaGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCAAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcg                                                         >  2:756978‑M2/1‑148 (MQ=255)
                                                                        aGCATAGCCGAATAGCGCCGGGTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        <  2:1208768‑M1/150‑53 (MQ=255)
                                                                        aGCATAGCCGAATAGCCCCGGCTGCTTCTTGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        >  1:1639627‑M1/1‑98 (MQ=255)
                                                                        aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                        >  2:778923‑M2/1‑147 (MQ=255)
                                                                        aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        <  1:1702974‑M1/150‑53 (MQ=255)
                                                                        aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        <  2:865581‑M1/150‑53 (MQ=255)
                                                                        aGCATAGACGAATAGCCCCGGCTGCTGGTGGATCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        <  2:1442348‑M1/150‑53 (MQ=255)
                                                                         gCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcggg                                                       >  1:177651‑M2/1‑146 (MQ=255)
                                                                          cATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        >  1:693839‑M1/1‑96 (MQ=255)
                                                                          cATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacatgtgggaacgg                                                        <  2:693839‑M1/148‑53 (MQ=255)
                                                                            tAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCaa                                                                       >  1:583149/1‑131 (MQ=255)
                                                                            tAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCaa                                                                       <  2:583149/131‑1 (MQ=255)
                                                                                     aGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcgga                                           >  1:216992‑M2/1‑134 (MQ=255)
                                                                                     aGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcgga                                           >  1:1858944‑M2/1‑134 (MQ=255)
                                                                                      gCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggag                                          >  1:1119137‑M2/1‑133 (MQ=255)
                                                                                       ccccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGGCGCGTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgggcaagcccgcggagg                                         <  2:1639627‑M1/150‑68 (MQ=255)
                                                                                         ccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgg                                                              <  1:815073‑M1/127‑47 (MQ=255)
                                                                                         ccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgg                                                              >  2:815073‑M1/1‑81 (MQ=255)
                                                                                          cGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggaggggc                                      >  1:286764‑M2/1‑129 (MQ=255)
                                                                                                 tggtggGTCTGTCTGTTCCTTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAAcgggcaagcccgcggaggg                                        >  2:2229309‑M2/1‑122 (MQ=255)
                                                                                                 tggtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAAcgggcaagcccgcggaggg                                        <  1:2229309‑M2/141‑20 (MQ=255)
                                                                                                   gtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaa                                                    >  2:1432518‑M2/1‑120 (MQ=255)
                                                                                                   gtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaa                                                    <  1:1432518‑M2/127‑8 (MQ=255)
                                                                                                   gtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAg                                >  2:487738/1‑147 (MQ=255)
                                                                                                   gtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAg                                <  1:487738/147‑1 (MQ=255)
                                                                                                     gggTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTTTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGcc                           <  1:895600/150‑1 (MQ=255)
                                                                                                                    tttCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        >  1:1829104‑M1/1‑54 (MQ=255)
                                                                                                                    tttCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                        <  2:1829104‑M1/106‑53 (MQ=255)
                                                                                                                        gatctggcGCCAAGCGGCCCCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGt        <  1:1317444/143‑1 (MQ=255)
                                                                                                                              gcGCCAAGCGGCCCCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCAAAAAAAGCGACATTGTTAgatt  >  2:1317444/2‑147 (MQ=255)
                                                                                                                                gCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgg                                                              <  2:1972956‑M1/88‑47 (MQ=255)
                                                                                                                                gCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgg                                                              >  1:1972956‑M1/1‑42 (MQ=255)
                                                                                                                                 ccAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGCACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaa                                                    <  2:796584‑M2/97‑8 (MQ=255)
                                                                                                                                 ccAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGCACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaa                                                    >  1:796584‑M2/1‑90 (MQ=255)
                                                                                                                                              gCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTg               <  2:2164122/121‑1 (MQ=255)
                                                                                                                                              gCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTg               >  1:2164122/1‑121 (MQ=255)
                                                                                                                                                 |                                                                                                                                  
TTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAG  >  NC_002947/2099856‑2100131

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: