Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 798,666 C→G 6.5% intergenic (+18/+9) fklB‑I → / ← PP_0685 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947798,6660CG 6.5% 892.4 / 40.5 275intergenic (+18/+9)fklB‑I/PP_0685FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (130/126);  new base G (7/11);  total (138/137)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 278 total aligned reads displayed.

CTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTTCGTC  >  NC_002947/798520‑798815
                                                                                                                                                  |                                                                                                                                                     
cTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGc                                                                                                                                                    >  1:6588331/1‑150 (MQ=255)
 tgatgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGcc                                                                                                                                                   <  2:3633873/150‑1 (MQ=255)
  gatgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCg                                                                                                                                                  >  1:1794969/1‑150 (MQ=255)
  gatgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCg                                                                                                                                                  <  1:1828503/150‑1 (MQ=255)
  gatgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCg                                                                                                                                                  >  1:3437217/1‑150 (MQ=255)
   atgaACGCTGGCAGCAACTGGCGCCTGTACGTGCCAAGCGAGCTGGCTTCCGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCgg                                                                                                                                                 <  2:6984042/150‑1 (MQ=255)
   atgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCgg                                                                                                                                                 >  1:1271526/1‑150 (MQ=255)
    tgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGc                                                                                                                                                <  1:6920682/150‑1 (MQ=255)
    tgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGc                                                                                                                                                <  1:2418083/150‑1 (MQ=255)
    tgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGc                                                                                                                                                <  1:113962/150‑1 (MQ=255)
     gaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGcc                                                                                                                                               >  1:7103294/1‑150 (MQ=255)
      aaCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCt                                                                                                                                              <  1:5047457/150‑1 (MQ=255)
       aCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTc                                                                                                                                             <  1:6660849/150‑1 (MQ=255)
       aCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTc                                                                                                                                             <  1:2914368/150‑1 (MQ=255)
        cGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCa                                                                                                                                            >  1:4308760/1‑150 (MQ=255)
         gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTcca                                                                                                                                           <  1:4548407/150‑3 (MQ=255)
         gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCaa                                                                                                                                           <  1:2690034/150‑1 (MQ=255)
         gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                           <  1:7260609/150‑3 (MQ=255)
              cagcaATGGGCGCCTGTCCGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCAGCCCGCCGCACCGCGTTCTGGTGTTCGACGTCGCGCTGCTCGAGGTTCTGTATTGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                      <  2:3151603/150‑1 (MQ=255)
              cagcaACTGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                      <  2:5571611/150‑1 (MQ=255)
              cagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                      <  2:4664941/150‑1 (MQ=255)
              cagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCGGCCGGCCTCAATTCCa                                                                                                                                      <  2:6262091/150‑1 (MQ=255)
              cagcaAATGGCGCCTGTACGGGCCAAGCGTGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCTTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                      <  2:2207683/150‑1 (MQ=255)
              cagcaAATGGCGCCTGAACGTGCCAAGCGAGCTGGCTTACGGCCCCCAAGGCGTTGGCACAAACCCCCCCCACAGCGTTCTGGTTTTCGACGTCGACCTGCTCGACGTTCTGTAATGCCCTGCAGCCCGGCCCGGCCGCCCACAATTCCa                                                                                                                                      <  2:2778587/150‑1 (MQ=255)
               agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAg                                                                                                                                     >  2:2515649/1‑150 (MQ=255)
               agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAg                                                                                                                                     >  1:750928/1‑150 (MQ=255)
               agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAg                                                                                                                                     >  1:2737158/1‑150 (MQ=255)
               agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTCGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAg                                                                                                                                     <  2:5074453/150‑1 (MQ=255)
                gcaAATGGCGCGTTTACGTGCCAAGCTAGCTGGCTTACGGCGCCAATGGCTTTGGCAGAATCCCGCCGCACAGCGTTGTGGTGTTCTACGTCGAGCTGCTCGAGTTTCTGTATTGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGt                                                                                                                                    <  1:532430/150‑1 (MQ=255)
                   attGGCGCCTGTCCGTGCCAAGCGAGCTGGCTTACGGCGCCCACGGCGTGGGCAGCATCCCGCCGCACCGCGTTCTGGTGTTCGACGTCGCGCTGCTCGACGTTCTGTAAGGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTg                                                                                                                                 <  2:6302608/148‑1 (MQ=255)
                    aTGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGt                                                                                                                                >  2:3314135/1‑150 (MQ=255)
                           cTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTACCGGCCTCAATTCCAGTCTGTCCCCCCa                                                                                                                         <  2:1605141/150‑1 (MQ=255)
                              taggTGCCAAGCGAGCTGGCTTACGGCGCCCATGGCGTTGGCAGCATCCCGCCGCACCGGGTTGTGGTTTTCGACGTCGCGCTGCTCGACGTTCTTTAATGCCCTGGTGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAggg                                                                                                                      <  1:4782198/147‑1 (MQ=255)
                                 gTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGgcgc                                                                                                                   >  2:4065957/1‑150 (MQ=255)
                                   gCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCaa                                                                                                                 >  1:5957915/1‑150 (MQ=255)
                                    ccAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAAc                                                                                                                >  2:3512858/1‑150 (MQ=255)
                                    ccAAGCGAGCTGGCTAACGGCCCCCAAGCCGTTGTCAGAATCCCGCCCCACAGCGTTCTGGTGTTCTACGTCCAGATGCTCCACGTTATGTAATCCCCTGTAGGCCGCCCCTGCCGCCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAAc                                                                                                                <  2:653269/150‑1 (MQ=255)
                                      aaGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGc                                                                                                              >  1:4877731/1‑150 (MQ=255)
                                        gCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGccc                                                                                                            <  2:3805152/150‑1 (MQ=255)
                                                gCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAg                                                                                                    >  2:1915001/1‑150 (MQ=255)
                                                  ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGca                                                                                                  >  1:2573997/1‑150 (MQ=255)
                                                  ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGca                                                                                                  >  2:7397973/1‑150 (MQ=255)
                                                  ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGca                                                                                                  <  2:720180/150‑1 (MQ=255)
                                                     cGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaa                                                                                               <  2:3052515/150‑1 (MQ=255)
                                                      ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCGGGCCGGCCTcca                                                                                                                                           >  1:2131630/1‑103 (MQ=255)
                                                      ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaac                                                                                              >  1:3497467/1‑150 (MQ=255)
                                                      ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaac                                                                                              >  1:928567/1‑150 (MQ=255)
                                                      ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                           <  2:2131630/105‑3 (MQ=255)
                                                        cgcCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCACCGCCCGCGCATAGcagaacag                                                                                            >  1:2539377/1‑150 (MQ=255)
                                                         gcCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaacaga                                                                                           <  2:6123486/150‑1 (MQ=255)
                                                           ccAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaa                                                                                         <  2:6004748/150‑1 (MQ=255)
                                                                gCGTTGCCAGCATCCCGCCGAACCGGGTAGTGGTGTCCCACGTCGCGCGGCCCGACGTTCTGTACTCCCCTGGCGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGAt                                                                                    <  1:1340073/150‑1 (MQ=255)
                                                                 ggTTGGCAGCATCCCGCCGCACCGCGTTGTGGTGTTCGACGTCGCGGGGCTCGACGTTCTGTACGGCCCTGGAGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGAtt                                                                                   <  1:6959248/149‑1 (MQ=255)
                                                                  gtgggCAGCATCCCGCCGCACCGGTTTCGGGTTTTCGCCGTCGAGCTGCTCGCGGTTCTGTAAGGCCCTGGAGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTc                                                                                  <  1:4073403/147‑1 (MQ=255)
                                                                  gTTGGCAGCACCCCGCCGCACAGCGTTGTGGTGTTCGACGCCGAGCTGCTCGACGTTCTGTATGGCCCTGGAGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTc                                                                                  <  1:1529038/150‑1 (MQ=255)
                                                                   ttGGCAGCATCCCGCCGCACAGCGTTCGGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTcc                                                                                 <  2:6920670/150‑1 (MQ=255)
                                                                      gcagcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTa                                                                              <  2:5634219/150‑1 (MQ=255)
                                                                      gcagaaTCCCGCCGCACAGCTTTGTGGTGTTCGAGGTGGTGCTGCTCGAGGTTCTGTATTGCCCTGGAGGCGGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTa                                                                              <  2:862165/150‑1 (MQ=255)
                                                                        agcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTAcc                                                                            >  1:6167415/1‑150 (MQ=255)
                                                                        agcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTAcc                                                                            >  2:2112422/1‑150 (MQ=255)
                                                                               cgccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCt                                                                                                                                              <  1:6370953/77‑1 (MQ=255)
                                                                                gccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCtt                                                                    >  2:1827457/1‑150 (MQ=255)
                                                                                   gcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAg                                                                 <  2:7107308/150‑1 (MQ=255)
                                                                                    cacaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCCGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGc                                                                >  1:3177296/1‑150 (MQ=255)
                                                                                    cacaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGc                                                                >  2:6832231/1‑150 (MQ=255)
                                                                                     acaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCa                                                               >  2:434765/1‑150 (MQ=255)
                                                                                            tCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTATTGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCgg                                                        <  2:235967/150‑1 (MQ=255)
                                                                                                   ttCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGa                                                                  >  2:4906944/1‑133 (MQ=255)
                                                                                                   ttCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTc                                                 >  2:2968503/1‑150 (MQ=255)
                                                                                                   ttCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTc                                                 <  2:4185663/150‑1 (MQ=255)
                                                                                                      gACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                              <  1:7355444/150‑1 (MQ=255)
                                                                                                      gACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                              >  1:19633/1‑150 (MQ=255)
                                                                                                        cGTCGCGCTGCTCGACGTTCTGTAATGCCCTGGTGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAAtt                                            <  2:903329/150‑1 (MQ=255)
                                                                                                        cGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAAtt                                            <  2:3437217/150‑1 (MQ=255)
                                                                                                        cGTCGAGCGGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAAtt                                            <  1:3877684/150‑1 (MQ=255)
                                                                                                            gAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGc                                        >  2:1062589/1‑150 (MQ=255)
                                                                                                            gAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTCCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGc                                        >  2:5389058/1‑150 (MQ=255)
                                                                                                              gctgctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGccc                                      >  1:5588122/1‑150 (MQ=255)
                                                                                                               ctgctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCg                                     >  2:6735554/1‑150 (MQ=255)
                                                                                                                 gctCGACGTTCTGTAATGCCGTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCTTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCt                                   <  1:3672177/150‑1 (MQ=255)
                                                                                                                  ctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTc                                  >  2:6700733/1‑150 (MQ=255)
                                                                                                                   tCGACGTTCTGTAATGCCCGGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCa                                 <  1:4116644/150‑1 (MQ=255)
                                                                                                                     gACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATg                               >  1:3076527/1‑150 (MQ=255)
                                                                                                                     gACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATg                               >  2:5191382/1‑150 (MQ=255)
                                                                                                                         ttCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCa                           >  2:4724847/1‑150 (MQ=255)
                                                                                                                          tCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                          <  2:1271526/150‑1 (MQ=255)
                                                                                                                              tAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCa                      >  2:2908827/1‑150 (MQ=255)
                                                                                                                               aaTGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCtt                                                                    <  1:1143781/103‑1 (MQ=255)
                                                                                                                               aaTGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCtt                                                                    >  2:1143781/1‑103 (MQ=255)
                                                                                                                                 tggaatGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCAccc                   >  1:6823875/6‑150 (MQ=255)
                                                                                                                                  gCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCt                  >  2:3240562/1‑150 (MQ=255)
                                                                                                                                    ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTgg                >  2:5880194/1‑150 (MQ=255)
                                                                                                                                    ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTgg                >  2:1772856/1‑150 (MQ=255)
                                                                                                                                    ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTgg                >  2:2678295/1‑150 (MQ=255)
                                                                                                                                     cTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCAccc                                                            >  2:1420137/1‑105 (MQ=255)
                                                                                                                                      tGGAGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTc              <  1:334151/150‑1 (MQ=255)
                                                                                                                                         aGGCCGGCAGGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACg           <  2:774300/150‑1 (MQ=255)
                                                                                                                                         aGGCCGGCAGGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACAAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACg           >  1:5784172/1‑150 (MQ=255)
                                                                                                                                             cggccCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGt       >  1:4684898/1‑150 (MQ=255)
                                                                                                                                                  cTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTtcgtc  <  2:2573997/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTTCGTC  >  NC_002947/798520‑798815

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: