Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,861,173 A→C 8.3% S68A (TCG→GCG)  gllD ← 4‑oxalomesaconate tautomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,861,1730AC 8.3% 216.7 / 10.3 82S68A (TCG→GCG) gllD4‑oxalomesaconate tautomerase
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base C (3/4);  total (37/45)
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.86e-01

GCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGGCGC  >  NC_002947/2861027‑2861322
                                                                                                                                                   |                                                                                                                                                     
gccgcTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAgc                                                                                                                                                    >  1:1728857/1‑150 (MQ=255)
  cgcTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGACCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAgcgc                                                                                                                                                  <  1:21551/150‑1 (MQ=255)
      ggcggcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgatg                                                                                                                                              <  2:1352610/150‑4 (MQ=255)
          gcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGatgatg                                                                                                                                          <  1:1182894/150‑1 (MQ=255)
            caccaGGCAGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGc                                                                                                                                        >  2:219540/1‑150 (MQ=255)
             accaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAgc                                                                                                                                                    <  1:477077/137‑1 (MQ=255)
             acaaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAgc                                                                                                                                                    >  2:477077/4‑137 (MQ=255)
                aGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACt                                                                                                                                    <  1:1346918/150‑1 (MQ=255)
                 ggCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACtt                                                                                                                                   <  2:1756950/150‑1 (MQ=255)
                  gCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACttt                                                                                                                                  >  2:217960/1‑150 (MQ=255)
                       gcTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTgg                                                                                                                             <  2:1073295/150‑1 (MQ=255)
                        cTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGt                                                                                                                            <  1:606795/150‑1 (MQ=255)
                        cTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGt                                                                                                                            <  1:1793662/150‑1 (MQ=255)
                        cTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGt                                                                                                                            <  2:955492/150‑1 (MQ=255)
                         tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgat                                                                                                                                               >  2:1787006/1‑128 (MQ=255)
                         tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgat                                                                                                                                               <  1:1787006/130‑3 (MQ=255)
                         tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTg                                                                                                                              <  1:366428/147‑1 (MQ=255)
                         tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTg                                                                                                                              >  2:366428/1‑147 (MQ=255)
                          cAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGa                                                                                                                          <  2:1083526/150‑1 (MQ=255)
                           aaGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAg                                                                                                                         >  1:1053002/1‑150 (MQ=255)
                           aaGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAg                                                                                                                         >  1:1052297/1‑150 (MQ=255)
                            aGGTCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGc                                                                                                                        <  2:724359/150‑1 (MQ=255)
                              ggCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGACTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGa                                                                                                                      >  1:525054/1‑150 (MQ=255)
                              ggCGAACGGCCCTACGCCGGCCAGAATGCTGATGGATTAGGTCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGa                                                                                                                      <  2:1589362/150‑1 (MQ=255)
                               gCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAg                                                                                                                     >  2:1003360/1‑150 (MQ=255)
                                   aCGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCgg                                                                                                                 >  2:1779335/1‑150 (MQ=255)
                                         cTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTg                                                                                                                              >  1:156580/1‑131 (MQ=255)
                                         cTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTg                                                                                                                              <  2:156580/131‑1 (MQ=255)
                                         cTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGcgccgc                                                                                                           <  1:999357/150‑1 (MQ=255)
                                            cGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGa                                                                                                        <  2:745145/150‑1 (MQ=255)
                                             gCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGaa                                                                                                                                                      <  2:1247241/103‑1 (MQ=255)
                                             gCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGaa                                                                                                                                                      >  1:1247241/1‑103 (MQ=255)
                                                 gCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGCCATCATCCGCGCTTCGCAGcgcgat                                                                                                                                               >  1:1786902/1‑105 (MQ=37)
                                                 gCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGCCATCATCCGCGCTTCGCAGcgcgat                                                                                                                                               <  2:1786902/107‑3 (MQ=37)
                                                      agtGTTGCCGCAGTTGTGCCCATAGTCCACCCGCGCCTCGGCCACCCCCACCTGAGCGAACAGGTAGTCCACGTCGGCGCCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTTTTGAGCGAGTCGGCGCCGCCGATGCCGTCAAt                                                                                              <  1:909848/148‑1 (MQ=255)
                                                      aaTGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAAt                                                                                              <  1:1142027/150‑1 (MQ=255)
                                                             ccGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCgg                                                                                       >  1:1639709/1‑150 (MQ=255)
                                                                    tttGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCgg                                                                                <  2:467769/150‑1 (MQ=255)
                                                                        cccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGa                                                                            <  2:1062716/150‑1 (MQ=255)
                                                                        cccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGa                                                                            <  1:1132597/150‑1 (MQ=255)
                                                                         ccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                           >  1:88833/1‑150 (MQ=255)
                                                                         ccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                           >  2:523287/1‑150 (MQ=255)
                                                                          cATAGCCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGc                                                                          >  2:1139396/1‑150 (MQ=255)
                                                                              gTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCAt                                                                      >  2:167581/1‑150 (MQ=255)
                                                                               tCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGCTGAAGGGCTCGTGGCGCGGGGGCTAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATc                                                                     >  1:431971/1‑150 (MQ=255)
                                                                                ccACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGCCGCCCGAAAACGGGCGACTTTGGTGGTGAGCGATTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCa                                                                    <  1:232401/150‑1 (MQ=255)
                                                                                   cccGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCCTCAAGAAGCGAAGCGCGTGGGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCg                                                                 >  1:1672831/1‑150 (MQ=255)
                                                                                    ccGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                           >  1:980905/1‑139 (MQ=255)
                                                                                    ccGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                           <  2:980905/139‑1 (MQ=255)
                                                                                    ccGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGc                                                                >  1:850346/1‑150 (MQ=255)
                                                                                      gcgcCTCTTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCca                                                              >  1:1492744/1‑150 (MQ=255)
                                                                                         ccTCGTCCACCACCACCTGAGCGACCAGGTAGTCCACGTAGGCGTCATCTCTCTTCGAAGCGCGGATGATGGCGGCTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCcagca                                                           <  1:1068434/150‑1 (MQ=255)
                                                                                         ccTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCcagca                                                           >  1:268909/1‑150 (MQ=255)
                                                                                          cTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCcagcag                                                          <  1:167581/150‑1 (MQ=255)
                                                                                          cTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCcagcag                                                          >  2:171197/1‑150 (MQ=255)
                                                                                                caccaccacCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCg                                                    <  1:182514/150‑1 (MQ=255)
                                                                                                caccaccacCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCg                                                    >  1:1827375/1‑150 (MQ=255)
                                                                                                 accaccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCgg                                                   <  2:1492744/150‑1 (MQ=255)
                                                                                                 accaccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCgg                                                   >  1:562033/1‑150 (MQ=255)
                                                                                                    accaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCa                                                <  1:1496401/150‑1 (MQ=255)
                                                                                                            aGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCgg                                        <  1:724457/150‑1 (MQ=255)
                                                                                                            aGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCgg                                        >  2:1713202/1‑150 (MQ=255)
                                                                                                             gCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCggg                                       <  1:346719/150‑1 (MQ=255)
                                                                                                               gAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGcc                                     <  1:1003360/150‑1 (MQ=255)
                                                                                                               gAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGcc                                     >  2:306034/1‑150 (MQ=255)
                                                                                                                aaCAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCa                                    >  1:722874/1‑150 (MQ=255)
                                                                                                                  cAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAgg                                  <  2:1639709/150‑1 (MQ=255)
                                                                                                                      tAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcg                              <  1:217960/150‑1 (MQ=255)
                                                                                                                       aGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcgg                             <  2:1728857/150‑1 (MQ=255)
                                                                                                                       aGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcgg                             >  1:879461/1‑150 (MQ=255)
                                                                                                                         tCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGtggcgggc                           >  1:1450067/1‑150 (MQ=255)
                                                                                                                         tCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcgggc                           >  1:1029011/1‑150 (MQ=255)
                                                                                                                           cACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAg                         <  1:587119/150‑1 (MQ=255)
                                                                                                                              gTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGgtc                      <  2:722874/150‑1 (MQ=255)
                                                                                                                                   catcatccgcgCTTCGCAGCGCGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGCCGTCAt                 >  2:715078/8‑150 (MQ=255)
                                                                                                                                   cGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCCGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCAt                 <  2:525054/150‑1 (MQ=255)
                                                                                                                                       aTCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCACGGGCGGGCAGGTCGTCATGCaa             >  2:958039/1‑150 (MQ=255)
                                                                                                                                        tCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCaaa            >  2:236086/1‑150 (MQ=255)
                                                                                                                                              cTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATa       <  1:1159290/150‑1 (MQ=255)
                                                                                                                                               tGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCa                          >  1:1596289/1‑130 (MQ=255)
                                                                                                                                               tGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCa                          <  2:1596289/130‑1 (MQ=255)
                                                                                                                                               tGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGACATTCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAg      <  1:118873/150‑1 (MQ=255)
                                                                                                                                                 cGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGgc    <  1:734253/150‑1 (MQ=255)
                                                                                                                                                 cGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGgc    <  2:1875446/150‑1 (MQ=255)
                                                                                                                                                   aaGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGgcgc  <  1:719244/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGGCGC  >  NC_002947/2861027‑2861322

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