Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,249,372 C→A 19.2% A213A (GCG→GCT yhaM ← putative L‑serine dehydratase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,249,3720CA19.2% 870.0 / 148.3 300A213A (GCG→GCTyhaMputative L‑serine dehydratase alpha chain
Reads supporting (aligned to +/- strand):  ref base C (112/130);  new base A (25/33);  total (137/163)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

GCGCCGCCCATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGT  >  CP009273/3249223‑3249521
                                                                                                                                                     |                                                                                                                                                     
gcgcCGCCCATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAc                                                                                                                                                       <  1:3918843/150‑1 (MQ=255)
 cgccgcCCATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACg                                                                                                                                                      >  1:2934207/1‑150 (MQ=255)
  gccgccCATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGc                                                                                                                                                     >  1:1089807/1‑150 (MQ=255)
   ccgccCATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGca                                                                                                                                                    <  2:2501582/150‑1 (MQ=255)
     gccCATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGcaca                                                                                                                                                  <  2:3097235/150‑1 (MQ=255)
      cccATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACaa                                                                                                                                                 >  1:612389/1‑150 (MQ=255)
        cATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAAtt                                                                                                                                               <  1:3460188/150‑1 (MQ=255)
        cATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAAtt                                                                                                                                               <  1:2516509/150‑1 (MQ=255)
          tACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTa                                                                                                                                             <  1:1534992/150‑1 (MQ=255)
          tACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTa                                                                                                                                             <  1:4374576/150‑1 (MQ=255)
          tACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTa                                                                                                                                             <  1:3993723/150‑1 (MQ=255)
             gcgcATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCt                                                                                                                                          <  2:2974706/150‑1 (MQ=255)
               gcATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTc                                                                                                                                        <  1:3015543/150‑1 (MQ=255)
                cATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTcg                                                                                                                                       >  2:3610732/1‑150 (MQ=255)
                 aTCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTcgc                                                                                                                                      >  1:2061708/1‑150 (MQ=255)
                       tGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGACTc                                                                                                                                >  2:2174283/1‑150 (MQ=255)
                        gCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCg                                                                                                                               <  2:354443/150‑1 (MQ=255)
                            cgcTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGaa                                                                                                                           >  1:3798825/1‑150 (MQ=255)
                             gcTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAAt                                                                                                                          >  2:881005/1‑150 (MQ=255)
                                ggTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATaaa                                                                                                                       >  2:2853470/1‑150 (MQ=255)
                                ggTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATaaa                                                                                                                       <  2:3967200/150‑1 (MQ=255)
                                  tACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGc                                                                                                                     >  2:4191461/1‑150 (MQ=255)
                                   aCGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCg                                                                                                                    <  1:964212/150‑1 (MQ=255)
                                      aaTCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGAt                                                                                                                 >  1:2004119/1‑150 (MQ=255)
                                       aTCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATc                                                                                                                <  2:3921419/150‑1 (MQ=255)
                                       aTCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATc                                                                                                                <  2:3454161/150‑1 (MQ=255)
                                          acaATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATcgcc                                                                                                             <  1:4659300/150‑1 (MQ=255)
                                          acaATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATcgcc                                                                                                             <  2:1297422/150‑1 (MQ=255)
                                            aaTGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATcgccgc                                                                                                           >  2:4340708/1‑150 (MQ=255)
                                                 aTGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACg                                                                                                      >  2:3676696/1‑150 (MQ=255)
                                                 aTGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACg                                                                                                      >  2:4715087/1‑150 (MQ=255)
                                                    aaGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACGTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGa                                                                                                   <  2:2320964/150‑1 (MQ=255)
                                                     agagagATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGAc                                                                                                  <  1:4080600/150‑1 (MQ=255)
                                                     agagagATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGAc                                                                                                  <  1:2468718/150‑1 (MQ=255)
                                                     agagagATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGACTCGAGAATAAAGCGGATCGCCGCAAACGGGAc                                                                                                  >  2:4664025/1‑150 (MQ=255)
                                                     agagagATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGAc                                                                                                  >  2:2212373/1‑150 (MQ=255)
                                                      gagagaTCTTTCGCCAGCTAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACt                                                                                                 <  2:2025588/150‑1 (MQ=255)
                                                       agagaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACtt                                                                                                <  1:4460897/150‑1 (MQ=255)
                                                       agagaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACtt                                                                                                <  1:865777/150‑1 (MQ=255)
                                                        gagaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTc                                                                                               <  2:2536235/150‑1 (MQ=255)
                                                         agaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCa                                                                                              >  2:4422602/1‑150 (MQ=255)
                                                         agaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCa                                                                                              >  2:4013327/1‑150 (MQ=255)
                                                          gaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCAt                                                                                             >  2:1704125/1‑150 (MQ=255)
                                                          gaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCAt                                                                                             <  1:1690018/150‑1 (MQ=255)
                                                           aTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCAtt                                                                                            <  1:451299/150‑1 (MQ=255)
                                                                tCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGa                                                                                       >  1:26704/1‑150 (MQ=255)
                                                                   ccAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACt                                                                                    >  1:3858260/1‑150 (MQ=255)
                                                                   ccAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACt                                                                                    <  2:1340417/150‑1 (MQ=255)
                                                                    cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACtt                                                                                   <  2:4157284/150‑1 (MQ=255)
                                                                    cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACtt                                                                                   <  1:1459416/150‑1 (MQ=255)
                                                                    cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACtt                                                                                   <  2:1107687/150‑1 (MQ=255)
                                                                      gCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCa                                                                                 <  2:725280/150‑1 (MQ=255)
                                                                       cAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAg                                                                                <  2:1001665/150‑1 (MQ=255)
                                                                        aaaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAgg                                                                               >  2:1163880/1‑150 (MQ=255)
                                                                         aaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGa                                                                              >  1:2525358/1‑150 (MQ=255)
                                                                             gcgcTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGAtctc                                                                          <  2:2618947/150‑1 (MQ=255)
                                                                             gcgcTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGAtctc                                                                          <  2:232630/150‑1 (MQ=255)
                                                                                cTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTcagc                                                                       <  2:4875144/150‑1 (MQ=255)
                                                                                 tCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACACTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTcagcc                                                                      <  2:3259922/150‑1 (MQ=255)
                                                                                      cTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGcgt                                                                 >  2:4577850/1‑150 (MQ=255)
                                                                                       tGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGcgtc                                                                <  2:1518615/150‑1 (MQ=255)
                                                                                       tGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGcgtc                                                                >  1:2331546/1‑150 (MQ=255)
                                                                                        gTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGcgtcg                                                               >  1:1221823/1‑150 (MQ=255)
                                                                                          ttttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGtt                                                             <  2:624579/150‑1 (MQ=255)
                                                                                           tttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTc                                                            >  2:1164420/1‑150 (MQ=255)
                                                                                              ccAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTgg                                                         >  2:2388089/1‑150 (MQ=255)
                                                                                                aGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGACTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGaa                                                       >  2:3337534/1‑150 (MQ=255)
                                                                                                aGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGaa                                                       >  1:2290161/1‑150 (MQ=255)
                                                                                                  cGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAg                                                     >  2:1949525/1‑150 (MQ=255)
                                                                                                  cGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAg                                                     >  2:750381/1‑150 (MQ=255)
                                                                                                    tCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCa                                                   <  2:2006187/150‑1 (MQ=255)
                                                                                                    tCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGACTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCa                                                   >  2:1062839/1‑150 (MQ=255)
                                                                                                      gcgcCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCAcc                                                 <  1:2176987/150‑1 (MQ=255)
                                                                                                        gcCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGt                                               >  2:1450426/1‑150 (MQ=255)
                                                                                                        gcCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGt                                               <  1:3837805/150‑1 (MQ=255)
                                                                                                         ccAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTc                                              <  2:2765709/150‑1 (MQ=255)
                                                                                                          cAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCa                                             <  1:2881960/150‑1 (MQ=255)
                                                                                                               tGCAGCCCCCACGTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCGTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGga                                        <  2:4494215/150‑1 (MQ=255)
                                                                                                                     ccccACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGActct                                  >  2:3835510/1‑150 (MQ=255)
                                                                                                                         aCTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCt                              >  1:74266/1‑150 (MQ=255)
                                                                                                                          cTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTc                             <  1:1926692/150‑1 (MQ=255)
                                                                                                                          cTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTc                             <  2:2298327/150‑1 (MQ=255)
                                                                                                                           ttACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCg                            >  2:1824604/1‑150 (MQ=255)
                                                                                                                             aCCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGcc                          >  2:3908061/1‑150 (MQ=255)
                                                                                                                              ccGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGccc                         <  2:2661417/150‑1 (MQ=255)
                                                                                                                                  tCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTg                     <  2:1866500/150‑1 (MQ=255)
                                                                                                                                  tCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTg                     <  2:4318833/150‑1 (MQ=255)
                                                                                                                                    aaaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGcc                   >  2:159793/1‑150 (MQ=255)
                                                                                                                                    aaaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGcc                   >  1:3047441/1‑150 (MQ=255)
                                                                                                                                     aaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCca                  <  1:3610732/150‑1 (MQ=255)
                                                                                                                                       ccTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCcaca                <  1:1164420/150‑1 (MQ=255)
                                                                                                                                       ccTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCcaca                >  2:3219601/1‑150 (MQ=255)
                                                                                                                                           ccTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACg            <  1:3240945/150‑1 (MQ=255)
                                                                                                                                           ccTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACg            <  1:492361/150‑1 (MQ=255)
                                                                                                                                            cTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGc           >  2:1409600/1‑150 (MQ=255)
                                                                                                                                              gCGATAAAGCACAATTTAGCGTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCct         <  2:4251936/150‑1 (MQ=255)
                                                                                                                                                gATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCctgc       >  1:2312513/1‑150 (MQ=255)
                                                                                                                                                  tAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCctgctg     >  2:2072643/1‑150 (MQ=255)
                                                                                                                                                  tAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCctgctg     >  1:2612540/1‑150 (MQ=255)
                                                                                                                                                     cGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGt  <  1:1819504/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCCGCCCATACGCGCATCGGATGCCGCGCTGGTACGAATCACAATGGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGT  >  CP009273/3249223‑3249521

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: