Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,129,417 G→T 14.0% V334L (GTA→TTA)  PP_4514 → alpha/beta hydrolase superfamily esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,129,4170GT14.0% 367.2 / 15.9 148V334L (GTA→TTA) PP_4514alpha/beta hydrolase superfamily esterase
Reads supporting (aligned to +/- strand):  ref base G (67/57);  new base T (3/21);  total (70/78)
Fisher's exact test for biased strand distribution p-value = 2.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

CGGCGTGACAGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGG  >  NC_002947/5129269‑5129566
                                                                                                                                                    |                                                                                                                                                     
cGGCGTGACAGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGt                                                                                                                                                      <  2:71938/150‑1 (MQ=255)
 ggCGTGACAGATTTCGGCCGTGCCCGCCACATGATCGACGCCGTTTACCGCGCCACTCGCCTGGTTGACCGGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGGGACGCAAGCCTGGCCGTGGCCCGCTGGCCGCGCCAGCGCACGCCGGTa                                                                                                                                                     <  1:1559036/150‑1 (MQ=255)
     tGACAGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTa                                                                                                                                                 >  2:3680739/1‑150 (MQ=255)
      gACAGATTTCGGCCGGGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTAc                                                                                                                                                <  1:4636357/150‑1 (MQ=255)
        cAGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCg                                                                                                                                              <  1:4357193/150‑1 (MQ=255)
         aGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGc                                                                                                                                             >  2:997153/1‑150 (MQ=255)
         aGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGc                                                                                                                                             >  2:4701666/1‑150 (MQ=255)
         aGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGc                                                                                                                                             >  2:3967928/1‑150 (MQ=255)
         aGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGc                                                                                                                                             >  2:289058/1‑150 (MQ=255)
           aTTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCaa                                                                                                                                           <  1:214116/150‑1 (MQ=255)
           aTTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCaa                                                                                                                                           >  1:3973640/1‑150 (MQ=255)
              tCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCTATcc                                                                                                                                        <  2:704651/150‑1 (MQ=255)
               cGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCg                                                                                                                                       <  2:1454545/150‑1 (MQ=255)
               cGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCg                                                                                                                                       >  2:5266037/1‑150 (MQ=255)
               cGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCg                                                                                                                                       <  1:673113/150‑1 (MQ=255)
               cGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCg                                                                                                                                       <  1:4685992/150‑1 (MQ=255)
                  ccGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGAt                                                                                                                                    <  2:1982513/150‑1 (MQ=255)
                        ccGCCACTTGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGaaaa                                                                                                                              <  1:1069955/150‑1 (MQ=255)
                        ccGCCACATGATCGACGCCGGTTACCGCGCCACTCGCGTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGaaaa                                                                                                                              <  1:2447944/150‑1 (MQ=255)
                        ccGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTCATTACCGCAATCCGGATCGaaaa                                                                                                                              >  2:1742457/1‑150 (MQ=255)
                        ccGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGaaaa                                                                                                                              >  2:2829352/1‑150 (MQ=255)
                        ccGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGaaaa                                                                                                                              >  1:1513468/1‑150 (MQ=255)
                         cGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAAc                                                                                                                             >  2:2406638/1‑150 (MQ=255)
                          gCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCGGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACg                                                                                                                            >  1:2780285/1‑150 (MQ=255)
                          gCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACg                                                                                                                            >  1:5263427/1‑150 (MQ=255)
                             acaTGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCTGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACt                                                                                                                         <  2:400924/150‑1 (MQ=255)
                                tGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCa                                                                                                                      >  2:3804288/1‑150 (MQ=255)
                                     gACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTc                                                                                                                 >  2:3186967/1‑150 (MQ=255)
                                           ggTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCgat                                                                                                           >  1:1381232/1‑150 (MQ=255)
                                           ggTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCgat                                                                                                           >  1:1381473/1‑150 (MQ=255)
                                                  gcgcCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGa                                                                                                    >  2:4452915/1‑150 (MQ=255)
                                                   cgcCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCGGCGCCAGCTCGTGGGCGACGCCCGCCTGGCCGGGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGTATCGAAAACGACTCCAAGGTCAGCGATGATGTGAt                                                                                                   <  2:5123553/150‑1 (MQ=255)
                                                    gcCACTCGCCTGCTTGACCGGCGCCGGGCCGCCCTGCGCCCGCTCGCGGGCGCCGCCCGCCTGGCGGGGGCCCGCGCGCCGCGCCAGCGCCCGCCGGTAATTACCGCACCCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATc                                                                                                  <  2:2602550/150‑1 (MQ=255)
                                                      cACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCg                                                                                                >  2:2180205/1‑150 (MQ=255)
                                                       aCTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGc                                                                                               >  2:2888300/1‑150 (MQ=255)
                                                        cTCGCCTGCTTGACCCGCGCCGGGCCGCCCTGCGCCAGCTCGAGGGCGCCGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAAGCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCt                                                                                              <  2:2250597/150‑1 (MQ=255)
                                                           gCCTGCTTGACCCGCGCGGGGCCGCCCGGCGCCTGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCAGGCGGGTAATTACCGCATTCCGGATCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGt                                                                                           <  2:2501478/150‑1 (MQ=255)
                                                           gCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGAGGTGATCCGCTCGt                                                                                           >  2:4896124/1‑150 (MQ=255)
                                                            ccTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCgatg                                                                                                          <  2:2631845/134‑1 (MQ=255)
                                                             cTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCCATCCGGATCGAAAACGACTCCAAGGTCCGCGAGGATGTGATCCGCTCGtat                                                                                         >  1:1949621/1‑150 (MQ=255)
                                                               gCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTatat                                                                                       >  1:974392/1‑150 (MQ=255)
                                                                   gACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATc                                                                                                                                   <  1:5003119/102‑1 (MQ=255)
                                                                   gACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATc                                                                                                                                   >  2:5003119/1‑102 (MQ=255)
                                                                             tGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCg                                                                         >  1:4764601/1‑150 (MQ=255)
                                                                             tGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCg                                                                         >  1:4764513/1‑150 (MQ=255)
                                                                             tGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCg                                                                         <  2:5272821/150‑1 (MQ=255)
                                                                               gccgccCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGAc                                                                       >  2:5230440/1‑150 (MQ=255)
                                                                               gccgccCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGAc                                                                       >  2:3275800/1‑150 (MQ=255)
                                                                                ccgccgGGCGCCCGGTGGGGGGCGCCGCCCGCGGGGCGGTGGCCCGCTCGCCGCGCCCGCGCAGGCGGTTAATTACCGCAAGCCGGTTCGACCACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACg                                                                      <  1:3938824/150‑1 (MQ=255)
                                                                                 cgccCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGc                                                                     >  2:1760990/1‑150 (MQ=255)
                                                                                 cgccCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGc                                                                     >  1:5272706/1‑150 (MQ=255)
                                                                                    ccTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCAcc                                                                  >  2:4474656/1‑150 (MQ=255)
                                                                                       gcgcCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCt                                                               <  2:2909994/150‑1 (MQ=255)
                                                                                        cgcCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTg                                                              >  1:1498189/1‑150 (MQ=255)
                                                                                         gcCAGCCCGTGGGCCACGCCAGCCTGGCGGGGGCCCGCTCGCCGCGCCAGCGCACGCCGTTATTTACCCAAATCCGGATCCAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTgg                                                             <  1:1157914/150‑1 (MQ=255)
                                                                                         cggcgtgcgctggcgcggcgagcgggccacggccaggctggcgtcgccAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTgg                                                             >  1:2028313‑M2/49‑150 (MQ=255)
                                                                                          ccAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGa                                                            >  2:5145278/1‑150 (MQ=255)
                                                                                             gtgcgctggcgcggcgagcgggccacggccaggctggcgtcgccAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACt                                                         >  2:748980‑M2/45‑150 (MQ=255)
                                                                                              cTCGTGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTg                                                        <  1:5379438/150‑1 (MQ=255)
                                                                                               tCGGGGGCGCCCCCCGCGTGGCCGGGGCCCGCGCGCGGCCCCCGCCCCCGCCTTTAATTACCCCAAGCCGCACCCAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTgg                                                       <  1:21137/150‑1 (MQ=255)
                                                                                               tCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTgg                                                       >  2:317863/1‑150 (MQ=255)
                                                                                                ggtgggCGACGCCTGCCGGGCGGTGCCCCGCGCGCCGCGCCCGCGCACGCCGTTAATTCCCGCAATCCGGAGCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGa                                                      <  1:5128816/147‑1 (MQ=255)
                                                                                                    GGCTACGCCAGCCTGGCCGTGGCCCGCCCGCCGCGCCTGCGCACGCCGTTAATTACCGAAACCCGGATCTAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGc                                                  >  2:1381473‑M1/1‑150 (MQ=255)
                                                                                                      CGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGTTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCt                                                >  2:4460583‑M1/1‑150 (MQ=255)
                                                                                                       GACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTg                                               <  1:2780449‑M1/150‑1 (MQ=255)
                                                                                                       GACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTg                                               <  2:2268518‑M1/150‑1 (MQ=255)
                                                                                                         cGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCa                                             <  2:2072063/150‑1 (MQ=255)
                                                                                                           ccAGCCTGGCGGTGGCCCGCTCGCCGCGCCCGCGCACGCCGGTAATTACCGCAATCCGGACCCAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGa                                           <  1:2543916/150‑1 (MQ=255)
                                                                                                            ccgcgtGGCGGGGGCCCGCCCGCCGCGCCAGCGCCCGCGGTTATTTACCGCAAGCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGAc                                          <  1:3735997/145‑1 (MQ=255)
                                                                                                            cAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGAGCGACGCACCGCTGGAACTGGACCGCCTGCAGAc                                          >  2:1404684/1‑150 (MQ=255)
                                                                                                               ccgggCCGTGGCCCGCTCGCCGCGCCCGCGCACGCCGGTAATTCCCGCAATCCGGAGCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGa                                       <  1:1321118/147‑1 (MQ=255)
                                                                                                                cTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGAc                                      >  1:3479235/1‑150 (MQ=255)
                                                                                                                 tGGCCGGGGCCCGCTCGCCGCGCCCGCGCACGCGGTTATTTACCGCAATCCGTAGCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACa                                     <  1:1245403/150‑1 (MQ=255)
                                                                                                                  ggCCGTGGCCCGGGCGCCGCGCCCGCGCCCGCCTTTAATTACCGAAATCCTTAGCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACAt                                    <  2:1949621/150‑1 (MQ=255)
                                                                                                                  ggCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGATCGCCTg                                               <  2:4816349/139‑1 (MQ=255)
                                                                                                                  ggCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGATCGCCTg                                               >  1:4816349/1‑139 (MQ=255)
                                                                                                                  ggCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACAt                                    <  2:4308645/150‑1 (MQ=255)
                                                                                                                    ggggggCCGGCTGGCCGCGCCTGCGCCGGCCGTTATTTACCGAAATCCGTATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATgg                                  <  1:3438776/146‑1 (MQ=255)
                                                                                                                    cgggggCCCGCGCGCCGCGCCCGCGCAGGCGGGTATTTCCCGAAAGCGGGAGCGACCACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATgg                                  <  1:1552709/146‑1 (MQ=255)
                                                                                                                    cggTGGCCCGCTCGCCGCGCCTGCGCACGCGGTTATTTACCGCAAGCCGGATCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATgg                                  <  1:5266037/148‑1 (MQ=255)
                                                                                                                    ccGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATgg                                  <  2:5263427/150‑1 (MQ=255)
                                                                                                                    ccGTGGCCCGCTCGCCGCGCCAGCCCACGCCGGTAATAACCGAAACCCGGACCGCAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATgg                                  <  1:289058/150‑1 (MQ=255)
                                                                                                                    ccGGGGCCCGCTCGCCGCGCCCGCGCACGCCGTTATTTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATgg                                  <  2:1498189/150‑1 (MQ=255)
                                                                                                                    ccGGGGCCCGCGCGCCCCGCCCGCGCACGCGGGTATTTACCGCAACCGGGATCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATgg                                  <  1:3263190/150‑1 (MQ=255)
                                                                                                                          cccGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAGGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGc                            >  2:821129/1‑150 (MQ=255)
                                                                                                                            cGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGACCTGGACCGCCTGCAGACCGACATGGGTACGCTg                          >  1:4062765/1‑150 (MQ=255)
                                                                                                                             gCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCgatga                                                                                                         <  2:3245034/70‑1 (MQ=255)
                                                                                                                             gCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCgatga                                                                                                         >  1:3245034/1‑70 (MQ=255)
                                                                                                                              cTCGCCGCGCCTGCGCACGCCGTTAATAACCGAAATCCGTATCCAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTa                        <  1:744343/150‑1 (MQ=255)
                                                                                                                                 gccgcGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATgg                     >  1:5439207/1‑150 (MQ=255)
                                                                                                                                 gccgcGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATgg                     >  1:3183138/1‑150 (MQ=255)
                                                                                                                                  ccgcGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGc                    >  1:5176196/1‑150 (MQ=255)
                                                                                                                                     cgccAGCGCACGCCGGTAATTACCGCACTCCGGATCTAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTg                 <  2:216165‑M2/146‑1 (MQ=255)
                                                                                                                                     cgccAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTg                 <  2:4364823‑M2/146‑1 (MQ=255)
                                                                                                                                       cccgcgcACGCCGTTAATTACCGAAATCCGTATCCAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGa               <  2:1196399/147‑1 (MQ=255)
                                                                                                                                       cccgcgcACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGa               <  2:994399/147‑1 (MQ=255)
                                                                                                                                             cACGCCGGTATTTCCCGCAACCCGGATCGAACACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACtt         <  2:5016523/150‑1 (MQ=255)
                                                                                                                                               cgcgGGTAATTACCGCAATCCGGATCCAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCg       <  2:1872746/150‑1 (MQ=255)
                                                                                                                                                   ggTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAg   >  1:3358978/1‑150 (MQ=255)
                                                                                                                                                    gTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAgg  <  1:4950634/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGCGTGACAGATTTCGGCCGTGCCCGCCACATGATCGACGCCGGTTACCGCGCCACTCGCCTGCTTGACCCGCGCCTGGCCGCCCTGCGCCAGCTCGAGGGCGACGCCAGCCTGGCCGTGGCCCGCTCGCCGCGCCAGCGCACGCCGGTAATTACCGCAATCCGGATCGAAAACGACTCCAAGGTCAGCGATGATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGG  >  NC_002947/5129269‑5129566

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