Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,404 C→G 7.1% intergenic (+174/+744) 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,4040CG 7.1% 267.3 / 9.3 84intergenic (+174/+744)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base C (42/36);  new base G (3/3);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAG  >  NC_002947/2087256‑2087551
                                                                                                                                                    |                                                                                                                                                   
cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTTAGGGGGCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc                                                                                                                                                    <  1:1677873/150‑1 (MQ=255)
      gTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                                              <  1:983705/150‑1 (MQ=255)
           ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAt                                                                                                                                         >  1:2188061/1‑150 (MQ=255)
           ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAt                                                                                                                                         <  1:1598375/150‑1 (MQ=255)
           ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAt                                                                                                                                         <  2:2188061/150‑1 (MQ=255)
           ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCCCCGGGGTCCCCTCGCTCCGGCCTTTCGCCCGAGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                         <  1:931587/150‑2 (MQ=17)
                aaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATctccc                                                                                                                                    >  1:1860963/1‑150 (MQ=255)
                    gAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCg                                                                                                                                >  1:131125/1‑150 (MQ=255)
                        cAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAt                                                                                                                            >  2:570598/1‑150 (MQ=255)
                         aaaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCTCGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                           <  2:1509047/134‑1 (MQ=255)
                         aaaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                           >  1:1509047/1‑134 (MQ=255)
                          aaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                         >  1:1163317/1‑134 (MQ=255)
                          aaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                         <  2:1163317/135‑2 (MQ=255)
                          aaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGc                                                                                                                          >  2:667325/1‑150 (MQ=255)
                             aCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTc                                                                                                                       <  1:1294406/150‑1 (MQ=255)
                                gTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCACCCTCCGCAATGCCTCCAc                                                                                                                    >  2:1534607/1‑150 (MQ=255)
                                  cGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTc                                                                                                                  <  1:667325/150‑1 (MQ=255)
                                   ggCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCg                                                                                                                 <  1:570598/150‑1 (MQ=255)
                                   ggCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCg                                                                                                                 <  2:704917/150‑1 (MQ=255)
                                   ggCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCg                                                                                                                 <  2:738258/150‑1 (MQ=255)
                                    gCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCgg                                                                                                                <  1:1300048/150‑1 (MQ=255)
                                        ccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGc                                                                                                                               >  2:2215654/1‑131 (MQ=255)
                                        ccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGc                                                                                                                               <  1:2215654/131‑1 (MQ=255)
                                              tcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAg                                                                                                      <  1:1960353/150‑1 (MQ=255)
                                                    cccTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCaa                                                                                                >  2:1214580/1‑150 (MQ=255)
                                                         cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                                                           >  1:899405/1‑150 (MQ=255)
                                                         cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                                                           <  2:1251984/150‑1 (MQ=255)
                                                               gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                                                     <  2:1860963/150‑1 (MQ=255)
                                                               gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                                                     >  1:1727295/1‑150 (MQ=255)
                                                                cTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGTGTCGCCt                                                                                    <  1:1214580/150‑1 (MQ=255)
                                                                cTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCt                                                                                    >  1:882896/1‑150 (MQ=255)
                                                                        tCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTAcc                                                                            <  2:1727295/150‑1 (MQ=255)
                                                                          cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                                                           <  2:596130/133‑1 (MQ=255)
                                                                          cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                                                           >  1:596130/1‑133 (MQ=255)
                                                                          cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc                                                                          >  2:1087507/1‑150 (MQ=255)
                                                                            cTCGCTCCGGCCTTGCTCCCGGGAGCACTCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgc                                                                        >  1:1280572/1‑150 (MQ=255)
                                                                            cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGCAGTCGCAATCTGCGCCGCATGGCCCACcgcgg                                                                        >  2:723038/1‑149 (MQ=17)
                                                                            cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgc                                                                        >  2:1281210/1‑150 (MQ=255)
                                                                              cGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCgcgcgc                                                                      <  2:131125/150‑1 (MQ=255)
                                                                                cTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCtt                                                                    <  2:899405/150‑1 (MQ=255)
                                                                                 tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTACCTCCGCAATGCCTCCACTCGGCCTCCTGAATTCGCAATCTGCGTCGCCTGGCCTACCGCGCGCtca                                                                   >  1:364645/1‑148 (MQ=255)
                                                                                   cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa                                                                 >  2:1576479/1‑150 (MQ=255)
                                                                                        ttGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCaa                                                            <  2:1527628/150‑1 (MQ=255)
                                                                                        ttGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCaa                                                            <  2:882896/150‑1 (MQ=255)
                                                                                            tCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                         <  1:28296/69‑2 (MQ=11)
                                                                                            tCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                         >  2:28296/1‑68 (MQ=11)
                                                                                            tCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATACCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATc                                                        >  2:2309073/1‑150 (MQ=255)
                                                                                              ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACt                                                                                                                   >  2:2050923/1‑89 (MQ=255)
                                                                                              ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACt                                                                                                                   <  1:2050923/89‑1 (MQ=255)
                                                                                              ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaa                                                      >  1:1998548/1‑150 (MQ=255)
                                                                                              ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTAGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGa                                                          >  1:2219519/1‑146 (MQ=255)
                                                                                              ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTAGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGa                                                          <  2:2219519/146‑1 (MQ=255)
                                                                                               cGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaa                                                     <  1:1087507/150‑1 (MQ=255)
                                                                                                 ggAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAg                                                   <  2:824615/150‑1 (MQ=255)
                                                                                                    gCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCga                                                <  2:1998548/150‑1 (MQ=255)
                                                                                                      aCGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCgaga                                              >  1:381365/1‑150 (MQ=255)
                                                                                                       cgcgCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                                                     <  1:1708421/110‑1 (MQ=255)
                                                                                                       cgcgCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                                                     >  2:1708421/1‑110 (MQ=255)
                                                                                                             tgtggCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAt                                       >  1:1034119/4‑150 (MQ=255)
                                                                                                             gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAt                                       >  1:113561/1‑150 (MQ=255)
                                                                                                                ggCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCg                                    <  2:1034119/150‑1 (MQ=255)
                                                                                                                    ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCgaga                                              <  2:1898053/136‑1 (MQ=255)
                                                                                                                    ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCgaga                                              >  1:1898053/1‑136 (MQ=255)
                                                                                                                      atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATc                              >  2:533995/1‑150 (MQ=255)
                                                                                                                       tccatGGGGCCTCAGCGCTTGACGCGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTCCACTCGGTCTCCTTAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcg                             <  2:777757/150‑1 (MQ=255)
                                                                                                                        ggatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcgg                            >  2:593719/3‑150 (MQ=255)
                                                                                                                           tGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggccg                         >  2:763551/1‑150 (MQ=255)
                                                                                                                              ggCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCa                      >  1:2189683/1‑150 (MQ=255)
                                                                                                                               gCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGaa                                                                >  2:991690/1‑107 (MQ=255)
                                                                                                                               gCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGaa                                                                <  1:991690/107‑1 (MQ=255)
                                                                                                                                cgtCAGCGCTGGACGGGCATCCAGGCCGGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGAATCCCGCGCGCTTAGAAGCAAGGTCAAGGGCGTGTGTGGGATTCGGTGTTCGGCGGGCAcc                    <  2:1780655/148‑1 (MQ=17)
                                                                                                                                 cTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCAccg                   >  1:2060875/1‑150 (MQ=255)
                                                                                                                                  tCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACAACCGCGCGCTTAGAAGCAAGATCAGCATCGAGAAGCGAATTCGGTGATCGGCCGGCAccgc                  >  2:186011/1‑150 (MQ=255)
                                                                                                                                    aGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCAGGCAccgcgg                >  1:1790607/1‑148 (MQ=255)
                                                                                                                                      cgcTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGAt              <  1:2091354/150‑1 (MQ=255)
                                                                                                                                       gcTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATg             >  2:576994/1‑150 (MQ=255)
                                                                                                                                         atGACGGGCATGGATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCa           >  2:312302/2‑150 (MQ=255)
                                                                                                                                            aCGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTc        >  2:466725/1‑150 (MQ=255)
                                                                                                                                             cGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCg       >  1:1245920/1‑150 (MQ=255)
                                                                                                                                               ggCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGc     <  2:1280572/150‑1 (MQ=255)
                                                                                                                                                gCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCg    >  2:2272126/1‑150 (MQ=255)
                                                                                                                                                 cATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGa   <  1:1662623/150‑1 (MQ=255)
                                                                                                                                                  aTCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAg  >  1:1346089/1‑150 (MQ=255)
                                                                                                                                                   tCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCt                                                                                             <  2:2337926/58‑1 (MQ=255)
                                                                                                                                                   tCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCt                                                                                             >  1:2337926/1‑58 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAG  >  NC_002947/2087256‑2087551

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: