Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,668,721 A→C 7.0% intergenic (+30/‑36) acnA‑II → / → prpF aconitate hydratase/aconitate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,668,7210AC 7.0% 372.7 / 6.6 129intergenic (+30/‑36)acnA‑II/prpFaconitate hydratase/aconitate isomerase
Reads supporting (aligned to +/- strand):  ref base A (51/67);  new base C (9/2);  total (60/69)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

GCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCC  >  NC_002947/2668573‑2668869
                                                                                                                                                    |                                                                                                                                                    
gCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAg                                                                                                                                                     <  2:288228/150‑1 (MQ=255)
 cGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGACGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAgg                                                                                                                                                    <  1:1084251/150‑1 (MQ=255)
  gAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAggg                                                                                                                                                   >  2:4719056/1‑150 (MQ=255)
   tgTGCCTTGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGTGTGCTGCAGCGCTTTGCCCAGGACTTCCTCTAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGt                                                                                                                                                  <  1:2473571/149‑1 (MQ=255)
     tGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCCGGGGAt                                                                                                                                                >  1:3371608/1‑150 (MQ=255)
      gCCTGTAAGTACCGGTGACCTCCCGTCTGGACACCGCCGAGGAAGTGTCGATCAACGAAGCGGGCGGGGTGCTGCAGCGATTTGCCCAGCACTTCCTCTAAGCGACCGCATGAGCATCAATCGGGGGCCGAAAAGCGGCCCCAgactcat                                                                                                                                               <  2:1922545/150‑7 (MQ=255)
        cTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTc                                                                                                                                             >  2:3688549/1‑150 (MQ=255)
         tGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCa                                                                                                                                            >  2:161744/1‑150 (MQ=255)
         tGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCa                                                                                                                                            >  2:161742/1‑150 (MQ=255)
           gAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCAcc                                                                                                                                          >  2:416416/1‑150 (MQ=255)
              gTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGa                                                                                                                                       <  2:2579166/150‑1 (MQ=255)
                  cGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCTATCTACGAAGCGGGCGGGTTGCTGCAGCTCTTTCCCCAGGACTTCCTCCAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTAt                                                                                                                                   <  1:1599996/150‑1 (MQ=255)
                  cGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCCGCCCCAGGGTATTTCAACGGACAAt                                                                                                                                   >  2:2433293/1‑150 (MQ=255)
                    gTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCa                                                                                                                                 <  1:161742/150‑1 (MQ=255)
                    gTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCa                                                                                                                                 <  2:2252368/150‑1 (MQ=255)
                    gTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCa                                                                                                                                 <  1:161744/150‑1 (MQ=255)
                      gACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCCCAAAGCGGCCCCAGGGTATTTAACCGGACAATCAgg                                                                                                                               >  1:4309109/1‑150 (MQ=255)
                         cTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCCAAGCGGCCCCCCGGGATTTCACCCGCCTATCAGGATg                                                                                                                            >  1:2241349/1‑150 (MQ=255)
                         cTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATg                                                                                                                            >  2:207156/1‑150 (MQ=255)
                          tGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGc                                                                                                                           <  2:735446/150‑1 (MQ=255)
                           gCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGcc                                                                                                                          <  2:4786466/150‑1 (MQ=255)
                           gCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGcc                                                                                                                          <  1:3534958/150‑1 (MQ=255)
                           gCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGcc                                                                                                                          <  1:3534959/150‑1 (MQ=255)
                                  ggACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCCAAGCGGCCCCAGGGTATTTAAACGGACAATAAGGATGCCCAACaaa                                                                                                                   >  1:460975/1‑150 (MQ=255)
                                   gACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGaaaa                                                                                                                  <  1:4347378/150‑1 (MQ=255)
                                    acacCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAAt                                                                                                                 <  1:1467268/150‑1 (MQ=255)
                                    acacCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAAt                                                                                                                 <  2:5364758/150‑1 (MQ=255)
                                    acacCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAAt                                                                                                                 <  1:282909/150‑1 (MQ=255)
                                      acCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGt                                                                                                               <  2:2426486/150‑1 (MQ=255)
                                       ccgccgAGGAAGTGTCGATCTACGAACCGGCCGCGGTGCTGCAGCGCACTCCCCACCACTTCCTCGAACCCACCGCCTGAGCATGAAACCGGGGCCGCAAAGCCCCCCCAGGGTATTTCACCGCACTATCACGATGCCCAAGAAAATgac                                                                                                              <  2:4605725/150‑3 (MQ=255)
                                        cgccgAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTca                                                                                                             <  2:4530857/150‑1 (MQ=255)
                                        cgccgAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTca                                                                                                             <  2:4318316/150‑1 (MQ=255)
                                          ccgAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTcaca                                                                                                           <  1:5127070/150‑1 (MQ=255)
                                                aaGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCCCAAAGCGGCCCCACGGTATTTCACCGGACTAACAGGATGCCCAAGAAAAAGTCAAATGTAcc                                                                                                     >  1:1282070/1‑150 (MQ=255)
                                                  gtgtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTAcccc                                                                                                   <  2:5302274/150‑1 (MQ=255)
                                                   tgtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGGCGCCAAGCGGCCCCCGGGGATTTGAACCGCCTATCAGGATGCCCAAGACAATGTCACCTGTGccccc                                                                                                  >  1:215418/1‑149 (MQ=255)
                                                    gtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCACGGTAATTAACCCGACAATCAGGATGCCCAAGAAAAAGTAACAAGTACCCCAg                                                                                                 >  2:492761/1‑150 (MQ=255)
                                                       gATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATc                                                                                              >  2:815841/1‑150 (MQ=255)
                                                          cTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCaaa                                                                                           <  2:2116031/150‑1 (MQ=255)
                                                             cGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATc                                                                                        <  2:3424377/150‑1 (MQ=255)
                                                              gAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGTCTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATAAAAAtac                                                                                       >  1:4054113/1‑148 (MQ=255)
                                                                aGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATcccc                                                                                     <  1:4116690/150‑1 (MQ=255)
                                                                  cgggcgggGTGCTGCAGCGCTTTGCCCTGGACTTCCTCGAAGCGACCGCCTGAGCTTGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGc                                                                                   <  1:464414/150‑1 (MQ=255)
                                                                    ggcgggGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCACGGTATTTCACCGGACTATCAGGATGCCCAAGAAAAAGTCACAAGAAACCCAGAACAAAAAACCCGccc                                                                                 >  2:4211728/1‑149 (MQ=255)
                                                                     gcgggGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTc                                                                                                              >  2:4587965/1‑120 (MQ=255)
                                                                          gTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACa                                                                           <  1:87747/150‑1 (MQ=255)
                                                                            gctgcAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCCAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATc                                                                         >  2:2813048/1‑150 (MQ=255)
                                                                             ctgcAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCAACCGACTATCAGGATGCCCAAGAAAATTTCACATGTACCCCAGATCAAAATCCCCCCCACCTACATcc                                                                        >  2:147933/1‑150 (MQ=255)
                                                                               gcAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCCAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCAACCTCcgc                                                                      >  2:4440516/1‑150 (MQ=255)
                                                                                cAGCGCTTTGCCCTGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGGCGGAAAGCGGGCCCAGGGTATTTGAGCGGTCTATGAGGATTCGCAAGTAAATGTGACATTTTAGCCAGTTCAAAATAGCCGCCACCTACATAcgcg                                                                     >  2:3874347/1‑150 (MQ=255)
                                                                                cAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCcgcg                                                                     <  2:2164757/150‑1 (MQ=255)
                                                                                    gcTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGcggcgg                                                                 <  1:4429028/150‑1 (MQ=255)
                                                                                    gcTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGcggcgg                                                                 <  1:3289049/150‑1 (MQ=255)
                                                                                       ttGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCAc                                                              <  2:1329536/150‑1 (MQ=255)
                                                                                          cccAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGc                                                                                   >  1:4996029/1‑126 (MQ=255)
                                                                                          cccAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGc                                                                                   <  2:4996029/126‑1 (MQ=255)
                                                                                           ccAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCCGGGTATTTCACCGGACTATCCCGCTGCCCAAGAAAAGGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGc                                                          >  2:1126774/1‑150 (MQ=255)
                                                                                            cAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCa                                                         >  2:3079477/1‑150 (MQ=255)
                                                                                                    ccTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCgtg                                                 <  2:768816/150‑1 (MQ=255)
                                                                                                               ccGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCt                                      <  1:2385584/150‑1 (MQ=255)
                                                                                                                   cTGAGCATGAATCGGGGGCCGCCAAGCGGCCCCAGGGTATTTCACCCGACTATCAGGATGCCCAAGAAAAAGTCACAAGTACCCCCGATCAAAATACCCCCCCCCAACAACCCCGGCGGCACCACCAAAGGCGTGTTATTACGCCAGCaa                                  >  1:2706128/1‑150 (MQ=255)
                                                                                                                      aGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGAc                               <  2:141847/150‑1 (MQ=255)
                                                                                                                      aGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGAc                               <  1:815841/150‑1 (MQ=255)
                                                                                                                          tGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGc                           <  1:4719056/150‑1 (MQ=255)
                                                                                                                           gAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGcc                          <  1:4156132/150‑1 (MQ=255)
                                                                                                                             ccctggggCCGCTTTGCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAATGCGTGTTCTTCCGCCTGCAAGACCTGc                           >  1:4824143/5‑146 (MQ=255)
                                                                                                                             ccctggggCCGCTTTGCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAATGCGTGTTCTTCCGCCTGCAAGACCTGc                           <  2:4824143/142‑1 (MQ=255)
                                                                                                                              cctggggCCGCTTTGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTca                                                                                                             <  1:842727/61‑1 (MQ=37)
                                                                                                                              cctggggCCGCTTTGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTca                                                                                                             >  2:842727/4‑64 (MQ=37)
                                                                                                                              cctggggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGa                       >  2:1141231/4‑150 (MQ=255)
                                                                                                                                 ggggCCGCTTTGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCa                    >  1:3583317/1‑150 (MQ=255)
                                                                                                                                 ggggCCGCAAAGCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGAGCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCa                    >  2:393917/1‑150 (MQ=255)
                                                                                                                                  gggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAg                   <  1:2921410/150‑1 (MQ=255)
                                                                                                                                  gggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAg                   >  2:5249013/1‑150 (MQ=255)
                                                                                                                                  gggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACAGGTCCCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCACAGGCGTGTTCTTCCGCCTGCACGACCTGCCCGAGCAg                   >  2:1196700/1‑150 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAgg                  >  2:2680754/1‑150 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAgg                  <  2:4942525/150‑1 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAgg                  <  2:4251662/150‑1 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAgg                  <  1:4580827/150‑1 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAgg                  <  2:4942494/150‑1 (MQ=255)
                                                                                                                                      cGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTcc                                          >  1:3760403/1‑123 (MQ=255)
                                                                                                                                      cGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGgcgc               >  1:1178340/1‑150 (MQ=255)
                                                                                                                                      cGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGgcgc               >  1:3281731/1‑150 (MQ=255)
                                                                                                                                       gCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCa              >  1:3760966/1‑150 (MQ=255)
                                                                                                                                       gCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCa              >  1:1697892/1‑150 (MQ=255)
                                                                                                                                         tttgCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCaaa            <  1:5130990/147‑1 (MQ=255)
                                                                                                                                         aaaGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCaaa            >  1:2307351/1‑150 (MQ=255)
                                                                                                                                         aaaGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCaaa            >  1:4545286/1‑150 (MQ=255)
                                                                                                                                            gCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTc         <  1:5249013/150‑1 (MQ=255)
                                                                                                                                             cGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTcc        <  2:3259786/150‑1 (MQ=255)
                                                                                                                                             cGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTcc        <  1:3688549/150‑1 (MQ=255)
                                                                                                                                              gggcccAGGGTATTTCTCCGGACTATCAGGATGCCCAAGAAAATGTCACATGTTCCCCAGATCAATATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAATTTCCa       <  1:4348005/147‑1 (MQ=255)
                                                                                                                                              ggCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCa       <  2:1346086/150‑1 (MQ=255)
                                                                                                                                              ggCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCa       <  1:1883513/150‑1 (MQ=255)
                                                                                                                                                ccccAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAgg     >  2:5441137/1‑150 (MQ=255)
                                                                                                                                                 cccAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAGTGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGt    >  1:3533971/1‑150 (MQ=255)
                                                                                                                                                 cccAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGt    >  2:559619/1‑150 (MQ=255)
                                                                                                                                                  ccAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGc                                    <  2:4185745/117‑1 (MQ=255)
                                                                                                                                                  ccAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGc                                    >  1:4185745/1‑117 (MQ=255)
                                                                                                                                                   cAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTcc  <  1:2433293/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCC  >  NC_002947/2668573‑2668869

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: