Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,177,743 C→A 10.8% 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,7430CA10.8% 826.6 / 48.2 230D85Y (GAT→TAT) thiDbifunctional hydroxy‑methylpyrimidine kinase/ hydroxy‑phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (104/101);  new base A (10/15);  total (114/116)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 236 total aligned reads displayed.

TGAAACCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGC  >  CP009273/2177594‑2177892
                                                                                                                                                     |                                                                                                                                                     
tGAAACCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                                                                                                                                                       >  2:2766355/1‑150 (MQ=255)
tGAAACCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                                                                                                                                                       <  2:1520172/150‑1 (MQ=255)
  aaaCCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAgg                                                                                                                                                     <  2:2301676/150‑1 (MQ=255)
    aCCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGtt                                                                                                                                                   <  1:674224/150‑1 (MQ=255)
     ccTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACAGCTTCAACAATATCGGttt                                                                                                                                                  >  2:2955522/1‑150 (MQ=255)
      cTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTc                                                                                                                                                 <  1:414666/150‑1 (MQ=255)
       tgtgGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc                                                                                                                                                <  1:3142602/150‑1 (MQ=255)
         tgGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGc                                                                                                                                              >  2:2175293/1‑150 (MQ=255)
          ggCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGcc                                                                                                                                             <  1:1762143/150‑1 (MQ=255)
           gCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCa                                                                                                                                            <  1:928585/150‑1 (MQ=255)
              ataatCGACTGCGCAGCGTAGCAACCGCCGAGGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCa                                                                                                                                         >  1:494581/1‑150 (MQ=255)
               taatCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCAt                                                                                                                                        >  1:354270/1‑150 (MQ=255)
               taatCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCAt                                                                                                                                        >  2:1927179/1‑150 (MQ=255)
                   cGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCg                                                                                                                                    >  1:3852816/1‑150 (MQ=255)
                    gACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGa                                                                                                                                   >  2:4028855/1‑150 (MQ=255)
                    gACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGa                                                                                                                                   <  1:2767201/150‑1 (MQ=255)
                      cTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGAtt                                                                                                                                 <  2:4124588/150‑1 (MQ=255)
                        gcgcAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGAtttt                                                                                                                               <  1:230932/150‑1 (MQ=255)
                         cgcAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTa                                                                                                                              <  1:266882/150‑1 (MQ=255)
                            aGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTg                                                                                                                           <  2:4564186/150‑1 (MQ=255)
                             gCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTgg                                                                                                                          <  1:567798/150‑1 (MQ=255)
                              cGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGt                                                                                                                         <  1:1309117/150‑1 (MQ=255)
                                tAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTAt                                                                                                                       >  2:1261342/1‑150 (MQ=255)
                                    aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGAt                                                                                                                   >  1:2418468/1‑150 (MQ=255)
                                          cgAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGcac                                                                                                             >  1:3663540/1‑150 (MQ=255)
                                          cgAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGcac                                                                                                             >  1:4560936/1‑150 (MQ=255)
                                            aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACAt                                                                                                           <  2:35661/150‑1 (MQ=255)
                                            aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACAt                                                                                                           >  2:4077920/1‑150 (MQ=255)
                                            aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACAt                                                                                                           >  2:439233/1‑150 (MQ=255)
                                                tGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCt                                                                                                       <  2:4610959/150‑1 (MQ=255)
                                                    agcagcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAac                                                                                                   <  1:3007757/150‑1 (MQ=255)
                                                       agcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACAcc                                                                                                <  1:280413/150‑1 (MQ=255)
                                                       agcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACAcc                                                                                                >  2:2145723/1‑150 (MQ=255)
                                                         cGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAAGAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGa                                                                                              <  2:2829585/150‑1 (MQ=255)
                                                           ggTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAAt                                                                                            <  2:3094084/150‑1 (MQ=255)
                                                            gTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATc                                                                                           <  1:750890/150‑1 (MQ=255)
                                                             tCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCg                                                                                          >  2:4332604/1‑150 (MQ=255)
                                                             tCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCg                                                                                          >  2:3211802/1‑150 (MQ=255)
                                                             tCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCg                                                                                          <  1:2045333/150‑1 (MQ=255)
                                                              cgccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGGTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGa                                                                                         >  1:1258588/1‑150 (MQ=255)
                                                               gccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAg                                                                                        >  1:3882955/1‑150 (MQ=255)
                                                               gccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAg                                                                                        >  2:1767320/1‑150 (MQ=255)
                                                                ccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGc                                                                                       <  1:10247/150‑1 (MQ=255)
                                                                 cgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCt                                                                                      <  2:1597174/150‑1 (MQ=255)
                                                                   cTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTgg                                                                                    <  1:1462622/150‑1 (MQ=255)
                                                                   cTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTgg                                                                                    >  2:792413/1‑150 (MQ=255)
                                                                      tttGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGggcg                                                                                 >  2:3736085/1‑150 (MQ=255)
                                                                        tGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGggcggc                                                                               >  1:3108745/1‑150 (MQ=255)
                                                                          ccAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGa                                                                             <  2:425606/150‑1 (MQ=255)
                                                                            aGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACa                                                                           >  1:2596257/1‑150 (MQ=255)
                                                                             gCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACaa                                                                          <  1:858217/150‑1 (MQ=255)
                                                                               aTAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACaaaa                                                                        <  2:3174494/150‑1 (MQ=255)
                                                                               aTAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACaaaa                                                                        >  1:3069377/1‑150 (MQ=255)
                                                                                tAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAAt                                                                       <  1:4634359/150‑1 (MQ=255)
                                                                                tAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAAt                                                                       >  1:3367448/1‑150 (MQ=255)
                                                                                   cGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAg                                                                    <  1:1836772/150‑1 (MQ=255)
                                                                                      tgtCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCt                                                                 <  1:4499465/150‑1 (MQ=255)
                                                                                       gtCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTc                                                                <  1:2175293/150‑1 (MQ=255)
                                                                                        tCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCa                                                               >  2:3364723/1‑150 (MQ=255)
                                                                                         cGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCaa                                                              >  2:391186/1‑150 (MQ=255)
                                                                                         cGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCaa                                                              <  1:1953297/150‑1 (MQ=255)
                                                                                          gAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAAt                                                             >  1:2985767/1‑150 (MQ=255)
                                                                                                cACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATa                                                       <  1:3711100/150‑1 (MQ=255)
                                                                                                  cGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATaca                                                     >  2:3359557/1‑150 (MQ=255)
                                                                                                  cGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATaca                                                     >  2:1092381/1‑150 (MQ=255)
                                                                                                   gTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATacac                                                    <  1:846219/150‑1 (MQ=255)
                                                                                                     tttGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCg                                                  >  2:1758854/1‑150 (MQ=255)
                                                                                                     tttGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCg                                                  >  2:4189539/1‑150 (MQ=255)
                                                                                                      ttGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGa                                                 >  1:1492313/1‑150 (MQ=255)
                                                                                                      ttGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGa                                                 <  2:1225473/150‑1 (MQ=255)
                                                                                                       tGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGAc                                                >  2:796051/1‑150 (MQ=255)
                                                                                                        ggATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACt                                               <  2:1269756/150‑1 (MQ=255)
                                                                                                             tGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACg                                          >  2:1814879/1‑150 (MQ=255)
                                                                                                             tGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACg                                          <  1:1820644/150‑1 (MQ=255)
                                                                                                                tAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCa                                       >  2:2704164/1‑150 (MQ=255)
                                                                                                                 aaCGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCAc                                      >  1:1808301/1‑150 (MQ=255)
                                                                                                                 aaCGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCAc                                      >  1:3962821/1‑150 (MQ=255)
                                                                                                                      tGCAACCGTTCTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTa                                 <  1:4642384/150‑1 (MQ=255)
                                                                                                                       gCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAt                                >  1:3957663/1‑150 (MQ=255)
                                                                                                                         aaCCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAttt                              <  1:430351/150‑1 (MQ=255)
                                                                                                                          aCCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTAtttt                             >  2:2425628/1‑150 (MQ=255)
                                                                                                                            cGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTgc                           >  1:2726399/1‑150 (MQ=255)
                                                                                                                             gTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTgcg                          >  2:2548122/1‑150 (MQ=255)
                                                                                                                              ttCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTgcgc                         >  2:1606763/1‑150 (MQ=255)
                                                                                                                                cTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTCTTTTGCGcca                       >  2:3716861/1‑150 (MQ=255)
                                                                                                                                cTGCCACCGCTTCAACAATATAGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGcca                       <  2:2019696/150‑1 (MQ=255)
                                                                                                                                   ccACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGccacca                    <  1:3481374/150‑1 (MQ=255)
                                                                                                                                   ccACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGccacca                    <  1:3473855/150‑1 (MQ=255)
                                                                                                                                    cACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAg                   <  2:2491859/150‑1 (MQ=255)
                                                                                                                                     aCCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGt                  >  2:530273/1‑150 (MQ=255)
                                                                                                                                          ttCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGt             <  2:1258683/150‑1 (MQ=255)
                                                                                                                                             aacaaTATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGtaat          >  1:2470228/1‑150 (MQ=255)
                                                                                                                                             aacaaTATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGtaat          >  2:3724798/1‑150 (MQ=255)
                                                                                                                                              acaaTATCGGTTTCCGCCAGCATAACTATTTTTGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCTCCACCAGTGCAGtaata         <  1:2955522/150‑1 (MQ=255)
                                                                                                                                                  tatCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTg     >  1:3873519/1‑150 (MQ=255)
                                                                                                                                                  tatCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTg     <  2:1192586/150‑1 (MQ=255)
                                                                                                                                                   atCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGa    >  2:2870220/1‑150 (MQ=255)
                                                                                                                                                    tCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAg   >  1:4030431/1‑150 (MQ=255)
                                                                                                                                                    tCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAg   >  1:1701087/1‑150 (MQ=255)
                                                                                                                                                     cGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGc  >  1:1004818/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAAACCTGTGGCAATAATCGACTGCGCAGCGTAGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCGATTTTAGTGGTATCGATTCGCACATCGCTGAACACCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGC  >  CP009273/2177594‑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: