Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,661,873 G→A 7.7% intergenic (‑43/+32) PP_3223 ← / ← PP_3224 ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,661,8730GA 7.7% 500.9 / 22.2 169intergenic (‑43/+32)PP_3223/PP_3224ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (82/73);  new base A (6/7);  total (88/81)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 170 total aligned reads displayed.

AAGGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3661724‑3662020
                                                                                                                                                      |                                                                                                                                                   
     gATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGc                                                                                                                                                >  1:5340212/1‑150 (MQ=255)
      aTGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGcc                                                                                                                                               <  1:5287691/150‑1 (MQ=255)
      aTGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGcc                                                                                                                                               >  1:54865/1‑150 (MQ=255)
       tGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCg                                                                                                                                              <  1:3257844/150‑1 (MQ=255)
        gTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCgg                                                                                                                                             <  2:4021539/150‑1 (MQ=255)
        gTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCgg                                                                                                                                             <  2:1568329/150‑1 (MQ=255)
           gAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAAc                                                                                                                                          >  2:467910/1‑150 (MQ=255)
           gAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCCGGTCGCCGGAAc                                                                                                                                          >  1:4696321/1‑150 (MQ=255)
           gAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCGCCTGTTCCGGCGAACTGGCCGGAAc                                                                                                                                          <  1:3721596/150‑1 (MQ=255)
           gAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCGCCTGTTCCGGCGAACTGGCCGGAAc                                                                                                                                          <  1:3721913/150‑1 (MQ=255)
             gggTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAg                                                                                                                                        >  1:491052/1‑150 (MQ=255)
                        gTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGc                                                                                                                             >  2:2239287/1‑150 (MQ=255)
                          tGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCt                                                                                                                           >  1:4993412/1‑150 (MQ=255)
                            aGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTg                                                                                                                         <  1:3946473/150‑1 (MQ=255)
                            aGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGGGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTg                                                                                                                         <  2:3189772/150‑1 (MQ=255)
                              gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  2:4355601/150‑1 (MQ=255)
                              gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  2:929268/150‑1 (MQ=255)
                                gttgtCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGagt                                                                                                                     >  2:3536194/1‑150 (MQ=255)
                                gttgtCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGagg                                                                                                                     >  2:3535889/1‑149 (MQ=255)
                                 tggtCGGCCGCGGAGCACGGGGGGCGGCCCAGGGCATGGGTGCGTTCGGGTCGCTGTGGGTGATTTGGGGGCATGGTTTTGATCCTTAATTGATGGGGGCTCCCATTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGagtca                                                                                                                   <  2:4155077/148‑1 (MQ=255)
                                       gccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCt                                                                                                              <  2:5382094/150‑1 (MQ=255)
                                       gccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCt                                                                                                              >  2:1337303/1‑150 (MQ=255)
                                        ccggcGAGCACGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTg                                                                                                             <  2:4973782/150‑1 (MQ=255)
                                        ccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAATGGTGCCTTGGAGTCAGTCCTg                                                                                                             <  2:3160201/150‑1 (MQ=255)
                                         cggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            <  1:1537978/150‑1 (MQ=255)
                                         cggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            >  1:4938582/1‑150 (MQ=255)
                                         cggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            >  1:4938514/1‑150 (MQ=255)
                                          ggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGt                                                                                                           >  2:813248/1‑150 (MQ=255)
                                          ggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGt                                                                                                           >  1:5246640/1‑150 (MQ=255)
                                          ggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGt                                                                                                           >  1:5246615/1‑150 (MQ=255)
                                             gAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAAt                                                                                                        <  1:4809520/150‑1 (MQ=255)
                                               gCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCg                                                                                                      <  1:3851170/150‑1 (MQ=255)
                                                 aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    >  1:356383/1‑150 (MQ=255)
                                                 aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCGCCTGTTCCGGCGAACTGGCCGGAACAGGc                                                                                                                                      >  2:330917/1‑116 (MQ=39)
                                                   gCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGc                                                                                                  <  1:3531019/150‑1 (MQ=255)
                                                   gCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGc                                                                                                  <  2:4095468/150‑1 (MQ=255)
                                                   gCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGc                                                                                                  >  2:4525857/1‑150 (MQ=255)
                                                     ggTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGTAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGccc                                                                                                >  1:1254492/1‑150 (MQ=255)
                                                      gTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCt                                                                                               >  1:441956/1‑150 (MQ=255)
                                                       tGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACa                                                                                                                                         <  2:197146/107‑1 (MQ=255)
                                                       tGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTc                                                                                              >  1:977069/1‑150 (MQ=255)
                                                           gCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGa                                                                                          >  2:1718638/1‑150 (MQ=255)
                                                            cccAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAAGGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTcaacat                                                                                         <  2:54865/150‑6 (MQ=255)
                                                             ccagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATc                                                                                        >  2:727161/1‑150 (MQ=255)
                                                             ccagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATc                                                                                        >  2:3751273/1‑150 (MQ=255)
                                                             ccagggccagggTGCCGATCGCGAGCAGATGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTGCGCCGGAACAGGCAAAAGGAGCCTTGCAGTCAGTCCTGGGAATCGAGGCCCTCGCGATc                                                                                        >  2:727274/1‑150 (MQ=255)
                                                              cagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATcc                                                                                       >  2:2563927/1‑150 (MQ=255)
                                                                gggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCaa                                                                                     <  1:422008/150‑1 (MQ=255)
                                                                 ggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGGCCTGGTAATCGAGGCCCTCGCGATCCAAc                                                                                    >  1:2363636/1‑150 (MQ=255)
                                                                  gccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACt                                                                                   >  1:1040914/1‑150 (MQ=255)
                                                                   ccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTg                                                                                  >  1:2712075/1‑150 (MQ=255)
                                                                     agggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGcc                                                                                <  1:1818213/150‑1 (MQ=255)
                                                                     agggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGcc                                                                                <  1:2239287/150‑1 (MQ=255)
                                                                     agggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGcc                                                                                <  2:2427602/150‑1 (MQ=255)
                                                                      gggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCg                                                                               <  1:813248/150‑1 (MQ=255)
                                                                        gTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGaa                                                                             >  2:30609/1‑150 (MQ=255)
                                                                           ccGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCa                                                                          <  2:4222164/150‑1 (MQ=255)
                                                                            cGATCGCGCGCAGTTGGTGCAGTGGGGGCATGGGAGTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGGCCTGGTAAGCGCGGCCCTCGCGATCCAACTGCCGACTCAg                                                                         >  2:4221439/1‑150 (MQ=255)
                                                                            cGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAg                                                                         <  1:2403369/150‑1 (MQ=255)
                                                                            cGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAg                                                                         <  2:207340/150‑1 (MQ=255)
                                                                                   gAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATcc                                                                  <  2:1117201/150‑1 (MQ=255)
                                                                                   gAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATcc                                                                  >  2:517028/1‑150 (MQ=255)
                                                                                        aGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAAt                                                              >  1:3445861/1‑150 (MQ=255)
                                                                                         gTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGTCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATg                                                             >  2:3763710/1‑150 (MQ=255)
                                                                                          ttGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAgg                                                                                                   >  2:4282216/1‑111 (MQ=255)
                                                                                          ttGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAgg                                                                                                   <  1:4282216/111‑1 (MQ=255)
                                                                                                 aGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGa                                                     <  2:1207937/150‑1 (MQ=255)
                                                                                                     ggggCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCtt                                                                                                                          >  1:997212/1‑77 (MQ=255)
                                                                                                     ggggCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAAt                                                                            >  2:360563/1‑123 (MQ=255)
                                                                                                           tGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTcccc                                           >  1:3781400/1‑150 (MQ=255)
                                                                                                                    tCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGtt                                                        <  1:1466741/128‑1 (MQ=255)
                                                                                                                    tCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGtt                                                        >  2:1466741/1‑128 (MQ=255)
                                                                                                                    tCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATc                                  <  2:1396963/150‑1 (MQ=255)
                                                                                                                     ccTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGGGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCg                                 <  2:3961024/150‑1 (MQ=255)
                                                                                                                     ccTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCg                                 >  2:2366312/1‑150 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  1:1839291/1‑150 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  1:2562248/1‑150 (MQ=255)
                                                                                                                                tgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCtt                      >  1:3349461/1‑150 (MQ=255)
                                                                                                                                   ggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGa                   <  2:4890609/150‑1 (MQ=255)
                                                                                                                                    gCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGAc                  <  1:2810621/150‑1 (MQ=255)
                                                                                                                                      tctcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCaa                                                               <  2:3862424/103‑1 (MQ=255)
                                                                                                                                      tctcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCaa                                                               >  1:3862424/1‑103 (MQ=255)
                                                                                                                                      tctcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTg                <  2:2788730/150‑1 (MQ=255)
                                                                                                                                         ccTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTc             >  1:1000443/1‑150 (MQ=255)
                                                                                                                                          cTGTTCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCa            >  1:4346916/1‑150 (MQ=255)
                                                                                                                                          cTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCa            >  1:3238897/1‑150 (MQ=255)
                                                                                                                                           tGTTCCGGCCAGTTCGCCGGACCAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCAt           >  2:2164379/1‑150 (MQ=255)
                                                                                                                                           tGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCAt           >  1:3485032/1‑150 (MQ=255)
                                                                                                                                           tGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCAt           >  2:2307046/1‑150 (MQ=255)
                                                                                                                                           tGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCAt           >  2:2307253/1‑150 (MQ=255)
                                                                                                                                             ttCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGc         >  2:2483997/1‑150 (MQ=255)
                                                                                                                                             ttCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGc         >  2:2725024/1‑150 (MQ=255)
                                                                                                                                                cGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGc                        >  2:3823528/1‑132 (MQ=255)
                                                                                                                                                 ggCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCt     <  2:491052/150‑1 (MQ=255)
                                                                                                                                                 ggCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCt     <  2:1839291/150‑1 (MQ=255)
                                                                                                                                                 ggCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCt     >  2:351157/1‑150 (MQ=255)
                                                                                                                                                  gCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCtt    <  1:1885587/150‑1 (MQ=255)
                                                                                                                                                  gCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCtt    <  1:4945972/150‑1 (MQ=255)
                                                                                                                                                   ccAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTg   >  1:2623336/1‑150 (MQ=255)
                                                                                                                                                    cAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGa  <  1:1525881/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AAGGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3661724‑3662020

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: