Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,344,121 A→C 5.4% C8G (TGC→GGC)  PP_1168 ← TRAP dicarboxylate transporter subunit DctQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,344,1210AC 5.4% 490.2 / 13.7 149C8G (TGC→GGC) PP_1168TRAP dicarboxylate transporter subunit DctQ
Reads supporting (aligned to +/- strand):  ref base A (80/61);  new base C (5/3);  total (85/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 152 total aligned reads displayed.

GCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGACACCGG  >  NC_002947/1343972‑1344270
                                                                                                                                                     |                                                                                                                                                     
gCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCa                                                                                                                                                       <  2:1238968/150‑1 (MQ=255)
  gggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCt                                                                                                                                                     <  1:84957/150‑1 (MQ=255)
  gggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACt                                                                                                                                                     >  1:3657624/1‑150 (MQ=255)
    gCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGc                                                                                                                                                   <  2:5304955/150‑1 (MQ=255)
     cGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGcc                                                                                                                                                  <  2:1259530/150‑1 (MQ=255)
     cGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGcc                                                                                                                                                  <  2:3038556/150‑1 (MQ=255)
      gAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCTCCGCACTGCCg                                                                                                                                                 >  2:3342689/1‑150 (MQ=255)
       aaGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGa                                                                                                                                                <  1:1359238/150‑1 (MQ=255)
        aGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAg                                                                                                                                               <  1:2768385/150‑1 (MQ=255)
          aTGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGCCgcttt                                                                                                                                             <  1:3053445/150‑5 (MQ=255)
            gCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCg                                                                                                                                           >  1:4115598/1‑150 (MQ=255)
             ccccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGa                                                                                                                                          >  1:1910832/1‑150 (MQ=255)
               ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                                        >  1:2939014/1‑150 (MQ=255)
                cACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTg                                                                                                                                       >  2:2884922/1‑150 (MQ=255)
                     gttgatCCGCGTCATGGCCTGGATCGACCGGCCCGCGATCCTGATGCAACTGGGGTACAGGGTTTCGTTCACGCTCAGGAAAAGCTATTTCATTGGCATCACCAAAGCGGCCGGGCGGGGGGCGCCGCACTGCCGAGGTCGATTGGGCaa                                                                                                                                  <  1:1720722/148‑1 (MQ=255)
                      atgatCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAg                                                                                                                                 >  1:3120674/1‑150 (MQ=255)
                      atgatCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAg                                                                                                                                 >  1:391335/1‑150 (MQ=255)
                         atCAGCGTCATGGCCAGGATCGACAGGACCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGctcgaaa                                                                                                                              >  2:3077539/1‑144 (MQ=255)
                          tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                                             >  1:1373211/1‑150 (MQ=255)
                          tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                                             >  2:378802/1‑150 (MQ=255)
                            cgCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGACGTTCACGCTCAGGAAAAGCGATTACATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                                           <  1:5347811/149‑1 (MQ=255)
                            aGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                                           >  2:3901748/1‑150 (MQ=255)
                            aGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                                           <  1:5347783/150‑1 (MQ=255)
                             gCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                                          <  2:1768956/150‑1 (MQ=255)
                              cGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTg                                                                                                                         >  1:2932254/1‑150 (MQ=255)
                                    ggCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATCGGCATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTg                                                                                                                   <  1:2224018/150‑1 (MQ=255)
                                    ggCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCACCGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTg                                                                                                                   <  1:4812677/150‑1 (MQ=255)
                                    ggCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTg                                                                                                                   >  2:2151508/1‑150 (MQ=255)
                                       cAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGAc                                                                                                                >  1:2696785/1‑150 (MQ=255)
                                       cAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGAc                                                                                                                >  1:956853/1‑150 (MQ=255)
                                       cAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGAc                                                                                                                >  1:956897/1‑150 (MQ=255)
                                            tCGACCGGCCCGCGATCCAGATGCAACTGCGGTACCGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAACGCGGCAGGGCGGGGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGt                                                                                                           <  1:5246264/150‑1 (MQ=255)
                                            tCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGt                                                                                                           <  2:1910832/150‑1 (MQ=255)
                                             cGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTc                                                                                                          <  2:4490594/150‑1 (MQ=255)
                                                  ggCCCGCGATCCAGAGGCAAGTGGGGTACCGCGTGTCGTTCACGCTCAGGAAAAGCGCTTTCATCGGCATCACCAACGCGGCCGGGCGGGGGGCGCCGCAGTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATAc                                                                                                     <  1:3653039/150‑1 (MQ=255)
                                                    cccGCGATCCTGATGCAACTGCGGTACTGCGTGTCGTTCACGCTCCGGAAAAGCTTTTTCATAGGCATCACCAACGCGGCTGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                                   <  1:2550738/150‑1 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                                   >  1:3775104/1‑150 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                                   <  1:3180740/150‑1 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGAGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                                   >  1:3570633/1‑150 (MQ=255)
                                                      cgcgATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCtt                                                                                                 >  1:3505431/1‑150 (MQ=255)
                                                      cgcgATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCtt                                                                                                 >  2:2604935/1‑150 (MQ=255)
                                                           tCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTtcatca                                                                                            <  1:4920550/150‑1 (MQ=255)
                                                            ccAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAg                                                                                           >  1:575937/1‑150 (MQ=255)
                                                            ccAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAg                                                                                           >  2:2833220/1‑150 (MQ=255)
                                                             cAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAgg                                                                                          >  1:3273441/1‑150 (MQ=255)
                                                              aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGt                                                                                         >  1:2833393/1‑150 (MQ=255)
                                                                aTGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCg                                                                                       >  1:1044508/1‑150 (MQ=255)
                                                                 tGCAACTGGGGTACTGGGTGTCGTTCACGCTCAGGAAAAGCTATTTCATAGGCATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCgg                                                                                      <  1:3370618/150‑1 (MQ=255)
                                                                      cTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTAc                                                                                 >  2:3511756/1‑150 (MQ=255)
                                                                         cgttACAGCTTGTCGTTCACGCTCAGTAAAAGCGATTTCATAGGAATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgc                                                                              <  1:481991/147‑1 (MQ=255)
                                                                         cGGTACCGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgc                                                                              <  1:1724067/150‑1 (MQ=255)
                                                                             accgggtgtCGTCCACGCTCAGGAAAAGCGATTCCATCGGCACCACCAACGCGGCTGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgcgc                                                                          <  1:3575/145‑1 (MQ=255)
                                                                             aCAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAGTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgcgc                                                                          >  2:3590769/1‑150 (MQ=255)
                                                                             aCAGCGTGTCGTTCACGCTCAGGAAAAGCGATGTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCGGGGACTGGACGGCGTCGATACGCTTCCTCAGGTCGGCGTCCTgcgcgc                                                                          >  2:3590790/1‑150 (MQ=255)
                                                                                   tgtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATAc                                                                                                     >  2:1804533/1‑117 (MQ=255)
                                                                                    gtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCTGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAttt                                                                   >  2:3240564/1‑150 (MQ=255)
                                                                                       gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTc                                                                <  2:260218/150‑1 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcg                                                               >  1:2888531/1‑150 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcg                                                               >  2:1101424/1‑150 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcg                                                               >  2:2394523/1‑150 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcg                                                               >  1:824926/1‑150 (MQ=255)
                                                                                         tCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgc                                                              <  1:3240564/150‑1 (MQ=255)
                                                                                               tCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCg                                                        >  2:1041869/1‑150 (MQ=255)
                                                                                               tCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCg                                                        >  1:5337290/1‑150 (MQ=255)
                                                                                               tCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCg                                                        >  2:1471260/1‑150 (MQ=255)
                                                                                               tCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAGTGGGCAAGGGCGGCGGGTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCg                                                        >  2:4317981/1‑150 (MQ=255)
                                                                                                 aGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGaa                                                      <  1:1882104/150‑1 (MQ=255)
                                                                                                   gAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTCCTGCGCGCCGTATTTTTCGCGTACCGAAGc                                                    >  1:5351127/1‑150 (MQ=255)
                                                                                                       aGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTg                                                <  1:3336737/150‑1 (MQ=255)
                                                                                                        gCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTCCTGCGCGCCGTATTTTTCGCGTCCCGAAGCCGTgg                                               >  1:1628647/1‑150 (MQ=255)
                                                                                                          gATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcg                                             >  2:2488713/1‑150 (MQ=255)
                                                                                                           aTTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgt                                            <  1:3511756/150‑1 (MQ=255)
                                                                                                           aTTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgt                                            >  2:1029606/1‑150 (MQ=255)
                                                                                                               cATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTCCTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTa                                        >  2:599632/1‑150 (MQ=255)
                                                                                                               cATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTa                                        >  2:599204/1‑150 (MQ=255)
                                                                                                                aTAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAg                                       >  2:4957356/1‑150 (MQ=255)
                                                                                                                      aTCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGTTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGc                                 >  1:3750937/1‑150 (MQ=255)
                                                                                                                       tCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGct                                >  1:4915402/1‑150 (MQ=255)
                                                                                                                         aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCATACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttc                              >  2:1630612/1‑150 (MQ=255)
                                                                                                                          ccAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtcg                                          >  2:5126518/1‑137 (MQ=255)
                                                                                                                           cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                            <  2:4392126/150‑1 (MQ=255)
                                                                                                                           cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                            >  2:2258529/1‑150 (MQ=255)
                                                                                                                             cctcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt                          >  1:1963660/4‑150 (MQ=255)
                                                                                                                             cctcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt                          >  2:1418214/4‑150 (MQ=255)
                                                                                                                             aaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt                          <  2:1253414/150‑1 (MQ=255)
                                                                                                                               gCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCg                        >  1:2774991/1‑150 (MQ=255)
                                                                                                                                  gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg                     >  1:4317393/1‑150 (MQ=255)
                                                                                                                                   cAGTGCGGCGCCCGCCGCCCTGCCGCGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACgg                    >  1:3792689/1‑150 (MQ=255)
                                                                                                                                   cAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgc                                                                            >  1:990625/1‑94 (MQ=255)
                                                                                                                                    agtgcggcgCCCGCCGCCCTGCCGCGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGt                   >  1:527506/4‑150 (MQ=255)
                                                                                                                                       gcggcgGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAAt                >  2:1620191/1‑150 (MQ=255)
                                                                                                                                        cggcggGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATg               >  1:986464/1‑150 (MQ=255)
                                                                                                                                        cggcggGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATg               <  1:4532756/150‑1 (MQ=255)
                                                                                                                                           ggggcgCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAAc            <  1:188431/149‑1 (MQ=255)
                                                                                                                                            gggcgCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGCCGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACt           <  1:4682185/150‑1 (MQ=255)
                                                                                                                                            gggcgCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACt           <  1:1620191/150‑1 (MQ=255)
                                                                                                                                                 ccgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGaca      <  2:986464/150‑1 (MQ=255)
                                                                                                                                                  cgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                            <  1:2223914/127‑1 (MQ=255)
                                                                                                                                                  cgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGacac     >  1:1065042/1‑150 (MQ=255)
                                                                                                                                                  cgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGGAGAACGGCTTCTTGTCGACGGTAATGAACTCGacac     >  2:3398095/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGACACCGG  >  NC_002947/1343972‑1344270

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: ‑