Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,453 C→T 41.2% D72D (GAC→GAT talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,4530CT41.2% 62.6 / 130.9 102D72D (GAC→GATtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base T (16/26);  total (48/54)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 104 total aligned reads displayed.

GCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACG  >  CP009273/8304‑8598
                                                                                                                                                     |                                                                                                                                                 
gCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGAc                                                                                                                                                   <  1:787570/150‑1 (MQ=255)
 cAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGAcg                                                                                                                                                  >  1:157184/1‑150 (MQ=255)
     gAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgac                                                                                                                                              >  2:17371/1‑150 (MQ=255)
      aaGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgacc                                                                                                                                             <  2:5709/150‑1 (MQ=255)
      aaGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgacc                                                                                                                                             <  1:1122061/150‑1 (MQ=255)
            tATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACaaa                                                                                                                                       >  2:425519/1‑150 (MQ=255)
                aaCCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTg                                                                                                                                    <  1:1034559/149‑1 (MQ=255)
                 gcctcAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGc                                                                                                                                  >  2:904138/5‑150 (MQ=255)
                  cctcAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGcc                                                                                                                                 <  1:277255/147‑2 (MQ=255)
                  ccGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTgg                                                                                                                                   >  2:363146/1‑148 (MQ=255)
                     cAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTa                                                                                                                              <  2:706123/150‑1 (MQ=255)
                     cAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTa                                                                                                                              <  2:313671/150‑1 (MQ=255)
                       ggATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                                                            <  2:962763/150‑1 (MQ=255)
                       ggATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                                                            >  2:594738/1‑150 (MQ=255)
                        gATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                                                            >  1:967438/1‑149 (MQ=255)
                        gATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                                                            >  2:589407/1‑149 (MQ=255)
                         aTGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTgaaca                                                                                                                          >  2:1109305/1‑148 (MQ=255)
                         aTGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTgaaca                                                                                                                          >  2:738123/1‑148 (MQ=255)
                         aTGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAata                                                                                                                          <  1:979128/150‑1 (MQ=255)
                          tGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAacat                                                                                                                         <  1:610249/150‑1 (MQ=255)
                            ccACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTg                                                                                                                       >  1:965700/1‑150 (MQ=255)
                             caccacTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCgg                                                                                                                      >  2:442801/1‑150 (MQ=255)
                             cacaACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTgg                                                                                                                      >  2:582833/1‑150 (MQ=255)
                              accactaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGt                                                                                                                     <  2:182023/144‑1 (MQ=255)
                              acaACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGt                                                                                                                     >  2:248138/1‑150 (MQ=255)
                                 actaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCtc                                                                                                                  >  2:950846/4‑149 (MQ=255)
                                 actaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCtc                                                                                                                  >  1:443385/4‑149 (MQ=255)
                                      cccTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAAt                                                                                                             <  2:1137026/150‑1 (MQ=255)
                                      cccTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAAt                                                                                                             <  1:796092/150‑1 (MQ=255)
                                       ccTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATc                                                                                                            >  1:479670/1‑150 (MQ=255)
                                           ctctcATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTg                                                                                                         <  2:678858/149‑1 (MQ=255)
                                                aTTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAGCAGCAGAGCAACGATCGCGCGCAGCAGGTCGTGGGCGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAAc                                                                                                     <  2:442265/148‑1 (MQ=255)
                                                  tCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGt                                                                                                 >  2:123565/1‑150 (MQ=255)
                                                   ctcggcgcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTgat                                                                                                 <  2:398679/144‑3 (MQ=255)
                                                   cTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGt                                                                                                 <  2:1191317/149‑1 (MQ=255)
                                                      ggcgcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCg                                                                                             <  2:41865/148‑1 (MQ=255)
                                                      ggcgcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCg                                                                                             <  1:353151/148‑1 (MQ=255)
                                                      ggcgcggctcagATCCCTGAGTAACGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCg                                                                                             <  2:1057072/148‑1 (MQ=255)
                                                      aaCGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCg                                                                                             >  2:303597/1‑150 (MQ=255)
                                                       aCGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCgg                                                                                            <  1:1192639/150‑1 (MQ=255)
                                                          cggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGcc                                                                                         <  1:825723/145‑1 (MQ=255)
                                                             cgcagATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTa                                                                                      >  2:749744/1‑150 (MQ=255)
                                                              gcagATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTAt                                                                                     >  1:376412/1‑150 (MQ=255)
                                                                agATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTAtct                                                                                   <  1:155379/150‑1 (MQ=255)
                                                                 gATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTAtctc                                                                                  >  1:803817/1‑150 (MQ=255)
                                                                   ttCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCaa                                                                                >  1:347520/1‑150 (MQ=255)
                                                                      ctgagtACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCAccg                                                                             <  2:258949/145‑3 (MQ=255)
                                                                       ggAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTg                                                                             >  2:67474/1‑149 (MQ=255)
                                                                       ggAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGa                                                                            <  1:1126570/150‑1 (MQ=255)
                                                                         aaTACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAg                                                                          <  2:564387/150‑1 (MQ=255)
                                                                             ccGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGa                                                                      <  1:1007092/150‑1 (MQ=255)
                                                                             ccGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGa                                                                      >  1:388733/1‑150 (MQ=255)
                                                                                 aaGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATgcg                                                                  >  1:583578/1‑150 (MQ=255)
                                                                                 aaGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATgcg                                                                  >  2:379514/1‑150 (MQ=255)
                                                                                  agctgattgacGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACgcac                                                                 <  2:724435/147‑3 (MQ=255)
                                                                                  agctgattgacGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACgcac                                                                 <  1:976071/147‑3 (MQ=255)
                                                                                     tgattgacGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTc                                                              >  2:675492/1‑150 (MQ=255)
                                                                                       attgatGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCtt                                                            >  2:488927/1‑150 (MQ=255)
                                                                                         tgatgaTGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTc                                                          >  1:341799/1‑150 (MQ=255)
                                                                                         tgacgaCGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGtctctc                                                          >  1:201575/1‑150 (MQ=255)
                                                                                          gatgatGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTcc                                                         >  1:422122/1‑150 (MQ=255)
                                                                                             gatgCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTAt                                                      >  1:1095627/1‑150 (MQ=255)
                                                                                                 cTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGaca                                                  >  2:1190830/1‑150 (MQ=255)
                                                                                                    tCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCg                                               <  2:672247/150‑1 (MQ=255)
                                                                                                            gCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTc                                        >  1:237286/1‑149 (MQ=255)
                                                                                                            gCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCa                                       <  1:409302/150‑1 (MQ=255)
                                                                                                              gAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGa                                                    >  2:1026255/1‑135 (MQ=255)
                                                                                                               aaaCAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGac                                                   <  1:1026255/135‑1 (MQ=255)
                                                                                                                gtagccagAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTAtc                                    <  1:615151/144‑2 (MQ=255)
                                                                                                                gtagccagAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTAtc                                    <  1:277631/144‑2 (MQ=255)
                                                                                                                aaCAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTg                                   <  1:653806/150‑1 (MQ=255)
                                                                                                                 tagccagAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGc                                  >  2:854300/5‑150 (MQ=255)
                                                                                                                 aCAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGc                                  >  2:782393/1‑150 (MQ=255)
                                                                                                                   agcagAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCg                                 >  1:746843/1‑149 (MQ=255)
                                                                                                                     cagAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGa                                >  2:515396/1‑148 (MQ=255)
                                                                                                                     cagAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGaaa                              <  2:214340/150‑1 (MQ=255)
                                                                                                                      agagCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAg                             <  2:253692/150‑1 (MQ=255)
                                                                                                                       gagCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGc                            <  2:186267/150‑1 (MQ=255)
                                                                                                                          cagcGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGcc                           >  2:586931/4‑147 (MQ=255)
                                                                                                                          cagcGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCcaa                         <  1:824920/147‑1 (MQ=255)
                                                                                                                          cAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCaa                          >  1:272108/1‑149 (MQ=255)
                                                                                                                           agcGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCcaag                        <  1:512431/148‑2 (MQ=255)
                                                                                                                           aaCGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCaaaa                        >  1:1019845/1‑150 (MQ=255)
                                                                                                                           aaCGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCaaaa                        <  2:1042122/150‑1 (MQ=255)
                                                                                                                            aCGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAAc                       <  1:17371/150‑1 (MQ=255)
                                                                                                                             cGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACg                      >  2:837782/1‑150 (MQ=255)
                                                                                                                              gATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgc                     <  1:904138/150‑1 (MQ=255)
                                                                                                                              gATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCg                      >  2:1179583/1‑149 (MQ=255)
                                                                                                                               aTCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgca                    <  1:826078/150‑2 (MQ=255)
                                                                                                                                 cgcgcgCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCt                   >  2:824078/1‑149 (MQ=255)
                                                                                                                                 cgcgcgCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgcat                   >  2:758841/1‑147 (MQ=255)
                                                                                                                                 cgcgcgCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgcat                   >  1:56570/1‑147 (MQ=255)
                                                                                                                                  gcgcgcAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGa                 <  2:528165/150‑1 (MQ=255)
                                                                                                                                  gcgcgcAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgcatca                 >  1:667618/1‑146 (MQ=255)
                                                                                                                                   cgcgcAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgcatcat                <  2:326589/150‑6 (MQ=255)
                                                                                                                                      gcagcaGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCaa             <  2:418588/150‑1 (MQ=255)
                                                                                                                                       cagcagATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCaag            <  2:658907/150‑2 (MQ=255)
                                                                                                                                             aTCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTAc      <  1:303597/150‑1 (MQ=255)
                                                                                                                                                 tGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACg  <  2:341799/150‑1 (MQ=255)
                                                                                                                                                 tGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACg  <  2:422122/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACG  >  CP009273/8304‑8598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: