Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,601,793 C→A 6.5% S64Y (TCT→TAT)  tam → trans‑aconitate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,601,7930CA 6.5% 1137.3 / 32.8 297S64Y (TCT→TAT) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (132/145);  new base A (12/8);  total (144/153)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

TGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTC  >  CP009273/1601644‑1601942
                                                                                                                                                     |                                                                                                                                                     
 gAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtct                                                                                                                                                      >  1:2739467/1‑150 (MQ=255)
 gAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtct                                                                                                                                                      >  1:179577/1‑150 (MQ=255)
   aCGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTcc                                                                                                                                                    <  1:2295221/150‑1 (MQ=255)
    cGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCg                                                                                                                                                   <  1:922957/150‑1 (MQ=255)
    cGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCg                                                                                                                                                   <  1:4193550/150‑1 (MQ=255)
     gATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCgg                                                                                                                                                  >  2:212661/1‑150 (MQ=255)
     gATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCgg                                                                                                                                                  >  1:539285/1‑150 (MQ=255)
     gATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCgg                                                                                                                                                  >  2:4172850/1‑150 (MQ=255)
     gATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCgg                                                                                                                                                  >  2:1979257/1‑150 (MQ=255)
         gcgTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCgat                                                                                                                                              <  2:722809/150‑1 (MQ=255)
              cggcggTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTg                                                                                                                                         <  2:3516911/150‑1 (MQ=255)
                gcggTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCt                                                                                                                                       <  2:4075412/150‑1 (MQ=255)
                   gTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGaa                                                                                                                                    <  1:374219/150‑1 (MQ=255)
                     ggAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAgc                                                                                                                                  <  2:2527140/150‑1 (MQ=255)
                      gAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAgcg                                                                                                                                 <  1:4626461/150‑1 (MQ=255)
                      gAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAgcg                                                                                                                                 <  1:1981582/150‑1 (MQ=255)
                       aGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAgcgc                                                                                                                                <  1:3141287/150‑1 (MQ=255)
                       aGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAgcgc                                                                                                                                <  1:3141267/150‑1 (MQ=255)
                            ctTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGt                                                                                                                           >  1:2911216/1‑150 (MQ=255)
                            ctTGCCAGAGTGCCGCTGGAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTg                                                                                                                          >  2:3647961/1‑150 (MQ=255)
                             ttGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTg                                                                                                                          <  2:3763737/150‑1 (MQ=255)
                                 cAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGAAGATATCTGAAGCGCGCAGTGCttt                                                                                                                      >  1:1401539/1‑150 (MQ=255)
                                  agagTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTg                                                                                                                     >  2:675141/1‑150 (MQ=255)
                                   gagTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGc                                                                                                                    >  1:984451/1‑150 (MQ=255)
                                   gagTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGc                                                                                                                    >  1:3855214/1‑150 (MQ=255)
                                     gTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCa                                                                                                                  <  2:2620075/150‑1 (MQ=255)
                                          gCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTg                                                                                                             <  2:2744808/150‑1 (MQ=255)
                                            tGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGcc                                                                                                           >  1:2815883/1‑150 (MQ=255)
                                            tGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGcc                                                                                                           >  1:4046467/1‑150 (MQ=255)
                                                 aaTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTc                                                                                                      <  1:2351381/150‑1 (MQ=255)
                                                 aaTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTc                                                                                                      >  1:551016/1‑150 (MQ=255)
                                                 aaTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTc                                                                                                      <  1:2087079/150‑1 (MQ=255)
                                                  aTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCg                                                                                                     >  2:2200331/1‑150 (MQ=255)
                                                       gAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCCGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGaa                                                                                                >  2:2825200/1‑150 (MQ=255)
                                                        aaTATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAg                                                                                               <  1:1549716/150‑1 (MQ=255)
                                                           atGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAg                                                                                            >  1:669376/1‑150 (MQ=255)
                                                            tGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGa                                                                                           <  2:4257357/150‑1 (MQ=255)
                                                             gTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGAc                                                                                          >  1:3755665/1‑150 (MQ=255)
                                                             gTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGACAGGATAACAGGCATCGACTCGTCTCCGGAGATGATTTCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGAc                                                                                          >  1:1700463/1‑150 (MQ=255)
                                                              ttGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACa                                                                                         <  1:675161/150‑1 (MQ=255)
                                                                 ccGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATcc                                                                                      <  2:312660/150‑1 (MQ=255)
                                                                  cGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCg                                                                                     <  2:414856/150‑1 (MQ=255)
                                                                    aTCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCa                                                                                   >  1:2731301/1‑150 (MQ=255)
                                                                    aTCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCa                                                                                   >  2:2907536/1‑150 (MQ=255)
                                                                     tCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCaa                                                                                  >  2:3669022/1‑150 (MQ=255)
                                                                         ggCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTgg                                                                              >  2:3965491/1‑150 (MQ=255)
                                                                          gCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGc                                                                             <  2:3846530/150‑1 (MQ=255)
                                                                             gtgGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAAc                                                                          <  2:1631177/150‑1 (MQ=255)
                                                                              tgGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAAcc                                                                         <  2:1874794/150‑1 (MQ=255)
                                                                                  ccAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTa                                                                     <  2:4004941/150‑1 (MQ=255)
                                                                                   cAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTAc                                                                    >  2:3847506/1‑150 (MQ=255)
                                                                                    aGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                   <  1:2162291/150‑1 (MQ=255)
                                                                                    aGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                   >  1:707358/1‑150 (MQ=255)
                                                                                        aaCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCa                                                               <  2:4384043/150‑1 (MQ=255)
                                                                                         aCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCAc                                                              >  2:3197152/1‑150 (MQ=255)
                                                                                         aCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCAc                                                              >  1:3545117/1‑150 (MQ=255)
                                                                                           aGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTc                                                            >  1:1230121/1‑150 (MQ=255)
                                                                                            gCACCGCCCTTCTACACCAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCg                                                           <  2:2279992/150‑1 (MQ=255)
                                                                                              aCCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGAt                                                         <  1:2518745/150‑1 (MQ=255)
                                                                                                 gCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTg                                                      <  2:1057177/150‑1 (MQ=255)
                                                                                                    cTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtt                                                   <  2:1580783/150‑1 (MQ=255)
                                                                                                      tCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                 <  1:244730/150‑1 (MQ=255)
                                                                                                       cTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttttt                                                <  1:3858954/150‑1 (MQ=255)
                                                                                                       cTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttttt                                                <  2:1390517/150‑1 (MQ=255)
                                                                                                          caacaaCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCt                                             >  1:342777/1‑150 (MQ=255)
                                                                                                           acccacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTa                                            <  1:3034397/146‑1 (MQ=255)
                                                                                                           aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTa                                            >  2:1522882/1‑150 (MQ=255)
                                                                                                             caacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAAt                                          >  1:2587743/1‑150 (MQ=255)
                                                                                                              aacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATg                                         >  2:2556982/1‑150 (MQ=255)
                                                                                                                cGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                       >  1:4270743/1‑150 (MQ=255)
                                                                                                                 gTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCt                                      >  1:260331/1‑150 (MQ=255)
                                                                                                                  ttGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTc                                     <  2:2668539/150‑1 (MQ=255)
                                                                                                                   tGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCa                                    >  1:2428561/1‑150 (MQ=255)
                                                                                                                    ggCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCAc                                   <  1:4505565/150‑1 (MQ=255)
                                                                                                                     gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACt                                  <  2:2890106/150‑1 (MQ=255)
                                                                                                                     gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACt                                  <  1:843911/150‑1 (MQ=255)
                                                                                                                     gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACt                                  <  2:4314710/150‑1 (MQ=255)
                                                                                                                           ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                            >  1:1041203/1‑150 (MQ=255)
                                                                                                                           ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                            <  2:2347421/150‑1 (MQ=255)
                                                                                                                            gCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATgg                           <  2:179577/150‑1 (MQ=255)
                                                                                                                             ccAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGc                          >  1:2835324/1‑150 (MQ=255)
                                                                                                                             ccAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGc                          <  2:4140388/150‑1 (MQ=255)
                                                                                                                              cAGGATAGCAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                         >  2:4245909/1‑150 (MQ=255)
                                                                                                                              cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                         <  1:3379646/150‑1 (MQ=255)
                                                                                                                              cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                         <  1:212661/150‑1 (MQ=255)
                                                                                                                              cAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                         >  2:53664/1‑150 (MQ=255)
                                                                                                                               aGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGACGCGCGCAGTGCTTTGCCAGACAGCCAGTTCGTGGAAGCAGACATCCGCAACAGGCAACCGGTACAGGCACTCGATCTGATTATTGCTAATGCATCACTGCCATGGCTg                        >  2:153671/1‑150 (MQ=255)
                                                                                                                               aGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTg                        >  2:3427542/1‑150 (MQ=255)
                                                                                                                               aGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTg                        <  2:4559835/150‑1 (MQ=255)
                                                                                                                                 gATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGcc                      <  2:3860727/150‑1 (MQ=255)
                                                                                                                                 gATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCACTGGCTGcc                      >  1:151028/1‑150 (MQ=255)
                                                                                                                                     aCAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAc                  >  1:565153/1‑150 (MQ=255)
                                                                                                                                        ggCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATCGCTGCCCGACCAc               >  2:3271951/1‑150 (MQ=255)
                                                                                                                                          cATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTCATGCCTCACTGCAATGGCTGCCCGACCACTa             >  2:339215/1‑150 (MQ=255)
                                                                                                                                           aTCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTAc            <  1:1169260/150‑1 (MQ=255)
                                                                                                                                               aCTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAt        >  2:3364868/1‑150 (MQ=255)
                                                                                                                                                 tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTg      <  1:3073164/150‑1 (MQ=255)
                                                                                                                                                 tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTg      <  2:4596805/150‑1 (MQ=255)
                                                                                                                                                 tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTg      <  2:1454292/150‑1 (MQ=255)
                                                                                                                                                   gTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGtt    >  1:2397469/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTC  >  CP009273/1601644‑1601942

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: