Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,866,119 G→T 5.1% 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,866,1190GT 5.1% 1326.5 / 41.8 389Q69K (CAG→AAG) gllD4‑oxalomesaconate tautomerase
Reads supporting (aligned to +/- strand):  ref base G (188/181);  new base T (8/12);  total (196/193)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 394 total aligned reads displayed.

GCTCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGG  >  NC_002947/2865970‑2866268
                                                                                                                                                      |                                                                                                                                                     
gCTCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTg                                                                                                                                                       <  1:6181988/150‑1 (MQ=255)
 cTCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGc                                                                                                                                                      >  1:4055110/1‑150 (MQ=255)
  tCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGGCATCGCGCTGCg                                                                                                                                                     >  2:3275009/1‑150 (MQ=255)
     cgccgcTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAg                                                                                                                                                  <  2:3653948/150‑1 (MQ=255)
        cgcTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAgcgc                                                                                                                                               <  2:5385576/150‑1 (MQ=255)
        cgcTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAgcgc                                                                                                                                               <  2:4953957/150‑1 (MQ=255)
        cgcTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTAGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCCCCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAgcgc                                                                                                                                               <  2:3719229/150‑1 (MQ=255)
           tGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCTGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGga                                                                                                                                            <  2:3998066/150‑1 (MQ=255)
           tGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgat                                                                                                                                            <  2:2974558/150‑3 (MQ=255)
           tGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGga                                                                                                                                            <  1:3164213/150‑1 (MQ=255)
            ggcggcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGgat                                                                                                                                           <  2:1508077/150‑1 (MQ=255)
            ggcggcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGgat                                                                                                                                           <  2:2367329/150‑1 (MQ=255)
              cggcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgatgac                                                                                                                                         <  2:6987428/150‑6 (MQ=255)
               ggcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGgatgat                                                                                                                                        >  2:5191390/1‑150 (MQ=255)
               ggcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGgatgat                                                                                                                                        >  2:1202312/1‑150 (MQ=255)
               ggcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGgatgat                                                                                                                                        >  1:3283188/1‑150 (MQ=255)
                gcCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGCCATCATCCGCGCTTCGCAGcgcgatgacg                                                                                                                                        <  1:6383738/150‑7 (MQ=255)
                 ccaccaGGCCGCGCTCTTGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCGATAGTCCACCCTCGCCTCGTCCACCACCACCTGAGCGAACAGGTTGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATgg                                                                                                                                      <  1:3839175/150‑1 (MQ=255)
                 ccaccaGGCCGCGCTCAAGGGCGAACGGCACTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATgg                                                                                                                                      <  1:7848713/150‑1 (MQ=255)
                   accaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCg                                                                                                                                    >  2:2447546/1‑150 (MQ=255)
                    ggtggCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGa                                                                                                                                   >  2:4366448/4‑150 (MQ=255)
                           gcgcTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCt                                                                                                                            <  2:7161222/150‑1 (MQ=255)
                           gcgcTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCt                                                                                                                            <  1:6651248/150‑1 (MQ=255)
                               tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTg                                                                                                                        >  2:7258827/1‑150 (MQ=255)
                                cAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGa                                                                                                                       <  2:6953511/150‑1 (MQ=255)
                                        aaCGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCg                                                                                                               <  1:1202312/150‑1 (MQ=255)
                                          cGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGgc                                                                                                             <  1:8239469/150‑1 (MQ=255)
                                          cGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGgc                                                                                                             >  2:129908/1‑150 (MQ=255)
                                             cccTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGcgcc                                                                                                          >  1:3730515/1‑150 (MQ=255)
                                                tACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCgccgcc                                                                                                       >  2:290157/1‑150 (MQ=255)
                                                tACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCgccgcc                                                                                                       <  2:6208319/150‑1 (MQ=255)
                                                tACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCgccgcc                                                                                                       >  2:5292805/1‑150 (MQ=255)
                                                 aCGCCGGCCAGACTGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGACGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGccgccg                                                                                                      >  2:7801280/1‑150 (MQ=255)
                                                   gCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGAt                                                                                                    <  1:4213445/150‑1 (MQ=255)
                                                   gCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGAt                                                                                                    >  1:8568540/1‑150 (MQ=255)
                                                       gCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCg                                                                                                <  1:6578786/150‑1 (MQ=255)
                                                          aGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCa                                                                                             >  2:4750842/1‑150 (MQ=255)
                                                           gAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCaa                                                                                            <  1:2947290/150‑1 (MQ=255)
                                                            aaTGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAAt                                                                                           <  1:2663632/150‑1 (MQ=255)
                                                               gTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCgccgcc                                                                                                       >  1:5078687/1‑135 (MQ=255)
                                                                  gCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCg                                                                                     <  1:8811422/150‑1 (MQ=255)
                                                                    cGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCggg                                                                                   >  2:5839365/1‑150 (MQ=255)
                                                                    cGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCggg                                                                                   >  2:1679763/1‑150 (MQ=255)
                                                                    cGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCggg                                                                                   >  1:1980033/1‑150 (MQ=255)
                                                                     gCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGc                                                                                  >  2:1749978/1‑150 (MQ=255)
                                                                        gTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATc                                                                               >  2:2403124/1‑150 (MQ=255)
                                                                         ttttGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCg                                                                              >  2:3152371/1‑150 (MQ=255)
                                                                          tttGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCgg                                                                             <  1:1777389/150‑1 (MQ=255)
                                                                            tGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGc                                                                           <  2:5976486/150‑1 (MQ=255)
                                                                             gCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCg                                                                          >  1:1484244/1‑150 (MQ=255)
                                                                              cccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGa                                                                         >  1:4543447/1‑150 (MQ=255)
                                                                               ccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                        >  1:2378015/1‑150 (MQ=255)
                                                                               ccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAg                                                                        >  1:7089093/1‑150 (MQ=255)
                                                                                cATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGc                                                                       >  2:3892733/1‑150 (MQ=255)
                                                                                 aTAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgat                                                                                                                                            >  2:2826190/1‑78 (MQ=255)
                                                                                 aTAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTTCGCAGcgcgat                                                                                                                                            <  1:2826190/80‑3 (MQ=255)
                                                                                 aTAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGcc                                                                      >  1:6939229/1‑150 (MQ=255)
                                                                                     tCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATc                                                                  <  2:3390297/150‑1 (MQ=255)
                                                                                       cACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCAc                                                                >  1:6855737/1‑150 (MQ=255)
                                                                                       cACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCAc                                                                >  1:480886/1‑150 (MQ=255)
                                                                                        aCCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCAcc                                                               >  2:8497563/1‑150 (MQ=255)
                                                                                        aCCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCAcc                                                               <  1:1600135/150‑1 (MQ=255)
                                                                                          ccGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGc                                                             >  1:3972026/1‑150 (MQ=255)
                                                                                                 tCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCagcagc                                                      >  2:1957586/1‑150 (MQ=255)
                                                                                                   gTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCAc                                                    <  1:135542/150‑1 (MQ=255)
                                                                                                   gTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCAc                                                    >  2:1375869/1‑150 (MQ=255)
                                                                                                     ccaccaccaCCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCAccc                                                  >  2:8399785/1‑150 (MQ=255)
                                                                                                      caccaccacCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCg                                                 <  1:7629310/150‑1 (MQ=255)
                                                                                                      caccaccacCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATATGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCg                                                 >  1:4380392/1‑150 (MQ=255)
                                                                                                         caccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTc                                              >  1:3075549/1‑150 (MQ=255)
                                                                                                         caccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTc                                              <  1:353924/150‑1 (MQ=255)
                                                                                                         caccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTc                                              >  1:8009591/1‑150 (MQ=255)
                                                                                                          accaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCa                                             <  2:3531781/150‑1 (MQ=255)
                                                                                                          accaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCa                                             <  2:8588068/150‑1 (MQ=255)
                                                                                                           ccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCAc                                            >  1:8624281/1‑150 (MQ=255)
                                                                                                           ccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCAc                                            <  1:2288650/150‑1 (MQ=255)
                                                                                                              ccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCa                                         <  2:767396/150‑1 (MQ=255)
                                                                                                                tGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGc                                       >  2:654974/1‑150 (MQ=255)
                                                                                                                tGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGc                                       >  1:1045046/1‑150 (MQ=255)
                                                                                                                     gAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGcc                                  >  2:667590/1‑150 (MQ=255)
                                                                                                                     gAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGcc                                  >  1:5247063/1‑150 (MQ=255)
                                                                                                                        cAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAgg                               <  2:4547234/150‑1 (MQ=255)
                                                                                                                         aGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAggg                              <  2:2578867/150‑1 (MQ=255)
                                                                                                                         aGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAggg                              >  1:5712549/1‑150 (MQ=255)
                                                                                                                         aGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAggg                              <  2:8133974/150‑1 (MQ=255)
                                                                                                                          ggTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAgggg                             <  2:5688745/150‑1 (MQ=255)
                                                                                                                            tAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcg                           <  1:2159443/150‑1 (MQ=255)
                                                                                                                                 cACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAg                      <  1:7447127/150‑1 (MQ=255)
                                                                                                                                  aCGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAgg                     >  2:8695137/1‑150 (MQ=255)
                                                                                                                                  aCGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAACACCCGGTCACGCAGCGGGCCAGGGGCGGGCAgg                     <  2:2700010/150‑1 (MQ=255)
                                                                                                                                       ggCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGgtcgtc                <  2:4055110/150‑1 (MQ=255)
                                                                                                                                       ggCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGgtcgtc                <  1:5840879/150‑1 (MQ=255)
                                                                                                                                           tcatccgcgCTTCGCAGCGCGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGc            >  2:6341447/6‑150 (MQ=255)
                                                                                                                                             atccgcgCTTCGCAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCa           >  1:349425/4‑150 (MQ=255)
                                                                                                                                                  cgcTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCaaaaaa      <  2:5495026/150‑1 (MQ=255)
                                                                                                                                                   gcTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAAt     <  1:5537010/150‑1 (MQ=255)
                                                                                                                                                    cTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATa    <  1:7249681/150‑1 (MQ=255)
                                                                                                                                                     tGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGTCAGGTCGTCATGCAAAAAATAg   <  2:7349698/150‑1 (MQ=255)
                                                                                                                                                     tGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAg   <  2:7305072/150‑1 (MQ=255)
                                                                                                                                                      gCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAgg  <  2:2236279/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCTCGCGCCGCTGGCGGCCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGG  >  NC_002947/2865970‑2866268

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