Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,307,447 C→G 9.1% A59A (GCG→GCC livK ← branched‑chain amino acids ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,307,4470CG 9.1% 270.1 / 15.8 88A59A (GCG→GCClivKbranched‑chain amino acids ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (43/37);  new base G (4/4);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

GCTGGCTGGGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAGTACCA  >  NC_002947/1307301‑1307592
                                                                                                                                                  |                                                                                                                                                 
gctggctgGGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCcgcg                                                                                                                                                >  1:944028/1‑150 (MQ=255)
  tggctggGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCgcggc                                                                                                                                              >  1:823531/1‑150 (MQ=255)
    gctggGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGt                                                                                                                                            <  1:640138/150‑1 (MQ=255)
        ggTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCTTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGat                                                                                                                                        <  1:461088/150‑1 (MQ=255)
        ggTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCAACGCCGCGGGCGGCGTGGat                                                                                                                                        <  2:145939/150‑1 (MQ=255)
          tggagctggagcACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCAACGCCGCGGGCGGCGTGGat                                                                                                                                        <  1:573832/148‑1 (MQ=255)
          tggagctggagcAAAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCAACGCCGCGGGCGGCGTGGat                                                                                                                                        >  2:573832/1‑148 (MQ=255)
              gctggagcACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCtt                                                                                                                                  >  2:1454313/1‑150 (MQ=255)
              gctggagcACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCtt                                                                                                                                  >  1:621521/1‑150 (MQ=255)
                tggagcACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCTTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGa                                                                                                                                <  1:631645/150‑1 (MQ=255)
                  gagcACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATg                                                                                                                              >  2:1702471/1‑150 (MQ=255)
                        aGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGa                                                                                                                                <  1:322604/142‑1 (MQ=255)
                        aGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGa                                                                                                                                >  2:322604/1‑142 (MQ=255)
                          gTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGgcttgctt                                                                                                                      >  2:159287/1‑150 (MQ=255)
                          gTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGgcttcctt                                                                                                                      >  1:249549/1‑150 (MQ=255)
                             gCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGgc                                                                                                                   >  1:1817834/1‑150 (MQ=255)
                               cGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGgcgc                                                                                                                 <  1:2184957/150‑1 (MQ=255)
                                 aCCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCGGGCGGCGt                                                                                                                                            <  2:1522500/121‑1 (MQ=255)
                                 aCCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCGGGCGGCGt                                                                                                                                            >  1:1522500/1‑121 (MQ=255)
                                  ccACGAACTTGACGCCGTCGTTGACCACTCTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTCGCTTGGCACCGa                                                                                                              >  2:10235/1‑150 (MQ=255)
                                   cACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGAt                                                                                                             >  1:1671207/1‑150 (MQ=255)
                                   cACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGAt                                                                                                             >  2:1825581/1‑150 (MQ=255)
                                   cACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGAt                                                                                                             >  1:716080/1‑150 (MQ=255)
                                          ttGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGc                                                                                                      <  1:1496852/150‑1 (MQ=255)
                                          ttGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGc                                                                                                      <  1:318486/150‑1 (MQ=255)
                                           tGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCa                                                                                                     <  1:1825581/150‑1 (MQ=255)
                                             aCGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATg                                                                                                   <  1:2103478/150‑1 (MQ=255)
                                                   tcgtTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCg                                                                                             <  2:31235/150‑1 (MQ=255)
                                                      ttGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATgg                                                                                                                             <  2:729816/115‑1 (MQ=255)
                                                      ttGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATgg                                                                                                                             >  1:729816/1‑115 (MQ=255)
                                                         aCCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTgg                                                                                       >  1:584986/1‑150 (MQ=255)
                                                         aCCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTgg                                                                                       <  2:1817834/150‑1 (MQ=255)
                                                            aCTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCCTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTc                                                                                    <  2:2070937/150‑1 (MQ=255)
                                                             cTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCa                                                                                   >  1:665838/1‑150 (MQ=255)
                                                              tttGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTGCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCAc                                                                                  >  1:1336455/1‑150 (MQ=255)
                                                                tgttgGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACgg                                                                                <  2:1400561/150‑1 (MQ=255)
                                                                 gttgGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACggg                                                                               >  1:788407/1‑150 (MQ=255)
                                                                 gttgGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACggg                                                                               <  1:1499485/150‑1 (MQ=255)
                                                                  ttgGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGa                                                                              >  2:35952/1‑150 (MQ=255)
                                                                       gACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGt                                                                         <  2:2032984/150‑1 (MQ=255)
                                                                          tGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGCTGTCTTTCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCtt                                                                      <  1:1454313/150‑1 (MQ=255)
                                                                          tGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACGTGGTCACGGGACCCGTCtt                                                                      >  2:1535376/1‑150 (MQ=255)
                                                                           gagaCGGCCTGTTGAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCTCCCGCGGCGTTGATATCCTTGATTGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCCCGGGACCCGTCTTg                                                                     <  2:249549/147‑1 (MQ=255)
                                                                           gCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTg                                                                     <  2:1609480/150‑1 (MQ=255)
                                                                           gCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTg                                                                     <  1:1656200/150‑1 (MQ=255)
                                                                            ccACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTg                                                                                        <  1:2030351/130‑1 (MQ=255)
                                                                            ccACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTg                                                                                        >  2:2030351/1‑130 (MQ=255)
                                                                            ccACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTgg                                                                    <  2:944028/150‑1 (MQ=255)
                                                                            ccACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTgg                                                                    >  1:148112/1‑150 (MQ=255)
                                                                              aCGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGa                                                                  >  1:1673458/1‑150 (MQ=255)
                                                                                 gCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCg                                                               <  2:148112/150‑1 (MQ=255)
                                                                                   cTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGc                                                             <  2:665838/150‑1 (MQ=255)
                                                                                    tGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCg                                                            >  2:118376/1‑150 (MQ=255)
                                                                                           ggTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGa                                                     >  1:77553/1‑150 (MQ=255)
                                                                                           ggTCGCAGGCGTCGTCGTATTCATTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGa                                                     <  2:1973228/150‑1 (MQ=255)
                                                                                            gTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGat                                                                                                                                        >  2:1864951/1‑66 (MQ=255)
                                                                                            gTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGat                                                                                                                                        <  1:1864951/66‑1 (MQ=255)
                                                                                               gCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCtt                                                 >  2:715540/1‑150 (MQ=255)
                                                                                                 aGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGa                                               <  1:715540/150‑1 (MQ=255)
                                                                                                  ggCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGAt                                              <  1:1702471/150‑1 (MQ=255)
                                                                                                   gCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATg                                             <  2:788407/150‑1 (MQ=255)
                                                                                                         tcgtATTCCTTGGCTTCGAGCATCTTGCCATCCACGCGGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATcg                                       >  1:1892650/1‑150 (MQ=255)
                                                                                                         tcgtATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATcg                                       >  1:1147717/1‑150 (MQ=255)
                                                                                                          cgtATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTCGCTATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATcgg                                      <  1:512433/150‑1 (MQ=255)
                                                                                                            tATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGAccc                                                                           <  1:1464795/111‑1 (MQ=255)
                                                                                                            tATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGAccc                                                                           >  2:1464795/1‑111 (MQ=255)
                                                                                                             aTTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGa                                                                                                              >  1:1237222/1‑75 (MQ=255)
                                                                                                             aTTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGa                                                                                                              <  2:1237222/75‑1 (MQ=255)
                                                                                                                  ttGGCTTCGAGCATCTTGCCATCCACGCCGCCCGAGGAGTTGCTCTAACTGCACGCTTGCTTGGCGCCGGGCAACGGCATGTCGCCGTACTGGGTCAGGGGAACCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAAc                              >  1:1212786/1‑150 (MQ=255)
                                                                                                                     gCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGaa                           >  2:8189/1‑150 (MQ=255)
                                                                                                                      cTTCTAGCAGCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAAt                          >  2:1203243/1‑150 (MQ=255)
                                                                                                                       ttCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAAtg                         <  1:159287/150‑1 (MQ=255)
                                                                                                                       ttCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAAtg                         >  1:806281/1‑150 (MQ=255)
                                                                                                                       ttCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAAtg                         >  2:1635565/1‑150 (MQ=255)
                                                                                                                         cGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAAtggc                       <  1:118376/150‑1 (MQ=255)
                                                                                                                              aTCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATggctggct                  >  1:1929124/1‑150 (MQ=255)
                                                                                                                              aTCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATggctggct                  >  1:228942/1‑150 (MQ=255)
                                                                                                                                  tGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCgg                                                              <  1:469045/102‑1 (MQ=255)
                                                                                                                                  tGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCgg                                                              >  2:469045/1‑102 (MQ=255)
                                                                                                                                   gataTCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTAcccc             >  1:2093390/4‑150 (MQ=255)
                                                                                                                                     cATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGTACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTAccc              <  2:1957842/147‑1 (MQ=255)
                                                                                                                                     cATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTAccc              >  1:1957842/1‑147 (MQ=255)
                                                                                                                                      aTCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGCTGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCgcc          >  2:1974252/1‑148 (MQ=255)
                                                                                                                                       tcaaCGCCGCGGGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGAt                                              >  2:2336102/4‑113 (MQ=255)
                                                                                                                                       tcaaCGCCGCGGGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGAt                                              <  1:2336102/110‑1 (MQ=255)
                                                                                                                                       tCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGcc         >  1:1259775/1‑150 (MQ=255)
                                                                                                                                       tCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGCGCTCCATCAGGGGACCGGCGATGCCGAGCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGcc         >  1:145251/1‑150 (MQ=255)
                                                                                                                                         aaCGCCGCGGGCGGCGTTGATATCCTTGATGGCTGGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCAAGGGACCCGGCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAg       >  1:1355129/2‑150 (MQ=255)
                                                                                                                                            gccgccCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAGTAc    <  2:1892650/150‑1 (MQ=255)
                                                                                                                                             ccgcccgcGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATcggcg                                    >  1:1433742/1‑117 (MQ=255)
                                                                                                                                             ccgcccgcGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATcggcg                                    <  2:1433742/117‑1 (MQ=255)
                                                                                                                                              cgcccgcGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAGTACCa  >  1:957744/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGGCTGGGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAGTACCA  >  NC_002947/1307301‑1307592

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: