Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,254,793 C→A 17.9% 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,7930CA17.9% 705.6 / 102.5 234D329Y (GAT→TAT) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base C (92/100);  new base A (30/12);  total (122/112)
Fisher's exact test for biased strand distribution p-value = 6.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 235 total aligned reads displayed.

ATGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGA  >  CP009273/3254645‑3254940
                                                                                                                                                    |                                                                                                                                                   
aTGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCa                                                                                                                                                    >  2:3622328/1‑150 (MQ=255)
 tGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATaaa                                                                                                                                                   >  2:3299736/1‑150 (MQ=255)
  ggTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCaaa                                                                                                                                                  <  2:3755767/150‑1 (MQ=255)
    tCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAAtt                                                                                                                                                >  2:2840683/1‑150 (MQ=255)
    tCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAAtt                                                                                                                                                >  2:2840671/1‑150 (MQ=255)
     cAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTc                                                                                                                                               >  1:1337110/1‑150 (MQ=255)
          taggttCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCCAATAGCGGAAGGAGGTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGCTTGGGTCGCCGGAGAACAGCGAATCAAATTCCGgtg                                                                                                                                          <  2:178954/150‑1 (MQ=255)
               tCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGc                                                                                                                                     <  2:810114/150‑1 (MQ=255)
                ccGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGACTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCa                                                                                                                                    >  2:768045/1‑150 (MQ=255)
                  ggTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAgg                                                                                                                                  >  1:976141/1‑150 (MQ=255)
                   gTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGa                                                                                                                                 >  2:2178884/1‑150 (MQ=255)
                   gTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGa                                                                                                                                 <  2:4087978/150‑1 (MQ=255)
                    tGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGaa                                                                                                                                >  2:943075/1‑150 (MQ=255)
                      ccGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGc                                                                                                                              <  2:589562/150‑1 (MQ=255)
                       cGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCg                                                                                                                             <  1:547943/150‑1 (MQ=255)
                       cGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCg                                                                                                                             >  1:27083/1‑150 (MQ=255)
                        ggCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGt                                                                                                                            >  2:291503/1‑150 (MQ=255)
                        ggCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGt                                                                                                                            >  2:3626672/1‑150 (MQ=255)
                          ccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTAc                                                                                                                          <  1:1683258/150‑1 (MQ=255)
                          ccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTAc                                                                                                                          <  2:3391648/150‑1 (MQ=255)
                           cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTAcc                                                                                                                         >  2:2943040/1‑150 (MQ=255)
                           cccAGAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTAcc                                                                                                                         <  1:1494068/150‑1 (MQ=255)
                             cATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCAt                                                                                                                       >  1:1148958/1‑150 (MQ=255)
                              aTAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATa                                                                                                                      <  2:764501/150‑1 (MQ=255)
                               tAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATAc                                                                                                                     <  2:2697000/150‑1 (MQ=255)
                               tAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATAc                                                                                                                     >  1:3463103/1‑150 (MQ=255)
                                aGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACg                                                                                                                    <  2:2877655/150‑1 (MQ=255)
                                 gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACgg                                                                                                                   <  2:184548/150‑1 (MQ=255)
                                   gtgCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGAt                                                                                                                 >  1:2463938/1‑150 (MQ=255)
                                      cAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCtt                                                                                                              <  2:104825/150‑1 (MQ=255)
                                      cAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCtt                                                                                                              <  2:2039286/150‑1 (MQ=255)
                                         ggTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCat                                                                                                           <  1:344009/150‑1 (MQ=255)
                                         ggTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCat                                                                                                           <  2:2467486/150‑1 (MQ=255)
                                            gtgCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCatgat                                                                                                        <  1:379516/150‑1 (MQ=255)
                                                  aTAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTg                                                                                                  >  1:3923409/1‑150 (MQ=255)
                                                   tAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTga                                                                                                 <  1:1018153/150‑1 (MQ=255)
                                                   tAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTga                                                                                                 >  1:4043338/1‑150 (MQ=255)
                                                    aGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgat                                                                                                >  1:3903810/1‑150 (MQ=255)
                                                     gCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatc                                                                                               <  1:1902336/150‑1 (MQ=255)
                                                       ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcga                                                                                             >  2:2975815/1‑150 (MQ=255)
                                                       ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcga                                                                                             <  1:1576566/150‑1 (MQ=255)
                                                       ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcga                                                                                             >  1:1525806/1‑150 (MQ=255)
                                                         aaGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCTTACGGATCTTCATGATGAAGTgatcgatc                                                                                           <  1:3472189/150‑1 (MQ=255)
                                                         aaGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGATGTgatcgatc                                                                                           <  2:1366129/150‑1 (MQ=255)
                                                         aaGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTgatcgatc                                                                                           <  1:2538767/150‑1 (MQ=255)
                                                           ggAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAg                                                                                         >  1:1884834/1‑150 (MQ=255)
                                                           ggAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAg                                                                                         >  1:3026145/1‑150 (MQ=255)
                                                             aGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGtt                                                                                       >  2:1133759/1‑150 (MQ=255)
                                                              gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTc                                                                                      <  1:3672501/150‑1 (MQ=255)
                                                                 tttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCtg                                                                                   <  1:3746381/150‑1 (MQ=255)
                                                                    ggTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGACGTGATCGATCAGTTCCTGTGc                                                                                >  1:1516647/1‑150 (MQ=255)
                                                                       cACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctg                                                                             >  2:3744755/1‑150 (MQ=255)
                                                                         ccAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgct                                                                           <  1:773752/150‑1 (MQ=255)
                                                                         ccAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgct                                                                           <  2:559841/150‑1 (MQ=255)
                                                                         ccAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCctgct                                                                           <  2:3852735/150‑1 (MQ=255)
                                                                                aCGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCAtt                                                                    <  2:426879/150‑1 (MQ=255)
                                                                                aCGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCAtt                                                                    <  2:3933526/150‑1 (MQ=255)
                                                                                aCGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAAGCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGGTGAAGTGATCGATCAGTGCCTGTGCCTGCTGCTctgt                                                                    >  2:3819564/1‑147 (MQ=255)
                                                                                  gACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGa                                                                  >  1:2128190/1‑150 (MQ=255)
                                                                                     cGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTa                                                               >  1:513981/1‑150 (MQ=255)
                                                                                        ccAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGc                                                            >  1:620285/1‑150 (MQ=255)
                                                                                         cAGCCCCATCCCGACGATCACTTCCGTCGACCAGATTGGGTCGCCGGAGAACAGCGAATCAATTTCAGGTGTACGCAGGAAGCGTACCATACGGATCTTCATTATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGcc                                                           <  1:547743/150‑1 (MQ=255)
                                                                                             cccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCt                                                       >  1:2422272/1‑150 (MQ=255)
                                                                                              ccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCtt                                                      <  2:2418447/150‑1 (MQ=255)
                                                                                              ccATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCtt                                                      <  1:2327249/150‑1 (MQ=255)
                                                                                               cATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCttt                                                     <  2:2640527/150‑1 (MQ=255)
                                                                                               cATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCttt                                                     >  1:2375746/1‑150 (MQ=255)
                                                                                               cATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCttt                                                     <  1:2714714/150‑1 (MQ=255)
                                                                                                 tCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTaa                                                   >  2:3018979/1‑150 (MQ=255)
                                                                                                  cccGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTaaa                                                  >  2:2361831/1‑150 (MQ=255)
                                                                                                    cgccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGt                                                >  1:2231158/1‑150 (MQ=255)
                                                                                                     gccgATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTc                                               <  2:3167085/150‑1 (MQ=255)
                                                                                                        gATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTcgcg                                            >  1:3007941/1‑150 (MQ=255)
                                                                                                         aTCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATAAAATTCCAGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCgggc                                           <  2:2769450/150‑1 (MQ=255)
                                                                                                          tCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCt                                          <  1:1786800/150‑1 (MQ=255)
                                                                                                             cTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCaa                                       <  1:2603647/150‑1 (MQ=255)
                                                                                                              ttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAAt                                      >  1:2206592/1‑150 (MQ=255)
                                                                                                              ttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAAt                                      <  1:3844602/150‑1 (MQ=255)
                                                                                                                     gCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGata                               >  2:1769917/1‑150 (MQ=255)
                                                                                                                     gCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGata                               >  2:1221268/1‑150 (MQ=255)
                                                                                                                      cccAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCCGCTTTAAAGTCGCGCTCAATGTAGatat                              >  2:68177/1‑150 (MQ=255)
                                                                                                                      cccAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGatat                              <  1:2324032/150‑1 (MQ=255)
                                                                                                                        cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCg                            <  1:2577402/150‑1 (MQ=255)
                                                                                                                        cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCg                            >  1:2890717/1‑150 (MQ=255)
                                                                                                                         aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGa                           >  1:2094963/1‑150 (MQ=255)
                                                                                                                          gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAg                          <  2:2579997/150‑1 (MQ=255)
                                                                                                                                  cGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATg                  >  2:174424/1‑150 (MQ=255)
                                                                                                                                       gagaACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTg             >  1:2800770/1‑150 (MQ=255)
                                                                                                                                        agaACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGc            <  1:3414792/150‑1 (MQ=255)
                                                                                                                                         gaACAGCGAATAAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCg           >  1:3222480/1‑150 (MQ=255)
                                                                                                                                          aaCAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCgg          >  2:4026388/1‑150 (MQ=255)
                                                                                                                                            cAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGcc        >  1:3985228/1‑150 (MQ=255)
                                                                                                                                            cAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGcc        <  2:3380362/150‑1 (MQ=255)
                                                                                                                                              gCGACTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGGTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCa      <  2:366174/150‑1 (MQ=255)
                                                                                                                                                gAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGc    >  2:1806834/1‑150 (MQ=255)
                                                                                                                                                gAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGc    >  2:272105/1‑150 (MQ=255)
                                                                                                                                                 aaTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCg   >  2:956715/1‑150 (MQ=255)
                                                                                                                                                 aaTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCg   >  1:3461177/1‑150 (MQ=255)
                                                                                                                                                 aaTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCg   >  2:2510392/1‑150 (MQ=255)
                                                                                                                                                  aTCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGa  >  2:3228127/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTGTGCCTGCTGCTCATTGAGTACGCCAGCTTTAAAGTCGCGCTCAATGTAGATATCGAGGAACGATGCCGTGCGGCCCAGCGA  >  CP009273/3254645‑3254940

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: