Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,366 Δ1 bp 6.9% intergenic (+136/+782) PP_1865 → / ← aroH transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,087,3660C. 6.9% 316.8 / 17.9 101intergenic (+136/+782)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base C (44/50);  new base . (3/4);  total (47/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 105 total aligned reads displayed.

CCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATT  >  NC_002947/2087217‑2087515
                                                                                                                                                     |                                                                                                                                                      
ccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGc                                                                                                                                                        >  1:20829/1‑150 (MQ=255)
 caaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCa                                                                                                                                                       >  1:501390/1‑150 (MQ=255)
  aaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCCCTGCAg                                                                                                                                                      >  2:2377712/1‑150 (MQ=255)
    aaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGc                                                                                                                                                    <  1:2029481/150‑1 (MQ=255)
      gACCTTGAAAGGGATGGCCCTTACGGCGGGACCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGccc                                                                                                                                                  >  1:1746046/1‑150 (MQ=255)
        ccTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCACATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCcca                                                                                                                                                >  2:573826/1‑150 (MQ=255)
         cTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccat                                                                                                                                               <  1:2274717/150‑1 (MQ=255)
          ttGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccatc                                                                                                                                              <  2:830080/150‑1 (MQ=255)
          ttGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccatc                                                                                                                                              >  2:359791/1‑150 (MQ=255)
          ttGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAAAAGTCGGCTCCAAGCGTCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccatc                                                                                                                                              >  1:295358/1‑150 (MQ=17)
            gAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccatcca                                                                                                                                            <  1:2075458/150‑1 (MQ=255)
              aaGGGATGGCCCTTACGGCGGGACCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATg                                                                                                                                          >  1:1829435/1‑150 (MQ=255)
                   aTGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATGGATGGGGcc                                                                                                                                     <  2:501390/150‑1 (MQ=255)
                       cccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAg                                                                                                                                >  1:1405718/1‑150 (MQ=255)
                          ttACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                             <  2:320921/150‑1 (MQ=17)
                            aCGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCtt                                                                                                                           <  2:471344/150‑1 (MQ=255)
                              ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCt                                                                                                                                    >  2:1069968/1‑140 (MQ=255)
                              ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCt                                                                                                                                    <  1:1069968/140‑1 (MQ=255)
                              ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                         <  1:929151/150‑1 (MQ=255)
                              ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                         <  1:15272/150‑1 (MQ=255)
                                 gggTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCACCTGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                      <  1:23584/150‑1 (MQ=255)
                                      ccTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                         >  2:250100/1‑127 (MQ=255)
                                      ccTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                         <  1:250100/127‑1 (MQ=255)
                                       cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCCTGGGGCCTCAGCGCTTGACGGGCAtgt                                                                                                                <  2:1842356/150‑3 (MQ=255)
                                       cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCAtgg                                                                                                                >  2:1709854/1‑148 (MQ=255)
                                           ttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc                                                                                                            <  2:1746046/150‑1 (MQ=255)
                                              tCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCtgcag                                                                                                                                <  2:1633720/127‑5 (MQ=17)
                                              tCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCtgcag                                                                                                                                >  1:1633720/1‑123 (MQ=17)
                                                gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc                                                                                                            >  2:756065/1‑145 (MQ=255)
                                                gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc                                                                                                            <  1:756065/145‑1 (MQ=255)
                                                gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                       <  1:2231728/150‑1 (MQ=255)
                                                  ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                     <  1:1709854/150‑2 (MQ=255)
                                                  ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAt                                                                                                     <  1:940498/150‑1 (MQ=255)
                                                    cAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAtct                                                                                                   <  1:1897327/150‑1 (MQ=255)
                                                    cAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAtct                                                                                                   >  2:2092965/1‑150 (MQ=255)
                                                    cAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAtct                                                                                                   >  1:1213889/1‑150 (MQ=255)
                                                     aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                         <  2:1562278/112‑1 (MQ=255)
                                                     aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                         >  1:1562278/1‑112 (MQ=255)
                                                      aaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATctcc                                                                                                 >  1:131959/1‑150 (MQ=255)
                                                          ggAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATctccctcc                                                                                             <  1:2092965/150‑1 (MQ=255)
                                                            aaCCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGc                                                                                           >  2:1255357/1‑150 (MQ=255)
                                                              ccAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCa                                                                                                      <  2:777176/137‑1 (MQ=255)
                                                              ccAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCa                                                                                                      >  1:777176/1‑137 (MQ=255)
                                                               cAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCACCCATGGGTCCTCCGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAt                                                                                        <  2:2393389/150‑1 (MQ=255)
                                                                aaaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATg                                                                                       >  2:595881/1‑150 (MQ=255)
                                                                    aCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTc                                                                                   <  1:96099/150‑1 (MQ=255)
                                                                      cGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                         <  2:1528140/95‑1 (MQ=255)
                                                                      cGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                         >  1:1528140/1‑95 (MQ=255)
                                                                       gTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCAc                                                                                <  1:397582/150‑1 (MQ=255)
                                                                         cGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTc                                                                              <  1:1841152/150‑1 (MQ=255)
                                                                          ggCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCg                                                                             <  2:295358/150‑1 (MQ=255)
                                                                            ctgccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTGGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCCTGGGGCCTCAGCGCTTGACGGGCAGCCATGCCCGTCATCTCCGTCCGCAATGCGTCCACTCGGc                                                                           >  1:61204/4‑150 (MQ=17)
                                                                              cccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcct                                                                         <  2:20829/150‑1 (MQ=255)
                                                                               ccATCATTCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTGCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCAATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctc                                                                        <  1:2301361/150‑1 (MQ=255)
                                                                               ccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctc                                                                        >  1:1262258/1‑150 (MQ=255)
                                                                                catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc                                                                       >  2:2262881/1‑150 (MQ=255)
                                                                                catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc                                                                       >  1:1607098/1‑150 (MQ=255)
                                                                                catcatCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc                                                                       >  1:2047720/1‑150 (MQ=255)
                                                                                 atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                         <  2:1157847/99‑1 (MQ=255)
                                                                                 atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                         >  1:1157847/1‑99 (MQ=255)
                                                                                 atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcct                                                                      >  2:2341772/1‑150 (MQ=255)
                                                                                    atcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGaa                                                                   <  1:2025463/150‑1 (MQ=255)
                                                                                    atcCGGCCCCTACGCTGCGCTCCGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                             <  2:244185/91‑1 (MQ=21)
                                                                                    atcCGGCCCCTACGCTGCGCTCCGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                             >  1:244185/1‑91 (MQ=25)
                                                                                       cGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccat                                                                                                                                               <  2:71103/72‑1 (MQ=255)
                                                                                       cGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccat                                                                                                                                               >  1:71103/1‑72 (MQ=255)
                                                                                       cGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTc                                                                >  1:1430947/1‑150 (MQ=255)
                                                                                               aCGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTg                                                        <  2:1213889/150‑1 (MQ=255)
                                                                                                cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                       <  1:573826/150‑1 (MQ=255)
                                                                                                cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                       <  1:595881/150‑1 (MQ=255)
                                                                                                cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                       <  2:1262258/150‑1 (MQ=255)
                                                                                                cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                       <  2:1910776/150‑1 (MQ=255)
                                                                                                   tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTTCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTgc                                                                                                                         <  2:1791613/81‑2 (MQ=18)
                                                                                                   tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgcgtgc                                                                                                                         >  1:1791613/1‑77 (MQ=18)
                                                                                                      gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                 <  2:2047717/150‑1 (MQ=255)
                                                                                                      gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                 >  1:2169917/1‑150 (MQ=255)
                                                                                                          cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                      >  1:2227855/1‑77 (MQ=255)
                                                                                                          cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCt                                                         >  2:575274/1‑138 (MQ=255)
                                                                                                          cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCt                                                         <  1:575274/138‑1 (MQ=255)
                                                                                                                 cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGACTAATGAAGTCGCAATCTTCCTCACCTGTACTACcgc                                      >  2:1344011/1‑150 (MQ=17)
                                                                                                                   cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTCAcgcgt                                    >  2:2177394/1‑149 (MQ=255)
                                                                                                                   cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgc                                    >  1:1517530/1‑150 (MQ=255)
                                                                                                                      gCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                 >  1:895716/1‑150 (MQ=255)
                                                                                                                        tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATGTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                               >  1:762006/1‑150 (MQ=255)
                                                                                                                        tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                               <  1:884328/150‑1 (MQ=255)
                                                                                                                        tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                               >  2:438724/1‑150 (MQ=255)
                                                                                                                          cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa                             >  2:1436668/1‑150 (MQ=255)
                                                                                                                              cTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGAc                                            <  2:2216325/131‑1 (MQ=255)
                                                                                                                              cTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGAc                                            >  1:2216325/1‑131 (MQ=255)
                                                                                                                               ttGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCaa                        <  2:131959/150‑1 (MQ=255)
                                                                                                                                 gCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGa                      <  1:275681/150‑1 (MQ=255)
                                                                                                                                 gCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGa                      <  2:1517530/150‑1 (MQ=255)
                                                                                                                                   tCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATc                    >  2:2397585/1‑150 (MQ=255)
                                                                                                                                     ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGCATACCGAGCGCt                                 <  2:1844259/135‑1 (MQ=255)
                                                                                                                                     ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                 >  1:1844259/1‑135 (MQ=255)
                                                                                                                                         gAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGc              >  2:1962274/1‑150 (MQ=255)
                                                                                                                                                   tGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa    <  2:1607098/150‑1 (MQ=255)
                                                                                                                                                     cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGc       <  2:2168907/145‑1 (MQ=255)
                                                                                                                                                     cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGc       >  1:2168907/1‑145 (MQ=255)
                                                                                                                                                     cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAtt  <  2:1829435/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATT  >  NC_002947/2087217‑2087515

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: