Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,151,973 A→G 8.9% Q231R (CAG→CGG) ‡ edd → phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,151,9730AG 8.9% 225.6 / 8.7 87Q231R (CAG→CGG) ‡eddphosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (4/4);  total (43/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 = 9.88e-01

TTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGAAA  >  NC_002947/1151827‑1152122
                                                                                                                                                  |                                                                                                                                                     
ttCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGtt                                                                                                                                                    <  2:1386074/150‑1 (MQ=255)
  ccGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTgg                                                                                                                                                  >  2:1380320/1‑150 (MQ=255)
   cGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGt                                                                                                                                                 >  2:273484/1‑150 (MQ=255)
   cGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGt                                                                                                                                                 >  1:442793/1‑150 (MQ=255)
   cGGCATCTCCAACAAGCAGAAGGCCGATGTGCGACAGCGCTATCGTGCCGGCAAGGCCAGCCGCGAAGAGCGGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGt                                                                                                                                                 >  1:1774884/1‑150 (MQ=255)
    ggCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTg                                                                                                                                                >  1:792561/1‑150 (MQ=255)
     gCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCCGATGGTGa                                                                                                                                               >  2:849738/1‑150 (MQ=255)
     gCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                               >  2:642119/1‑150 (MQ=255)
       aTCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATg                                                                                                                                             >  1:1417362/1‑150 (MQ=255)
           ccAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTgtt                                                                                                                                              >  2:1738213/1‑143 (MQ=255)
           ccAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTgtt                                                                                                                                              <  1:1738213/145‑3 (MQ=255)
            caacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTATGGCACCGCCAACACCAACCAGTTGGTGATGGAAGt                                                                                                                                        >  2:1598698/1‑150 (MQ=255)
            caacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGt                                                                                                                                        <  2:695163/150‑1 (MQ=255)
            caacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGt                                                                                                                                        >  2:1038357/1‑150 (MQ=255)
                aaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                    >  1:1458470/1‑150 (MQ=255)
                aaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                    <  2:926982/150‑1 (MQ=255)
                  gCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATggg                                                                                                                                  <  1:1321387/150‑1 (MQ=255)
                     gAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCt                                                                                                                               >  1:1090128/1‑150 (MQ=255)
                      aaGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTg                                                                                                                              >  1:1924691/1‑150 (MQ=255)
                               gtgCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCg                                                                                                                     <  1:1318422/150‑1 (MQ=255)
                               gtgCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCg                                                                                                                     <  1:1318370/150‑1 (MQ=255)
                                tgCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCgg                                                                                                                    <  2:1881512/150‑1 (MQ=255)
                                  cgcCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGgc                                                                                                                  >  1:437159/1‑150 (MQ=255)
                                         gcTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgt                                                                                                           >  1:1184101/1‑150 (MQ=255)
                                           tATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTATGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                         <  1:1598698/150‑1 (MQ=255)
                                           tATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                         <  2:47376/150‑1 (MQ=255)
                                           tATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                         <  1:626143/150‑1 (MQ=255)
                                             tGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttcgt                                                                                                       <  2:1628088/150‑1 (MQ=255)
                                               ccGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTgttg                                                                                                                                             >  1:1283427/1‑107 (MQ=255)
                                               ccGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTGTTg                                                                                                                                             <  2:1283427/110‑1 (MQ=255)
                                                cGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCaa                                                                                                    >  1:359986/1‑150 (MQ=255)
                                                cGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGCCCTAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGAGCTGCACCTGCCGGGCGCCTCGTTCGTCaa                                                                                                    >  1:626840/1‑150 (MQ=255)
                                                  aaggcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTgttg                                                                                                                                             >  1:1812562/1‑104 (MQ=255)
                                                  aaggcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTGTTg                                                                                                                                             <  2:1812562/107‑1 (MQ=255)
                                                    ggcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCACGGCACCTGCCCCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCg                                                                                                >  2:366741/1‑150 (MQ=255)
                                                     gcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGt                                                                                               <  1:1740526/150‑1 (MQ=255)
                                                     gcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGt                                                                                               <  1:1739996/150‑1 (MQ=255)
                                                           ccAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACAcgc                                                                                         <  1:1331910/150‑1 (MQ=255)
                                                                   gAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGt                                                                                 <  1:727766/150‑1 (MQ=255)
                                                                      gagCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGAc                                                                              <  2:792561/150‑1 (MQ=255)
                                                                       agCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACg                                                                             >  1:23434/1‑150 (MQ=255)
                                                                       agCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACg                                                                             <  2:674381/150‑1 (MQ=255)
                                                                           gTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACt                                                                         <  1:570766/150‑1 (MQ=255)
                                                                              ggAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAc                                                                      <  1:550232/150‑1 (MQ=255)
                                                                               gAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAcc                                                                     <  1:917589/150‑1 (MQ=255)
                                                                                    cGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTCCGTGACGCACTCACCGCCGa                                                                >  2:578485/1‑150 (MQ=255)
                                                                                     gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAg                                                               >  2:478561/1‑150 (MQ=255)
                                                                                          gAATTCCTACCACAGCGCCGGCACCAGCACCTTCTACGGCACCGCCAACACCAACCGGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTgc                                                          <  2:712109/150‑1 (MQ=255)
                                                                                          gAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGa                                                                >  1:646064/1‑144 (MQ=255)
                                                                                          gAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGa                                                                <  2:646064/144‑1 (MQ=255)
                                                                                           aaTTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTgcg                                                         <  2:1924691/150‑1 (MQ=255)
                                                                                           aaTGCCTACCACAGACCCGGAACACGCACCTTCTACGGCCCCGACAAGACCACCAAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTgcg                                                         <  2:1841893/150‑1 (MQ=255)
                                                                                             ttCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgca                                                       >  2:1460008/1‑150 (MQ=255)
                                                                                              tCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgcag                                                      >  1:1611940/1‑150 (MQ=255)
                                                                                                 tACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGcagcag                                                   >  2:329605/1‑150 (MQ=255)
                                                                                                 tACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGACCTCACCGCCGAGACTGCgcatgag                                                   >  2:1543942/1‑146 (MQ=255)
                                                                                                  aCCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAgg                                                  <  2:1774884/150‑1 (MQ=255)
                                                                                                   ccACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGt                                                 <  1:1380320/150‑1 (MQ=255)
                                                                                                     acaGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGa                                               >  1:21674/1‑150 (MQ=255)
                                                                                                          cccGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAg                                                               <  2:318570/129‑1 (MQ=255)
                                                                                                          cccGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAg                                                               >  1:318570/1‑129 (MQ=255)
                                                                                                          cccGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACgcgc                                          <  2:1724017/150‑1 (MQ=255)
                                                                                                           ccGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCa                                         >  1:288279/1‑150 (MQ=255)
                                                                                                            cGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCAt                                        <  2:1806883/150‑1 (MQ=255)
                                                                                                              gcacctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACgcgcgcga                                      <  2:1169183/150‑5 (MQ=255)
                                                                                                              gcacctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGa                                      <  1:966629/150‑1 (MQ=255)
                                                                                                                  ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGt                                                                                 >  1:1080951/1‑103 (MQ=255)
                                                                                                                  ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGt                                                                                 <  2:1080951/103‑1 (MQ=255)
                                                                                                                  ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGAccca                                  >  1:268245/1‑148 (MQ=255)
                                                                                                                  ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCaa                                  <  1:1426561/150‑1 (MQ=255)
                                                                                                                   tgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                                 <  1:1490711/150‑1 (MQ=255)
                                                                                                                    gcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTac                                                                                             >  1:187562/1‑89 (MQ=255)
                                                                                                                    gcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTac                                                                                             <  2:187562/89‑1 (MQ=255)
                                                                                                                       cctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCa                             >  2:1505176/1‑150 (MQ=255)
                                                                                                                             aCGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCa                       >  2:1450635/1‑150 (MQ=255)
                                                                                                                               ggCACCGCCCACACCAACCAGTTGGCGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGc                     >  2:642138/1‑150 (MQ=255)
                                                                                                                                  aCCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCGGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTc                  >  1:323667/1‑150 (MQ=255)
                                                                                                                                  aCCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTc                  >  2:1118317/1‑150 (MQ=255)
                                                                                                                                    cGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCAt                <  2:1900262/150‑1 (MQ=255)
                                                                                                                                      ccAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGc              >  2:247595/1‑150 (MQ=255)
                                                                                                                                      ccAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGc              >  2:686213/1‑150 (MQ=255)
                                                                                                                                      ccAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGc              <  2:1458470/150‑1 (MQ=255)
                                                                                                                                      ccAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGc              <  2:442793/150‑1 (MQ=255)
                                                                                                                                      ccAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTAATGc              <  2:1169936/150‑1 (MQ=255)
                                                                                                                                        atcacCAACTGGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCg            >  2:633196/3‑150 (MQ=255)
                                                                                                                                         acacCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGc           <  2:94579/150‑1 (MQ=255)
                                                                                                                                         acacCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGc           <  1:1118317/150‑1 (MQ=255)
                                                                                                                                           accaactgGTTGGTGATGGAAGTGATGGGGCTGCAGCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTg         <  1:633196/150‑1 (MQ=255)
                                                                                                                                           accaaccaGTTGGTGATGGAGGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTg         >  1:1832182/1‑150 (MQ=255)
                                                                                                                                           accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGc                     <  1:73902/138‑1 (MQ=255)
                                                                                                                                           accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGc                     >  2:73902/1‑138 (MQ=255)
                                                                                                                                           accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCAGCGGCGTgtacg                                                                            <  2:1405988/83‑5 (MQ=255)
                                                                                                                                           accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCAGCGGCGTgtacg                                                                            >  1:1405988/1‑79 (MQ=255)
                                                                                                                                            ccaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTgg        <  2:1184101/150‑1 (MQ=255)
                                                                                                                                              aaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGc      <  2:437159/150‑1 (MQ=255)
                                                                                                                                               accaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCg     >  1:616743/1‑150 (MQ=255)
                                                                                                                                                 caGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCTaa   <  1:732747/150‑1 (MQ=255)
                                                                                                                                                 caGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGaa   >  2:1498416/1‑150 (MQ=255)
                                                                                                                                                  aGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGaaa  >  2:1906861/1‑150 (MQ=255)
                                                                                                                                                  aGTTGGTGATGGAAGTGAATCGTCTGCACCCGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGaaa  <  1:1038357/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGAAA  >  NC_002947/1151827‑1152122

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