Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 49,840 G→T 88.2% A6A (GCG→GCT)  folA → dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927349,8400GT88.2% 693.6 / 85.8 251A6A (GCG→GCT) folAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (22/8);  new base T (113/108);  total (135/116)
Fisher's exact test for biased strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

CGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCC  >  CP009273/49694‑49986
                                                                                                                                                  |                                                                                                                                                  
cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgctgcg                                                                                                                                                 <  2:1780231/150‑1 (MQ=255)
cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgctgcg                                                                                                                                                 <  1:3017283/150‑1 (MQ=255)
    tCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAg                                                                                                                                             >  1:247004/1‑150 (MQ=255)
      cAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCg                                                                                                                                           >  2:1019421/1‑150 (MQ=255)
      cAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCg                                                                                                                                           <  1:3597362/150‑1 (MQ=255)
      cAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCg                                                                                                                                           <  1:1787400/150‑1 (MQ=255)
      cAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCg                                                                                                                                           >  2:4056461/1‑150 (MQ=255)
             agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATc                                                                                                                                    >  2:1614198/1‑150 (MQ=255)
                 cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGt                                                                                                                                >  2:1276856/1‑150 (MQ=255)
                 cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGt                                                                                                                                >  1:1209071/1‑150 (MQ=255)
                  gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGtt                                                                                                                               <  1:3272950/150‑1 (MQ=255)
                   tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTa                                                                                                                              >  1:3954445/1‑150 (MQ=255)
                   tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTa                                                                                                                              <  1:2858638/150‑1 (MQ=255)
                    ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTAt                                                                                                                             >  1:409314/1‑150 (MQ=255)
                       cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCgg                                                                                                                          >  2:1633559/1‑150 (MQ=255)
                         cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCa                                                                                                                        >  2:4095424/1‑150 (MQ=255)
                            cagcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATgg                                                                                                                     <  2:2342001/150‑1 (MQ=255)
                               cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaa                                                                                                                  >  1:1048700/1‑150 (MQ=255)
                                agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaaa                                                                                                                 >  1:2387829/1‑150 (MQ=255)
                                 gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAAc                                                                                                                >  2:2010797/1‑150 (MQ=255)
                                   atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGc                                                                                                              >  1:2519424/1‑150 (MQ=255)
                                    tataAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGcc                                                                                                             >  1:909241/1‑150 (MQ=255)
                                      taAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCAt                                                                                                           <  2:3085136/150‑1 (MQ=255)
                                         aaTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGcc                                                                                                        >  2:3282847/1‑150 (MQ=255)
                                         aaTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGcc                                                                                                        <  1:3801201/150‑1 (MQ=255)
                                            tttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTg                                                                                                     >  2:3151961/1‑150 (MQ=255)
                                            tttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTg                                                                                                     <  1:1236070/150‑1 (MQ=255)
                                                cTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTGTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAc                                                                                                 >  2:2832261/1‑150 (MQ=255)
                                                   aaCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg                                                                                              >  2:2842378/1‑150 (MQ=255)
                                                   aaCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg                                                                                              >  2:3952036/1‑150 (MQ=255)
                                                     catcatCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcc                                                                                            <  1:501885/150‑1 (MQ=255)
                                                     catcatCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcc                                                                                            >  1:1732427/1‑150 (MQ=255)
                                                      atcatcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcct                                                                                           >  2:1221260/1‑150 (MQ=255)
                                                        catcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAACTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACctgcctgc                                                                                         >  1:1703066/1‑150 (MQ=255)
                                                        catcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACctgcctgc                                                                                         >  2:2290603/1‑150 (MQ=255)
                                                         ctcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctgcc                                                                                        <  1:817326/149‑1 (MQ=255)
                                                           ccTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGGTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGa                                                                                      >  1:2224509/1‑150 (MQ=255)
                                                             tCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                                                                    <  1:2095027/150‑1 (MQ=255)
                                                             tCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGAGTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                                                                    <  2:131357/150‑1 (MQ=255)
                                                                   cAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCt                                                                              <  1:22804/150‑1 (MQ=255)
                                                                   cAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCt                                                                              <  1:3619335/150‑1 (MQ=255)
                                                                   cAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCt                                                                              <  1:1018109/150‑1 (MQ=255)
                                                                   cAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCt                                                                              <  1:4187143/150‑1 (MQ=255)
                                                                     gTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTgg                                                                            >  2:3047522/1‑150 (MQ=255)
                                                                       cgacgaCGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                                                          <  1:2605375/150‑1 (MQ=255)
                                                                       cgacgaCGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                                                          >  1:18067/1‑150 (MQ=255)
                                                                         acgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                                                        <  2:2723319/150‑1 (MQ=255)
                                                                            acgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                                                     <  2:409204/150‑1 (MQ=255)
                                                                             cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACg                                                                    >  1:4239409/1‑150 (MQ=255)
                                                                              ggTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTGTAAACGc                                                                   >  2:2503969/1‑150 (MQ=255)
                                                                               gTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                                                  <  1:1915387/150‑1 (MQ=255)
                                                                               gTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                                                  <  2:594111/150‑1 (MQ=255)
                                                                                  tACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAAGCAaca                                                               >  1:802087/1‑150 (MQ=255)
                                                                                     gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCTCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                            >  2:3660309/1‑150 (MQ=255)
                                                                                     gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                            <  1:887294/150‑1 (MQ=255)
                                                                                      cTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                                           <  2:2264766/150‑1 (MQ=255)
                                                                                          aCGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                                                       <  2:746734/150‑1 (MQ=255)
                                                                                            gTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                                     >  2:1480468/1‑150 (MQ=255)
                                                                                              ataGTGGCGACATTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                                   <  2:3510673/150‑1 (MQ=255)
                                                                                                aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                                 <  1:2212909/150‑1 (MQ=255)
                                                                                                  tGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                                               <  2:2177006/150‑1 (MQ=255)
                                                                                                   ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                              >  2:1151261/1‑150 (MQ=255)
                                                                                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                            <  2:2453877/150‑1 (MQ=255)
                                                                                                      gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTATCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAtt                                           >  2:1094085/1‑150 (MQ=255)
                                                                                                      gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAtt                                           <  2:2545667/150‑1 (MQ=255)
                                                                                                       aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                          >  1:4031557/1‑150 (MQ=255)
                                                                                                        cAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                         <  1:3115535/150‑1 (MQ=255)
                                                                                                        cAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                         <  2:3619642/150‑1 (MQ=255)
                                                                                                         aaTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATg                                        >  1:270644/1‑150 (MQ=255)
                                                                                                         aaTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATg                                        >  1:3359684/1‑150 (MQ=255)
                                                                                                           ttttttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                      <  1:2701462/150‑1 (MQ=255)
                                                                                                            tttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                     >  1:1693902/1‑150 (MQ=255)
                                                                                                             ttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTGAAATAAACCCGTGATTATGGgcc                                    >  2:1675295/1‑150 (MQ=255)
                                                                                                                 ttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                <  2:3190241/150‑1 (MQ=255)
                                                                                                                  tATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                               <  2:4026538/150‑1 (MQ=255)
                                                                                                                   aTCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                              <  2:254966/150‑1 (MQ=255)
                                                                                                                    tCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAc                             <  1:964626/150‑1 (MQ=255)
                                                                                                                      gggAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                           <  1:3081779/150‑1 (MQ=255)
                                                                                                                       ggAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                          >  2:823074/1‑150 (MQ=255)
                                                                                                                       ggAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                          <  2:1828042/150‑1 (MQ=255)
                                                                                                                       ggAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                          >  2:4008916/1‑150 (MQ=255)
                                                                                                                        gAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTgg                         >  2:1082222/1‑150 (MQ=255)
                                                                                                                        gAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTgg                         >  2:1084001/1‑150 (MQ=255)
                                                                                                                           aTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaa                      >  1:1156546/1‑150 (MQ=255)
                                                                                                                           aTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaa                      >  2:440871/1‑150 (MQ=255)
                                                                                                                             ctcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTCAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                    >  2:540479/1‑150 (MQ=255)
                                                                                                                               cAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa                  >  1:3470602/1‑150 (MQ=255)
                                                                                                                                 aTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc                >  2:2835274/1‑150 (MQ=255)
                                                                                                                                  tGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCg               >  2:151839/1‑150 (MQ=255)
                                                                                                                                   gATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg              <  1:2076109/150‑1 (MQ=255)
                                                                                                                                   gATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg              <  1:1251154/150‑1 (MQ=255)
                                                                                                                                    aTCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt             >  2:2024154/1‑150 (MQ=255)
                                                                                                                                    aTCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt             <  1:3990534/150‑1 (MQ=255)
                                                                                                                                    aTCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt             >  1:2867317/1‑150 (MQ=255)
                                                                                                                                    aTCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAAGCACTCGgt             >  2:2921176/1‑150 (MQ=255)
                                                                                                                                      cAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcg           >  1:2510981/1‑150 (MQ=255)
                                                                                                                                          cTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg       <  2:2976376/150‑1 (MQ=255)
                                                                                                                                            gATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt     <  1:1665892/150‑1 (MQ=255)
                                                                                                                                               tgctgcGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGcc  <  1:1507921/150‑1 (MQ=255)
                                                                                                                                               tgctgcGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGcc  >  2:696707/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCC  >  CP009273/49694‑49986

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: ‑