Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 353,105 T→G 6.1% intergenic (+69/+49) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,1050TG 6.1% 624.8 / 2.3 235intergenic (+69/+49)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base T (120/93);  new base G (8/12);  total (130/105)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTT  >  NC_002947/352956‑353254
                                                                                                                                                      |                                                                                                                                                     
   aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                                                                    <  2:8047921/150‑1 (MQ=255)
    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  1:5515544/150‑1 (MQ=255)
    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  1:2208234/150‑1 (MQ=255)
    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  1:2859159/150‑1 (MQ=255)
     cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTggg                                                                                                                                                  >  2:3391668/1‑150 (MQ=255)
      aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGa                                                                                                                                                 <  1:6439449/150‑1 (MQ=255)
        cTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGc                                                                                                                                               <  1:6250839/150‑1 (MQ=255)
         ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCg                                                                                                                                              >  2:8793207/1‑150 (MQ=255)
          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                                             <  2:6441256/150‑1 (MQ=255)
          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                                             <  1:5908294/150‑1 (MQ=255)
          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                                             <  2:6264179/150‑1 (MQ=255)
           ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCggg                                                                                                                                            <  1:6662542/150‑1 (MQ=255)
            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                                                           <  1:7288185/150‑1 (MQ=255)
               tttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTg                                                                                                                                        >  2:6113334/1‑150 (MQ=255)
               tttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTg                                                                                                                                        >  1:9063286/1‑150 (MQ=255)
               tttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTg                                                                                                                                        >  2:3708062/1‑150 (MQ=255)
                         acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCGGCTCCATTGACCGGGGCTGCTTGGCAGCCCTTCGCGGGCACGCCCGCTCCCCCCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                                                              <  2:8233853/150‑1 (MQ=39)
                           acaGTGCGGGTGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACGGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGACCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGc                                                                                                                            <  2:4797031/150‑1 (MQ=255)
                           acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGc                                                                                                                            >  2:7376031/1‑150 (MQ=255)
                               tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcg                                                                                                                        >  2:4136235/1‑150 (MQ=255)
                               tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcg                                                                                                                        >  2:8742085/1‑150 (MQ=255)
                                 ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGaa                                                                                                                      >  2:3777872/1‑150 (MQ=255)
                                      ggCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCCACAGGTTTGCGGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                                                                 >  1:4960106/1‑150 (MQ=255)
                                        ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCCGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCACTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAAAGCGGccc                                                                                                               <  2:6629571/150‑1 (MQ=255)
                                        ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCACCCACACGGTTGCGTAGTTATCTGTGGGTGCGGGCGTGCCCGCGAAGGGCCCCGCAGCGGccc                                                                                                               >  1:4130980/1‑150 (MQ=255)
                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCCCGAAGGGCCGCGAAGCGGcccc                                                                                                              >  1:695241/1‑150 (MQ=255)
                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCGCCCCCAGGGTTGCTGAGATAACTGTGGGAGCGGGCGAGCCCCCGAAAGGCCCCCCAGCCGcccc                                                                                                              >  1:5758855/1‑150 (MQ=37)
                                               cTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCTACCACAGGTTTTCTGAGATCTCTGTGGGGTCGGGCGTGCCCGCGAAGGGCCCCGAAGCGGCCCCAAATGc                                                                                                        >  2:8550713/1‑150 (MQ=38)
                                                  cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGACGCCCCCTCCCACACGGTTTGTGTTATTAATGTGGGAGCGGGGGTGGCCGCCAATGGGCGCGAAGCGGGCCCaaattatta                                                                                                     >  1:2348683/1‑143 (MQ=21)
                                                     ggCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCtcccac                                                                                                  >  2:6672051/1‑145 (MQ=255)
                                                     ggCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGAGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAATGATTAAAc                                                                                                  >  2:7926650/1‑150 (MQ=255)
                                                     ggCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCGCCCCCACGGGGGGGGTGATGTGTGGGGGGGCGGGCGGGCCCGCGCCGGGGCGCGGCGCGGCCCCCCTTGCGGCAAc                                                                                                  >  1:255772/1‑150 (MQ=255)
                                                      gCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAcc                                                                                                 <  2:4542140/150‑1 (MQ=255)
                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                                                               >  2:8548965/1‑150 (MQ=255)
                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGACCCCCCGGTGGCGGTGATCGCTGTGGGGGCGGGCGGGGCCGCCAAGGGCCGCGCAGCGGCCCCCATTTCGTAAACCgg                                                                                               >  1:4238342/1‑149 (MQ=19)
                                                                ggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTg                                                                                       <  2:1096964/150‑1 (MQ=255)
                                                                ggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATATCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTg                                                                                       <  2:3240727/150‑1 (MQ=255)
                                                                 gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCACCGGGTTGCGGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGACGCGGCCCCCAATGCTTCAACCCTGGTTCGTGc                                                                                      >  1:2138992/1‑150 (MQ=38)
                                                                   tcgtcgTTCGCTGATCCATTGACCGGGCCTCCTTTCCAGCCCTGCCCGGGCACCCCCCCGCCCCCAGGGTGGCTGAGAGCGCGGGGGGAGCGGGCGGCCCCCCCCAGGGCCGCCAACCGGCCCCCATGGCTACAACCGTGGTACTTGCCg                                                                                    >  2:8985540/1‑150 (MQ=255)
                                                                    cgtcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGGGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                                                                   <  1:586614/150‑1 (MQ=255)
                                                                       cgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                                                                <  1:6005029/150‑1 (MQ=255)
                                                                       cgtTCGCTGATCCATTGACCGGGGCGGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                                                                <  1:2060884/150‑1 (MQ=255)
                                                                        gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCACTGTGGGAGCGGGCCTGCCCGCCAAGGGCCGCGAAGCGGCCCCAATAGCTACAACCGTGGTTCTTGCCGagcag                                                                               >  2:5643619/1‑150 (MQ=38)
                                                                        gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcag                                                                               >  2:2060060/1‑150 (MQ=255)
                                                                         ttCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                                                                              >  2:2243173/1‑150 (MQ=255)
                                                                                   cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                    <  2:5666109/150‑1 (MQ=255)
                                                                                   cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCGGGGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCACTTGCGTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                    <  2:7776267/150‑1 (MQ=38)
                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAg                                                                   >  1:7492372/1‑150 (MQ=255)
                                                                                     ttGACGGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGTTTTGCTGAGATATCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGc                                                                  <  1:5718621/150‑1 (MQ=255)
                                                                                     ccgcgaagggCTGCAAAGCAGCCCTTCGCGGGCACGCCCCCGCCCACAGGTTTGCGTAGTAGTCTGTGGGGGCGGGCCTGCCCCCGAAGGGCCCCGAAGCGGCCCCAAATTCTTAAACCGGGGTGATTGCCCCGCCGCGCGTGGTAAAGc                                                                  >  1:117180/8‑150 (MQ=255)
                                                                                       gACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATGACTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCTTCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCtg                                                                >  1:6810193/1‑149 (MQ=38)
                                                                                        attggggCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                                             >  1:6398327/4‑72 (MQ=17)
                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                               <  1:9059562/150‑1 (MQ=255)
                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                               >  2:4350556/1‑150 (MQ=255)
                                                                                             ggCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATc                                                          >  2:6939983/1‑150 (MQ=255)
                                                                                              gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcc                                                         <  1:8920626/150‑1 (MQ=255)
                                                                                              gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcc                                                         <  2:7199678/150‑1 (MQ=255)
                                                                                              gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcc                                                         <  2:5807254/150‑1 (MQ=255)
                                                                                               ctgctTTCCCGCCCTTCGCGGTCACGCCCTCTCCAAACGGTTTGCTGAGATCGCTGTGGTCGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc                                                        <  1:324708/150‑1 (MQ=38)
                                                                                                tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc                                                       >  2:7586937/1‑150 (MQ=255)
                                                                                                tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc                                                       <  2:4811527/150‑1 (MQ=255)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  2:8218152/150‑1 (MQ=255)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  2:3605307/150‑1 (MQ=255)
                                                                                                  cgttGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  1:6762522/148‑1 (MQ=255)
                                                                                                       cAGCCCTTCGCGGGCACGCCCCCTCCCACAGGGTTGCGGGGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCCAATTCTTCAAACCTGGTTCTTTCCGCGCAGCGCGTGGGAAAGCTCGGTATCCCCCAGTAta                                                >  1:2885643/1‑149 (MQ=25)
                                                                                                            ctgcgcgGGCACGCCCGCCCCCACCGGTTTGCTGAGAGCTCGGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccac                                           <  1:4682731/147‑1 (MQ=39)
                                                                                                                gcgGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCaccacc                                       >  1:4514372/1‑150 (MQ=255)
                                                                                                                gcgGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCaccacc                                       <  1:7376030/150‑1 (MQ=255)
                                                                                                                  gggCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCt                                     >  2:1415417/1‑150 (MQ=255)
                                                                                                                   ggCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTg                                    <  1:3707879/150‑1 (MQ=255)
                                                                                                                     cACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTc                                  <  2:8096565/150‑1 (MQ=255)
                                                                                                                     cACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTc                                  >  2:8523814/1‑150 (MQ=255)
                                                                                                                      aCGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                 <  2:2138992/150‑1 (MQ=255)
                                                                                                                           cGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATc                            >  2:2970120/1‑150 (MQ=255)
                                                                                                                           cGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATc                            >  2:5664322/1‑150 (MQ=255)
                                                                                                                            gCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                           >  1:866925/1‑150 (MQ=255)
                                                                                                                              tCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTg                         <  1:707281/150‑1 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       >  2:4713932/1‑150 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       >  1:3109875/1‑150 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       >  2:7479199/1‑150 (MQ=255)
                                                                                                                                 cacaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTCTCCTgcag                      >  2:5042480/1‑150 (MQ=255)
                                                                                                                                       gTTGGCTGAGATCCCGGTGCGAGCCGGCGCGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCCTCAACCGCGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCa                 <  1:1022820/150‑1 (MQ=255)
                                                                                                                                          tGCTGAGATCTCTGTGGGATAGGCCGAGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTTCTTCAACCGTGGAACTTTACGACCAGCGCGTGGTAAAGCTCTGTATCCCCGAGTATCCCAACCACCCTGTCGATATCCTGCAGCACAAGCt              >  2:2189493/1‑150 (MQ=255)
                                                                                                                                           gCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt             <  1:7566333/150‑1 (MQ=255)
                                                                                                                                            cTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttg            >  2:2657750/1‑150 (MQ=255)
                                                                                                                                            cTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttg            >  1:8282567/1‑150 (MQ=255)
                                                                                                                                             tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  2:5300165/1‑150 (MQ=255)
                                                                                                                                             tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  2:3856213/1‑150 (MQ=255)
                                                                                                                                             tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  2:603859/1‑150 (MQ=255)
                                                                                                                                               agaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgttg         >  2:8589499/1‑150 (MQ=255)
                                                                                                                                                gaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGc        >  1:7811471/1‑150 (MQ=255)
                                                                                                                                                 aTCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGcc       <  1:5116821/150‑1 (MQ=255)
                                                                                                                                                  tctcTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCg      <  1:8742085/150‑1 (MQ=255)
                                                                                                                                                  tctcTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCg      >  1:7636223/1‑150 (MQ=255)
                                                                                                                                                  tctcTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCg      >  2:6075909/1‑150 (MQ=255)
                                                                                                                                                    tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt    >  1:8167199/1‑150 (MQ=255)
                                                                                                                                                     ctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCACCCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt   <  2:2797227/150‑1 (MQ=255)
                                                                                                                                                     ctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt   >  2:1014018/1‑150 (MQ=255)
                                                                                                                                                     ctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt   >  1:2260392/1‑150 (MQ=255)
                                                                                                                                                      tgtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGttt  <  1:3777872/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTT  >  NC_002947/352956‑353254

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: