Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 497,714 A→C 5.9% D528A (GAC→GCC) ‡ PP_0409 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947497,7140AC 5.9% 593.1 / 17.4 179D528A (GAC→GCC) ‡PP_0409two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (90/78);  new base C (5/6);  total (95/84)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

CCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGCGGGCCGAGCTGG  >  NC_002947/497565‑497860
                                                                                                                                                     |                                                                                                                                                   
ccccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGa                                                                                                                                                     <  1:2245883/150‑1 (MQ=255)
 cccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                                    <  1:4689239/150‑1 (MQ=255)
 cccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                                    >  2:2263275/1‑150 (MQ=255)
 cccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                                    >  1:5030253/1‑150 (MQ=255)
 cccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                                    >  1:3736389/1‑150 (MQ=255)
 cccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                                    >  2:3448496/1‑150 (MQ=255)
 cccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                                    >  1:3334668/1‑150 (MQ=255)
  ccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACt                                                                                                                                                   >  2:264777/1‑150 (MQ=255)
   cTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTg                                                                                                                                                  <  1:4636341/150‑1 (MQ=255)
      tatCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTc                                                                                                                                               <  2:1405702/150‑1 (MQ=255)
        tCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCt                                                                                                                                            <  1:4497020/150‑1 (MQ=255)
         cAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTg                                                                                                                                           <  1:4972570/150‑1 (MQ=255)
         cAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTg                                                                                                                                           <  2:2910482/150‑1 (MQ=255)
         cAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTg                                                                                                                                           <  2:4555901/150‑1 (MQ=255)
         cAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGGACCAGTCCTGGTcgctg                                                                                                                                           <  1:2671881/150‑5 (MQ=255)
            tGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCa                                                                                                                                        <  1:4771237/150‑1 (MQ=255)
             gAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAg                                                                                                                                       <  2:1528422/150‑1 (MQ=255)
               gggTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACGGGGTGCAGGACCTGGCCCTGGCCAGGACCAGTCCTGGTCGCTGGccagg                                                                                                                                     <  1:1063923/150‑5 (MQ=255)
               gggTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACt                                                                                                                                                   <  2:3740848/137‑1 (MQ=255)
               gggTAGCGAGGCTGGGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAggg                                                                                                                                     <  2:2741580/150‑1 (MQ=255)
                ggTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAgggg                                                                                                                                    >  2:1987793/1‑150 (MQ=255)
                ggTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAgggg                                                                                                                                    <  1:707183/150‑1 (MQ=255)
                    gCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAGGACCAGTCCTGGTCGCTGGCCAGGgccag                                                                                                                                <  2:570068/150‑4 (MQ=37)
                     cGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAGGACCAGTCCTGGTCGCTGGCCAGGgccagg                                                                                                                               <  2:1930808/150‑5 (MQ=37)
                     cGAGGCTGTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCAt                                                                                                                               >  2:4614984/1‑150 (MQ=255)
                      gAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATc                                                                                                                              >  2:2041702/1‑150 (MQ=255)
                      gAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATc                                                                                                                              >  2:4064488/1‑150 (MQ=255)
                       aGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATcc                                                                                                                             <  2:271815/150‑1 (MQ=255)
                            ttGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCAcc                                                                                                                        <  2:5338989/150‑1 (MQ=255)
                            ttGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCAcc                                                                                                                        <  1:527420/150‑1 (MQ=255)
                              ggtggGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCg                                                                                                                      >  1:817066/1‑150 (MQ=255)
                              ggtggGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCg                                                                                                                      >  1:3374312/1‑150 (MQ=255)
                                 gggTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCa                                                                                                                   <  2:1323/150‑1 (MQ=255)
                                       tGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGa                                                                                                             >  1:187426/1‑150 (MQ=255)
                                          aCCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACa                                                                                                          >  1:2525979/1‑150 (MQ=255)
                                          aCCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACa                                                                                                          >  1:4308778/1‑150 (MQ=255)
                                          aCCAAGCGAGATCCAGAAACAGCAGATGTGACTCCGAAAACTTGTGCACCACTAGGACCGTGATGAATATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACAAGGACTCTTCCTGGCATGGCAGCATCCTCACACGGCAGGGAGAACa                                                                                                          >  2:3263030/1‑150 (MQ=255)
                                           cccAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCACACCCGGCAGGGCGAACAg                                                                                                         <  2:4979954/150‑1 (MQ=255)
                                            ccAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAgc                                                                                                        >  2:4833793/1‑150 (MQ=255)
                                              aGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAgcgc                                                                                                      <  1:2275804/150‑1 (MQ=255)
                                                cGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACCGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTg                                                                                                    <  2:902226/150‑1 (MQ=255)
                                                 gagaTCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTgg                                                                                                   >  2:2535168/1‑150 (MQ=255)
                                                 gagaTCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTgg                                                                                                   <  1:1662832/150‑1 (MQ=255)
                                                   gaTCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGc                                                                                                 <  1:456573/150‑1 (MQ=255)
                                                    aTCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGcc                                                                                                >  1:3200004/1‑150 (MQ=255)
                                                    aTCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGcc                                                                                                >  2:790033/1‑150 (MQ=255)
                                                     tCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGCCAGGACCAGTCCTGGTCGCTGGccaggg                                                                                                                                    >  1:2360488/1‑107 (MQ=37)
                                                      ccAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGa                                                                                              >  2:4965220/1‑150 (MQ=255)
                                                      ccAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGa                                                                                              <  1:2087661/150‑1 (MQ=255)
                                                      ccAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGa                                                                                              <  2:4302844/150‑1 (MQ=255)
                                                       cAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGat                                                                                             >  1:2860929/1‑150 (MQ=255)
                                                       cAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGat                                                                                             >  1:5044371/1‑150 (MQ=255)
                                                            cccGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCaa                                                                                        <  2:4925910/150‑1 (MQ=255)
                                                            cccGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCaa                                                                                        <  2:5324675/150‑1 (MQ=255)
                                                            cccGCACATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCaa                                                                                        >  2:162535/1‑150 (MQ=255)
                                                             ccGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAg                                                                                       <  1:1469392/150‑1 (MQ=255)
                                                                       gCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACt                                                                                                                                                   >  2:4976258/1‑81 (MQ=255)
                                                                       gCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACt                                                                                                                                                   <  1:4976258/81‑1 (MQ=255)
                                                                       gCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGccc                                                                             >  1:4137450/1‑150 (MQ=255)
                                                                        ccTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCg                                                                            <  1:3448496/150‑1 (MQ=255)
                                                                          tGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGc                                                                          <  1:264777/150‑1 (MQ=255)
                                                                          tGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGc                                                                          >  2:4520654/1‑150 (MQ=255)
                                                                              cAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCt                                                                      <  1:1436879/150‑1 (MQ=255)
                                                                               aaCCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTg                                                                     >  1:1358362/1‑150 (MQ=255)
                                                                                  cTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGc                                                                  >  1:2037514/1‑150 (MQ=255)
                                                                                  cTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCACATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGc                                                                  >  2:977964/1‑150 (MQ=255)
                                                                                   tGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGc                                                                  <  2:2982378/149‑1 (MQ=255)
                                                                                   tGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGc                                                                  >  1:2982378/1‑149 (MQ=255)
                                                                                   tGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCa                                                                 >  2:5096479/1‑150 (MQ=255)
                                                                                   tGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCa                                                                 >  2:5096464/1‑150 (MQ=255)
                                                                                       gCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACgg                                                             >  2:1498255/1‑150 (MQ=255)
                                                                                             cGAGGACCGTGCTGATTATCCCAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGt                                                       <  1:4074418/150‑1 (MQ=255)
                                                                                             cGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGt                                                       >  1:108437/1‑150 (MQ=255)
                                                                                                ggACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAgt                                                    <  2:108437/150‑1 (MQ=255)
                                                                                                    cGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTgg                                                <  1:676321/150‑1 (MQ=255)
                                                                                                    cGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTgg                                                <  1:3411464/150‑1 (MQ=255)
                                                                                                         tGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCTGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACgg                                           >  1:522287/1‑150 (MQ=255)
                                                                                                         tGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACgg                                           >  1:3798615/1‑150 (MQ=255)
                                                                                                          gATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGgt                                          >  1:2075216/1‑150 (MQ=255)
                                                                                                          gATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGgt                                          <  1:304199/150‑1 (MQ=255)
                                                                                                           aTTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGgtg                                         >  1:4421736/1‑150 (MQ=255)
                                                                                                           aTTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGgtg                                         <  2:521261/150‑1 (MQ=255)
                                                                                                             tATCACAGGGGGCAGGACCTGGCCCTGGCCAGCGACCAGGCCTGGTCCTGGCGGGGGCGCATCCTCACCCGGCCGGGCGACCCGCGCTGGGCCGCGAGCACGGCCCGCGCCCGGCATCTGGGCAGCGGCCGGGTCGGGGGGGCCGGTGTa                                       >  2:341032/1‑150 (MQ=255)
                                                                                                                cacaGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCGGGTAGTGTGGGACGGTGTAGTc                                    <  1:1498255/150‑1 (MQ=255)
                                                                                                                  caGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTg                                  <  1:4457501/150‑1 (MQ=255)
                                                                                                                  caGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTg                                  >  1:890133/1‑150 (MQ=255)
                                                                                                                   aGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTgg                                 <  2:3163809/150‑1 (MQ=255)
                                                                                                                   aGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTgg                                 <  2:3646241/150‑1 (MQ=255)
                                                                                                                             tcctggccctggcCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATcacac                       >  2:2250387/2‑150 (MQ=255)
                                                                                                                             aCCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTCGTCTGGGACATcacac                       >  2:3498389/1‑150 (MQ=255)
                                                                                                                             aCCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATcac                         <  2:2711875/148‑1 (MQ=255)
                                                                                                                                 ggccctggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAggg                   >  1:3566069/1‑150 (MQ=255)
                                                                                                                                 ggccctggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAggg                   >  1:2441584/1‑150 (MQ=255)
                                                                                                                                 ggccctggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAggg                   <  1:990730/150‑1 (MQ=255)
                                                                                                                                       ctgccaggACCAGTCCTGGTCGCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGc              >  1:2474387/8‑150 (MQ=255)
                                                                                                                                        tgccaggACCAGTCCTGGTCGCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGCg             >  2:4112478/7‑150 (MQ=255)
                                                                                                                                                ccAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGCGGGCCGAgct    <  1:4079826/150‑1 (MQ=255)
                                                                                                                                                 cAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGCGGGCGGAgctg   >  1:1155956/1‑150 (MQ=255)
                                                                                                                                                 cAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGCGGGCCGAgctg   <  2:4308778/150‑1 (MQ=255)
                                                                                                                                                  aGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGCGGGCCGAgctgg  >  1:5361170/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGCGGGCCGAGCTGG  >  NC_002947/497565‑497860

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