Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,177,743 C→A 35.9% D85Y (GAT→TAT)  thiD ← bifunctional hydroxy‑methylpyrimidine kinase/ hydroxy‑phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,177,7430CA35.9% 446.4 / 250.7 270D85Y (GAT→TAT) thiDbifunctional hydroxy‑methylpyrimidine kinase/ hydroxy‑phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (105/67);  new base A (50/47);  total (155/115)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

GAAACCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGC  >  CP009273/2177595‑2177892
                                                                                                                                                    |                                                                                                                                                     
gAAACCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCg                                                                                                                                                      >  2:1279916/1‑150 (MQ=255)
 aaaCCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCgg                                                                                                                                                     <  2:4371919/150‑1 (MQ=255)
    ccTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGttt                                                                                                                                                  >  1:80683/1‑150 (MQ=255)
      tgtgGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTcc                                                                                                                                                >  2:1837757/1‑150 (MQ=255)
       gtgGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCg                                                                                                                                               <  2:1514456/150‑1 (MQ=255)
          gCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCa                                                                                                                                            <  1:629589/150‑1 (MQ=255)
          gCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCa                                                                                                                                            <  1:1894106/150‑1 (MQ=255)
            aataatCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGc                                                                                                                                          >  1:2620845/1‑150 (MQ=255)
             ataatCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCa                                                                                                                                         >  1:2111568/1‑150 (MQ=255)
             ataatCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCa                                                                                                                                         >  1:4616431/1‑150 (MQ=255)
              tactCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCAt                                                                                                                                        <  1:2818964/147‑1 (MQ=255)
               aatCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATa                                                                                                                                       >  1:4194874/1‑150 (MQ=255)
               aatCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATa                                                                                                                                       >  1:3862188/1‑150 (MQ=255)
                atCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATAc                                                                                                                                      >  1:154432/1‑150 (MQ=255)
                atCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTGCAACAATATAGGTTTCCGGCAGCATag                                                                                                                                      >  2:4614174/1‑149 (MQ=255)
                 tCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATAcc                                                                                                                                     <  1:101141/150‑1 (MQ=255)
                   gACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGa                                                                                                                                   <  1:4393375/150‑1 (MQ=255)
                   gACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGa                                                                                                                                   >  1:2740058/1‑150 (MQ=255)
                   gACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGa                                                                                                                                   >  1:2601922/1‑150 (MQ=255)
                    aCTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGAt                                                                                                                                  <  2:2718643/150‑1 (MQ=255)
                      tGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGAttt                                                                                                                                >  1:1444495/1‑150 (MQ=255)
                          cAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGt                                                                                                                            <  1:3685623/150‑1 (MQ=255)
                              gTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTa                                                                                                                        >  2:3200077/1‑150 (MQ=255)
                              gTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTa                                                                                                                        <  1:4429546/150‑1 (MQ=255)
                               tAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTAt                                                                                                                       >  1:224188/1‑150 (MQ=255)
                                aGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATc                                                                                                                      <  1:171094/150‑1 (MQ=255)
                                aGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATc                                                                                                                      >  1:2306829/1‑150 (MQ=255)
                                 gCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCg                                                                                                                     <  1:1056517/150‑1 (MQ=255)
                                 gCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCg                                                                                                                     <  1:1164263/150‑1 (MQ=255)
                                   aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGAt                                                                                                                   >  1:3548876/1‑150 (MQ=255)
                                      cgccgAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCg                                                                                                                >  1:2059845/1‑150 (MQ=255)
                                                gAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTg                                                                                                      >  2:2605894/1‑150 (MQ=255)
                                                 aaaGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGa                                                                                                     >  1:2996352/1‑150 (MQ=255)
                                                     cagcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAacac                                                                                                 >  2:1607691/1‑150 (MQ=255)
                                                     cagcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAacac                                                                                                 >  2:1607785/1‑150 (MQ=255)
                                                      agcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACAcc                                                                                                >  1:748880/1‑150 (MQ=255)
                                                        cGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGa                                                                                              <  1:3968326/150‑1 (MQ=255)
                                                            tCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCg                                                                                          >  1:1760275/1‑150 (MQ=255)
                                                            tCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCg                                                                                          >  1:4083659/1‑150 (MQ=255)
                                                              gccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAg                                                                                        >  1:516721/1‑150 (MQ=255)
                                                                cgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCCCCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCt                                                                                      <  2:1683405/150‑1 (MQ=255)
                                                                cgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCt                                                                                      >  2:3578580/1‑150 (MQ=255)
                                                                    ttttGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGggc                                                                                  >  1:2041818/1‑150 (MQ=255)
                                                                    ttttGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGggc                                                                                  <  1:3341183/150‑1 (MQ=255)
                                                                        gCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGgcggcg                                                                              <  1:2337670/150‑1 (MQ=255)
                                                                           aGCATAACGGTGTCGAGTACCACGTTTTGGNTCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACa                                                                           >  1:4280322/1‑150 (MQ=255)
                                                                              aTAACGGTGTCGAGTACCACGTTTTGGATCTGATACCGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACaaaa                                                                        <  2:2145963/150‑1 (MQ=255)
                                                                              aTAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACaaaa                                                                        <  2:2266553/150‑1 (MQ=255)
                                                                               tAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAAt                                                                       >  1:865996/1‑150 (MQ=255)
                                                                                     tgtCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCt                                                                 <  1:3236186/150‑1 (MQ=255)
                                                                                     tgtCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCt                                                                 >  2:3743951/1‑150 (MQ=255)
                                                                                         gAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAAt                                                             >  2:3989070/1‑150 (MQ=255)
                                                                                               cACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATa                                                       <  1:2202690/150‑1 (MQ=255)
                                                                                               cACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATa                                                       >  2:1152564/1‑150 (MQ=255)
                                                                                                 cGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATaca                                                     >  1:564312/1‑150 (MQ=255)
                                                                                                    tttGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGAGTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCg                                                  <  2:4048320/150‑1 (MQ=255)
                                                                                                     ttGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGa                                                 >  2:2490467/1‑150 (MQ=255)
                                                                                                     ttGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGa                                                 >  2:2696431/1‑150 (MQ=255)
                                                                                                      tGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGAc                                                >  1:1324533/1‑150 (MQ=255)
                                                                                                      tGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGAc                                                <  1:17306/150‑1 (MQ=255)
                                                                                                      gggAGCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGATTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGAc                                                <  2:2815728/149‑1 (MQ=255)
                                                                                                       ggATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACt                                               <  1:3020947/150‑1 (MQ=255)
                                                                                                       ggATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACt                                               >  2:2906562/1‑150 (MQ=255)
                                                                                                         aTCTGATAACGTTGCATCCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCGAGCATACAGATTTTTGTGGTATCGATTCGCACATCGCTGAACAGCGAATCGAGCTGGGCGGCGAAAAAATCAGGCTCAATGCGATACACCGACTGt                                             <  2:2379845/150‑1 (MQ=255)
                                                                                                         aTCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGt                                             >  2:357003/1‑150 (MQ=255)
                                                                                                         aTCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGt                                             <  2:2464139/150‑1 (MQ=255)
                                                                                                           cTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTAc                                           >  1:4548645/1‑150 (MQ=255)
                                                                                                              aTAACGTTGCAACCGTTCTGCCACCGCTTCAACAATCTCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGcc                                        <  2:2093916/150‑1 (MQ=255)
                                                                                                              aTAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGcc                                        <  1:2436121/150‑1 (MQ=255)
                                                                                                               tAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCa                                       <  1:3554217/150‑1 (MQ=255)
                                                                                                                aaCGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCAc                                      <  2:3419090/150‑1 (MQ=255)
                                                                                                                  cGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACgg                                    >  2:101607/1‑150 (MQ=255)
                                                                                                                  cGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACgg                                    <  2:714413/150‑1 (MQ=255)
                                                                                                                  cGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACgg                                    >  2:110921/1‑150 (MQ=255)
                                                                                                                   gTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACggg                                   >  2:3303745/1‑150 (MQ=255)
                                                                                                                   gTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACggg                                   >  1:146909/1‑150 (MQ=255)
                                                                                                                    ttGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGt                                  <  1:3886883/150‑1 (MQ=255)
                                                                                                                     tGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTa                                 <  1:4465850/150‑1 (MQ=255)
                                                                                                                     tGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTa                                 <  1:1292823/150‑1 (MQ=255)
                                                                                                                      gCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAt                                <  1:2460508/150‑1 (MQ=255)
                                                                                                                      gCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAt                                <  1:4370260/150‑1 (MQ=255)
                                                                                                                        aaCCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAttt                              >  1:2399734/1‑150 (MQ=255)
                                                                                                                        aaCCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAttt                              >  2:3825197/1‑150 (MQ=255)
                                                                                                                        aaCCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAtgt                              >  2:1748098/1‑148 (MQ=255)
                                                                                                                            gTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTgcg                          <  2:561983/150‑1 (MQ=255)
                                                                                                                               cTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGcca                       <  2:1984727/150‑1 (MQ=255)
                                                                                                                                  ccACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGccacca                    <  1:3109891/150‑1 (MQ=255)
                                                                                                                                    aCCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGt                  >  1:3488561/1‑150 (MQ=255)
                                                                                                                                      cGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGc                >  1:662323/1‑150 (MQ=255)
                                                                                                                                       gCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCa               >  2:2069663/1‑150 (MQ=255)
                                                                                                                                         ttCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGt             >  2:3226438/1‑150 (MQ=255)
                                                                                                                                         ttCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGt             >  2:4162090/1‑150 (MQ=255)
                                                                                                                                          tCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGta            >  1:1982078/1‑150 (MQ=255)
                                                                                                                                           caacaaTATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGtaa           <  1:3192205/150‑1 (MQ=255)
                                                                                                                                            aacaaTATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGtaat          >  2:3081272/1‑150 (MQ=255)
                                                                                                                                            aacaaTATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGtaat          >  2:2055083/1‑150 (MQ=255)
                                                                                                                                              caaTATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGtaataa        >  2:3112236/1‑150 (MQ=255)
                                                                                                                                              caaTATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGTGTATTTTGCGCCACCAGTGCAGtaataa        <  1:140841/150‑1 (MQ=255)
                                                                                                                                               aaTATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAAc       >  1:414261/1‑150 (MQ=255)
                                                                                                                                                  atCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGa    >  1:770713/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAACCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGC  >  CP009273/2177595‑2177892

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: