Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,861,170 G→T 7.9% Q69K (CAG→AAG)  gllD ← 4‑oxalomesaconate tautomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,861,1700GT 7.9% 247.4 / 7.3 83Q69K (CAG→AAG) gllD4‑oxalomesaconate tautomerase
Reads supporting (aligned to +/- strand):  ref base G (36/40);  new base T (3/4);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CTCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAG  >  NC_002947/2861022‑2861318
                                                                                                                                                     |                                                                                                                                                    
cTCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCCACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGc                                                                                                                                                     >  2:724457/1‑150 (MQ=255)
     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                                                                                                                        <  2:955492/150‑1 (MQ=255)
                             cTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGt                                                                                                                        <  1:606795/150‑1 (MQ=255)
                             cTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGt                                                                                                                        <  1:1793662/150‑1 (MQ=255)
                              tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgat                                                                                                                                           >  2:1787006/1‑128 (MQ=255)
                              tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgat                                                                                                                                           <  1:1787006/130‑3 (MQ=255)
                              tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTg                                                                                                                          >  2:366428/1‑147 (MQ=255)
                              tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTg                                                                                                                          <  1:366428/147‑1 (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                                                                                                                          <  2:156580/131‑1 (MQ=255)
                                              cTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTg                                                                                                                          >  1:156580/1‑131 (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                                                                        <  1:1132597/150‑1 (MQ=255)
                                                                             cccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGa                                                                        <  2:1062716/150‑1 (MQ=255)
                                                                              ccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                       >  2:523287/1‑150 (MQ=255)
                                                                              ccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                       >  1:88833/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                                                                       <  2:980905/139‑1 (MQ=255)
                                                                                         ccGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                       >  1:980905/1‑139 (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:1827375/1‑150 (MQ=255)
                                                                                                     caccaccacCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCg                                                <  1:182514/150‑1 (MQ=255)
                                                                                                      accaccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCgg                                               >  1:562033/1‑150 (MQ=255)
                                                                                                      accaccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCgg                                               <  2:1492744/150‑1 (MQ=255)
                                                                                                         accaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCa                                            <  1:1496401/150‑1 (MQ=255)
                                                                                                                 aGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCgg                                    >  2:1713202/1‑150 (MQ=255)
                                                                                                                 aGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCgg                                    <  1:724457/150‑1 (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                         >  1:879461/1‑150 (MQ=255)
                                                                                                                            aGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcgg                         <  2:1728857/150‑1 (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                      <  2:1596289/130‑1 (MQ=255)
                                                                                                                                                    tGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCa                      >  1:1596289/1‑130 (MQ=255)
                                                                                                                                                    tGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGACATTCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAg  <  1:118873/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAG  >  NC_002947/2861022‑2861318

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: