Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,254,793 C→A 38.3% D329Y (GAT→TAT)  tdcE ← pyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,254,7930CA38.3% 375.2 / 292.2 271D329Y (GAT→TAT) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base C (84/83);  new base A (49/55);  total (133/138)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 280 total aligned reads displayed.

AATGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGACA  >  CP009273/3254644‑3254942
                                                                                                                                                     |                                                                                                                                                     
  tGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATaaa                                                                                                                                                     >  2:629041/1‑150 (MQ=255)
     tCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAAtt                                                                                                                                                  <  1:3021637/150‑1 (MQ=255)
      cAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTc                                                                                                                                                 >  2:4080440/1‑150 (MQ=255)
       aggttaggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTcc                                                                                                                                                <  1:4634751/150‑1 (MQ=255)
        ggttaggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCg                                                                                                                                               <  2:3983513/150‑1 (MQ=255)
         gttaggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCgg                                                                                                                                              <  1:2301293/150‑1 (MQ=255)
          ttaggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGgt                                                                                                                                             >  1:2719873/1‑150 (MQ=255)
           taggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGgtg                                                                                                                                            <  1:4237761/150‑1 (MQ=255)
            aggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGgtgt                                                                                                                                           >  2:843834/1‑150 (MQ=255)
               ttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACg                                                                                                                                        >  2:3455108/1‑150 (MQ=255)
                 ccGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCa                                                                                                                                      >  1:4173889/1‑150 (MQ=255)
                  cGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAg                                                                                                                                     >  2:4628891/1‑150 (MQ=255)
                     tGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGACTCAAATTCGGTGTACGCAGGAAg                                                                                                                                 >  2:3866889/1‑150 (MQ=255)
                     tGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGaa                                                                                                                                  <  1:1841513/150‑1 (MQ=255)
                       ccGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGc                                                                                                                                >  1:1784968/1‑150 (MQ=255)
                            cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTAcc                                                                                                                           >  1:1316935/1‑150 (MQ=255)
                            cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATAACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTAcc                                                                                                                           >  1:2376944/1‑150 (MQ=255)
                               aTAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATa                                                                                                                        >  1:3089509/1‑150 (MQ=255)
                               aTAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATa                                                                                                                        >  2:3617305/1‑150 (MQ=255)
                                tAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATAc                                                                                                                       >  2:3715765/1‑150 (MQ=255)
                                tAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATAc                                                                                                                       <  2:4388057/150‑1 (MQ=255)
                                 aGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATAACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACg                                                                                                                      >  2:2356297/1‑150 (MQ=255)
                                    gtgCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGAt                                                                                                                   <  1:2218574/150‑1 (MQ=255)
                                    gtgCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGAt                                                                                                                   <  2:3241209/150‑1 (MQ=255)
                                       cAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCtt                                                                                                                >  2:2852685/1‑150 (MQ=255)
                                        aGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTc                                                                                                               <  2:287895/150‑1 (MQ=255)
                                        aGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTc                                                                                                               >  1:2893390/1‑150 (MQ=255)
                                         gggTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCa                                                                                                              >  2:3490715/1‑150 (MQ=255)
                                         gggTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCa                                                                                                              <  1:3876094/150‑1 (MQ=255)
                                            tgtgCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCatga                                                                                                           >  1:3051853/1‑150 (MQ=255)
                                             gtgCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCatgat                                                                                                          >  2:1761856/1‑150 (MQ=255)
                                                  aaTAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGt                                                                                                     >  1:1504630/1‑150 (MQ=255)
                                                    tAGCTGAAGGAGTTTTATGTCACCATCGTACGACCGACCAGACCCATCCCGCCGATCACTTCCTTCGCCCAGAGTCGGTCGCCGGAGAACAGCGAATCAAATTCCGCTGTACGCAGCAAGCGTACCAAACGCAACAACATCATCAAGTga                                                                                                   >  2:586155/1‑150 (MQ=255)
                                                     aGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgat                                                                                                  <  2:4373865/150‑1 (MQ=255)
                                                       cGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcg                                                                                                >  1:161320/1‑150 (MQ=255)
                                                          aaGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcgatc                                                                                             >  2:381171/1‑150 (MQ=255)
                                                           aGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCa                                                                                            <  1:225142/150‑1 (MQ=255)
                                                           aGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATAACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCa                                                                                            <  1:485814/150‑1 (MQ=255)
                                                             gAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCGTCATGATGAAGTGATCGATCAGt                                                                                          >  2:2719930/1‑150 (MQ=255)
                                                              aGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGtt                                                                                         >  2:1590302/1‑150 (MQ=255)
                                                               gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTc                                                                                        >  2:1722737/1‑150 (MQ=255)
                                                                  tttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATAACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCtg                                                                                     <  1:1543470/150‑1 (MQ=255)
                                                                    tGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCtgtg                                                                                   >  2:1117720/1‑150 (MQ=255)
                                                                      gTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGcc                                                                                 <  2:1884740/150‑1 (MQ=255)
                                                                       tCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCGTCATGATGAAGTGATCGATCAGTTCCTGTGCct                                                                                >  1:1142600/1‑150 (MQ=255)
                                                                       tCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCct                                                                                >  1:3919978/1‑150 (MQ=255)
                                                                        cACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctg                                                                               <  2:789221/150‑1 (MQ=255)
                                                                        cACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctg                                                                               <  2:2018236/150‑1 (MQ=255)
                                                                          ccAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgct                                                                             <  1:1709176/150‑1 (MQ=255)
                                                                           cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgctg                                                                            <  1:4078934/150‑1 (MQ=255)
                                                                           cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgctg                                                                            <  1:1932394/150‑1 (MQ=255)
                                                                           cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgctg                                                                            >  2:2043032/1‑150 (MQ=255)
                                                                             gCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTgctgct                                                                          <  1:941341/150‑1 (MQ=255)
                                                                              cGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTc                                                                         >  2:12954/1‑150 (MQ=255)
                                                                              cGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTc                                                                         <  2:2661605/150‑1 (MQ=255)
                                                                              cGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTc                                                                         <  1:1341790/150‑1 (MQ=255)
                                                                              cGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTc                                                                         <  1:2425164/150‑1 (MQ=255)
                                                                               gTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCa                                                                        >  1:3230784/1‑150 (MQ=255)
                                                                               gTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCa                                                                        <  1:1638098/150‑1 (MQ=255)
                                                                                tACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCAt                                                                       >  1:328157/1‑150 (MQ=255)
                                                                                    aCCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAg                                                                   <  2:2254030/150‑1 (MQ=255)
                                                                                      cGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTa                                                                 >  2:2489781/1‑150 (MQ=255)
                                                                                      cGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTa                                                                 >  2:2950482/1‑150 (MQ=255)
                                                                                       gTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTAc                                                                <  1:4602801/150‑1 (MQ=255)
                                                                                       gTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTAc                                                                <  2:1850137/150‑1 (MQ=255)
                                                                                          cAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGcc                                                             <  1:610661/150‑1 (MQ=255)
                                                                                             ccccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGc                                                          >  1:1076742/1‑150 (MQ=255)
                                                                                             ccccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGc                                                          >  1:1583113/1‑150 (MQ=255)
                                                                                              cccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCt                                                         >  1:2012274/1‑150 (MQ=255)
                                                                                                cATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCttt                                                       <  1:2547440/150‑1 (MQ=255)
                                                                                                 aTCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTa                                                      >  1:3279634/1‑150 (MQ=255)
                                                                                                    ccgccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAg                                                   <  1:3519559/150‑1 (MQ=255)
                                                                                                    ccgccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAg                                                   <  2:4444142/150‑1 (MQ=255)
                                                                                                      gccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTc                                                 >  2:1500568/1‑150 (MQ=255)
                                                                                                      gccgATCACTTCCGTCGCCCAGATAGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGAACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAgac                                                 >  1:3436488/1‑148 (MQ=255)
                                                                                                        cgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTcgc                                               >  1:1457893/1‑150 (MQ=255)
                                                                                                        cgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTcgc                                               >  1:2119014/1‑150 (MQ=255)
                                                                                                                tCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATg                                       <  1:3241485/150‑1 (MQ=255)
                                                                                                                tCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATg                                       <  1:4425251/150‑1 (MQ=255)
                                                                                                                   gTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAg                                    <  1:1083753/150‑1 (MQ=255)
                                                                                                                     cGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGat                                  >  2:3389941/1‑150 (MQ=255)
                                                                                                                       cccAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGatat                                <  1:1227188/150‑1 (MQ=255)
                                                                                                                          aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGa                             <  2:120961/150‑1 (MQ=255)
                                                                                                                            aTTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAgg                           <  1:1978583/150‑1 (MQ=255)
                                                                                                                              tGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGaa                         >  1:3778156/1‑150 (MQ=255)
                                                                                                                                ggTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACg                       >  2:4535370/1‑150 (MQ=255)
                                                                                                                                ggTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATAACGAGGAACg                       <  1:4334973/150‑1 (MQ=255)
                                                                                                                                 gTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATAAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGa                      >  1:2152884/1‑150 (MQ=255)
                                                                                                                                 gTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGa                      <  1:1276872/150‑1 (MQ=255)
                                                                                                                                   cGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATg                    <  1:1313548/150‑1 (MQ=255)
                                                                                                                                      cGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCg                 >  1:4222032/1‑150 (MQ=255)
                                                                                                                                        gagaACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTg               <  2:3071455/150‑1 (MQ=255)
                                                                                                                                        gagaACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGTTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTg               <  2:1068206/150‑1 (MQ=255)
                                                                                                                                         agaACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGACCGATGCCGTGc              >  1:748207/1‑150 (MQ=255)
                                                                                                                                           aaCAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCgg            >  2:3298802/1‑150 (MQ=255)
                                                                                                                                              aGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGccc         <  2:4213421/150‑1 (MQ=255)
                                                                                                                                              aGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGccc         >  2:914289/1‑150 (MQ=255)
                                                                                                                                               gCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCa        <  2:2088943/150‑1 (MQ=255)
                                                                                                                                                 gAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGc      >  2:2643274/1‑150 (MQ=255)
                                                                                                                                                    tCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGAc   <  1:3862834/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGACA  >  CP009273/3254644‑3254942

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: