Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,254,793 C→A 10.0% 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,7930CA10.0% 1161.3 / 62.0 312D329Y (GAT→TAT) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base C (140/140);  new base A (15/17);  total (155/157)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

AATGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGACA  >  CP009273/3254644‑3254942
                                                                                                                                                     |                                                                                                                                                     
aaTGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATc                                                                                                                                                       >  1:2583322/1‑150 (MQ=255)
     tCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCaaagt                                                                                                                                                  <  2:4218707/150‑3 (MQ=255)
      cAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTc                                                                                                                                                 <  2:3453224/150‑1 (MQ=255)
        ggttaggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCg                                                                                                                                               >  1:2368237/1‑150 (MQ=255)
         gttaggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCgg                                                                                                                                              >  2:781734/1‑150 (MQ=255)
          ttaggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGgt                                                                                                                                             >  2:2559938/1‑150 (MQ=255)
           taggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGgtg                                                                                                                                            >  2:1416023/1‑150 (MQ=255)
                  cGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAg                                                                                                                                     <  2:4127753/150‑1 (MQ=255)
                    gTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGCCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGa                                                                                                                                   <  1:1719253/150‑1 (MQ=255)
                     tGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGaa                                                                                                                                  >  2:32108/1‑150 (MQ=255)
                          gCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTa                                                                                                                             <  2:3148977/150‑1 (MQ=255)
                            cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTAcc                                                                                                                           >  2:1825779/1‑150 (MQ=255)
                             ccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa                                                                                                                          >  2:1046045/1‑150 (MQ=255)
                             ccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCa                                                                                                                          >  2:2652620/1‑150 (MQ=255)
                               aTAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATa                                                                                                                        >  1:51760/1‑150 (MQ=255)
                               aTAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATa                                                                                                                        <  1:3649821/150‑1 (MQ=255)
                               aTAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATa                                                                                                                        <  2:3586508/150‑1 (MQ=255)
                                tAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATAc                                                                                                                       <  1:3892301/150‑1 (MQ=255)
                                   tgtgCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGa                                                                                                                    >  2:2276689/1‑150 (MQ=255)
                                       cAGGGTGTGCAAATAGCGGAAGGAGTTGTTGGTCACCAGCGTACGACCGTCCAGCCCCATCACGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCtt                                                                                                                >  2:898983/1‑150 (MQ=255)
                                        aGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTc                                                                                                               >  2:281248/1‑150 (MQ=255)
                                        aGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTc                                                                                                               >  2:615661/1‑150 (MQ=255)
                                         gggTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCa                                                                                                              <  2:1842249/150‑1 (MQ=255)
                                            tgtgCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCatga                                                                                                           <  2:2254658/150‑1 (MQ=255)
                                             gtgCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCatgat                                                                                                          <  2:1231195/150‑1 (MQ=255)
                                             gtgCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCatgat                                                                                                          <  2:2356466/150‑1 (MQ=255)
                                                  aaTAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGt                                                                                                     <  2:3926601/150‑1 (MQ=255)
                                                  aaTAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGt                                                                                                     >  1:2020520/1‑150 (MQ=255)
                                                     aGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgat                                                                                                  >  1:4291336/1‑150 (MQ=255)
                                                       cGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcg                                                                                                <  1:838414/150‑1 (MQ=255)
                                                         gAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcgat                                                                                              >  1:1184603/1‑150 (MQ=255)
                                                           aGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTACGTCGCCCAGATTGGGTCGCAGGAGAACAGCGAATCAAATTCCGGAGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCa                                                                                            >  2:3501260/1‑150 (MQ=255)
                                                              aGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGtt                                                                                         <  2:3961618/150‑1 (MQ=255)
                                                                tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTcc                                                                                       <  1:1115453/150‑1 (MQ=255)
                                                                 ttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCt                                                                                      <  2:936912/150‑1 (MQ=255)
                                                                 ttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCt                                                                                      <  2:3032098/150‑1 (MQ=255)
                                                                  tttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGAGGACGTGATCGATCAGTTCCtg                                                                                     >  2:3839767/1‑150 (MQ=255)
                                                                   ttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCtgt                                                                                    <  1:412707/150‑1 (MQ=255)
                                                                    tGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCtgtg                                                                                   <  1:2914085/150‑1 (MQ=255)
                                                                         aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgc                                                                              <  1:3065336/150‑1 (MQ=255)
                                                                         aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgc                                                                              >  1:798810/1‑150 (MQ=255)
                                                                          ccAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgct                                                                             <  1:1637437/150‑1 (MQ=255)
                                                                           cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgctg                                                                            >  1:4369359/1‑150 (MQ=255)
                                                                             gCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTgctgct                                                                          <  2:2002080/150‑1 (MQ=255)
                                                                                  cGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTg                                                                     >  2:1163755/1‑150 (MQ=255)
                                                                                    aCCGTCCAGCCCCATCCCGCCGATAACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAg                                                                   <  2:1555558/150‑1 (MQ=255)
                                                                                     ccGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGt                                                                  <  2:1610305/150‑1 (MQ=255)
                                                                                      cGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCTTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTa                                                                 >  2:4466436/1‑150 (MQ=255)
                                                                                       gTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTAc                                                                <  1:3498411/150‑1 (MQ=255)
                                                                                          cAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGcc                                                             >  1:475940/1‑150 (MQ=255)
                                                                                           aGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCa                                                            <  1:1883691/150‑1 (MQ=255)
                                                                                             ccccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGc                                                          >  2:465177/1‑150 (MQ=255)
                                                                                             ccccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGc                                                          >  1:897505/1‑150 (MQ=255)
                                                                                             ccccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGc                                                          >  1:3265523/1‑150 (MQ=255)
                                                                                              cccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCt                                                         >  2:2582745/1‑150 (MQ=255)
                                                                                              cccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCt                                                         <  1:2677120/150‑1 (MQ=255)
                                                                                               ccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCtt                                                        >  2:3713820/1‑150 (MQ=255)
                                                                                                cATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCttt                                                       <  1:4344229/150‑1 (MQ=255)
                                                                                                  tCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTaa                                                     >  2:707199/1‑150 (MQ=255)
                                                                                                  tCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTaa                                                     <  1:3908334/150‑1 (MQ=255)
                                                                                                   cccGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTaaa                                                    <  1:3329082/150‑1 (MQ=255)
                                                                                                     cgccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGt                                                  <  1:3784664/150‑1 (MQ=255)
                                                                                                     cgccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGt                                                  >  2:2148325/1‑150 (MQ=255)
                                                                                                      gctgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTc                                                 <  1:2390739/147‑1 (MQ=255)
                                                                                                      gccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTc                                                 <  1:1942223/150‑1 (MQ=255)
                                                                                                      gccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTc                                                 <  2:3635029/150‑1 (MQ=255)
                                                                                                      gccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTc                                                 <  1:2065823/150‑1 (MQ=255)
                                                                                                      gccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTc                                                 >  2:2968179/1‑150 (MQ=255)
                                                                                                       ccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTcg                                                <  2:1911702/150‑1 (MQ=255)
                                                                                                         gATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTcgcg                                              >  2:2008465/1‑150 (MQ=255)
                                                                                                          aTCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCgcgc                                             >  1:4099579/1‑150 (MQ=255)
                                                                                                             aCTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCa                                          >  2:2904191/1‑150 (MQ=255)
                                                                                                                tCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATg                                       >  1:971304/1‑150 (MQ=255)
                                                                                                                tCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATg                                       >  1:732019/1‑150 (MQ=255)
                                                                                                                tCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATg                                       >  1:158726/1‑150 (MQ=255)
                                                                                                                    tCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGa                                   <  2:4573672/150‑1 (MQ=255)
                                                                                                                      gCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGata                                 <  2:2527050/150‑1 (MQ=255)
                                                                                                                      gCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGata                                 <  1:1465207/150‑1 (MQ=255)
                                                                                                                      gCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGata                                 <  1:1277474/150‑1 (MQ=255)
                                                                                                                        ccAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATc                               <  1:3738050/150‑1 (MQ=255)
                                                                                                                        ccAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATc                               >  1:397695/1‑150 (MQ=255)
                                                                                                                          aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGa                             >  1:287493/1‑150 (MQ=255)
                                                                                                                            aTTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAgg                           <  2:1685135/150‑1 (MQ=255)
                                                                                                                                  tCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGAt                     <  2:4350698/150‑1 (MQ=255)
                                                                                                                                   cGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATg                    <  1:185168/150‑1 (MQ=255)
                                                                                                                                   cGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATg                    >  1:1047560/1‑150 (MQ=255)
                                                                                                                                       ggAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGt                >  1:4539442/1‑150 (MQ=255)
                                                                                                                                           aaCAGCGAATCAAATTTCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCgg            >  1:2796944/1‑150 (MQ=255)
                                                                                                                                           aaCAGCGAATCAAATTTCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCgg            >  1:3188563/1‑150 (MQ=255)
                                                                                                                                             cAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGcc          >  1:2273587/1‑150 (MQ=255)
                                                                                                                                              aGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGccc         >  1:3502730/1‑150 (MQ=255)
                                                                                                                                                 gAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGc      >  2:1418732/1‑150 (MQ=255)
                                                                                                                                                  aaTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCg     <  1:2735413/150‑1 (MQ=255)
                                                                                                                                                  aaTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCg     <  1:1837418/150‑1 (MQ=255)
                                                                                                                                                  aaTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCg     <  1:1670075/150‑1 (MQ=255)
                                                                                                                                                    tCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGAc   <  2:2379943/150‑1 (MQ=255)
                                                                                                                                                    tCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGAc   >  2:665433/1‑150 (MQ=255)
                                                                                                                                                    tCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGAc   >  1:4169040/1‑150 (MQ=255)
                                                                                                                                                     cAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGACa  >  1:2168646/1‑150 (MQ=255)
                                                                                                                                                     cAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGACa  >  1:1087363/1‑150 (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: