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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927349,8400GT99.3% 984.0 / ‑1.8 288A6A (GCG→GCTfolAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (166/119);  total (168/120)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 293 total aligned reads displayed.

AGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGC  >  CP009273/49696‑49985
                                                                                                                                                |                                                                                                                                                 
aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGt                                                                                                                                               >  2:4607410/1‑149 (MQ=255)
 gTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGtt                                                                                                                                              >  1:1486474/1‑149 (MQ=255)
   ccAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAg                                                                                                                                            <  1:81901/149‑1 (MQ=255)
   ccAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAg                                                                                                                                            <  1:1124583/149‑1 (MQ=255)
    cAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGc                                                                                                                                           >  1:3149532/1‑149 (MQ=255)
    cAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCATTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGc                                                                                                                                           <  2:1189010/149‑1 (MQ=255)
      gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGACGACGACGGTTTACGCTTTACGTATAGTGGCGACAATATTTTATATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCgg                                                                                                                                         >  2:2183296/1‑149 (MQ=255)
       ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGt                                                                                                                                        >  1:924891/1‑149 (MQ=255)
        gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTa                                                                                                                                       <  2:4362283/149‑1 (MQ=255)
            gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATc                                                                                                                                   >  1:2378139/1‑149 (MQ=255)
                gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGt                                                                                                                               <  1:3381211/149‑1 (MQ=255)
                    aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGTTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATc                                                                                                                           <  1:495262/149‑1 (MQ=255)
                    aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTTTc                                                                                                                           >  1:3599042/1‑149 (MQ=255)
                     cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCg                                                                                                                          <  2:2395039/149‑1 (MQ=255)
                        gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCa                                                                                                                       >  2:3556508/1‑149 (MQ=255)
                           agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATgg                                                                                                                    >  2:1432571/1‑149 (MQ=255)
                           agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATgg                                                                                                                    >  2:111069/1‑149 (MQ=255)
                            gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                                                                                                                   >  1:3812476/1‑149 (MQ=255)
                            gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                                                                                                                   >  2:917853/1‑149 (MQ=255)
                                      aaaTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATg                                                                                                         >  1:2567522/1‑149 (MQ=255)
                                        aTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGcc                                                                                                       <  2:3691247/149‑1 (MQ=255)
                                              cTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                                                 >  1:2103685/1‑149 (MQ=255)
                                              cTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                                                 >  1:3526385/1‑149 (MQ=255)
                                                 aaCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcct                                                                                              <  1:4375862/149‑1 (MQ=255)
                                                 aaCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcct                                                                                              <  2:2804913/149‑1 (MQ=255)
                                                  aCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGACATCGCCATGCCGTGGAAcctg                                                                                             >  2:4171960/1‑149 (MQ=255)
                                                    atcatcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCCTGCCGTGGAAcctgcc                                                                                           >  2:3257505/1‑149 (MQ=255)
                                                    atcatcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcc                                                                                           >  2:1818644/1‑149 (MQ=255)
                                                     tcatcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcct                                                                                          >  1:1339288/1‑149 (MQ=255)
                                                     tcatcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcct                                                                                          <  1:591082/149‑1 (MQ=255)
                                                       atcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACctgcctgc                                                                                        >  1:1707232/1‑149 (MQ=255)
                                                        tcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctgcc                                                                                       >  1:2441732/1‑149 (MQ=255)
                                                           tCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAt                                                                                    >  1:1074120/1‑149 (MQ=255)
                                                            cGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                                                                   <  2:2906327/149‑1 (MQ=255)
                                                            cGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                                                                   >  1:3433389/1‑149 (MQ=255)
                                                            cGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                                                                   <  1:2988929/149‑1 (MQ=255)
                                                            cGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                                                                   >  2:3503953/1‑149 (MQ=255)
                                                             gCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtct                                                                                  >  2:2553234/1‑149 (MQ=255)
                                                               aCCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCg                                                                                >  1:4556728/1‑149 (MQ=255)
                                                                 cAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGcc                                                                              <  2:4676982/149‑1 (MQ=255)
                                                                 cAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGcc                                                                              <  2:2281359/149‑1 (MQ=255)
                                                                  aGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCt                                                                             <  2:297332/149‑1 (MQ=255)
                                                                  aGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGGTATCGGCATGGAAAACGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCt                                                                             >  2:2651215/1‑149 (MQ=255)
                                                                      gacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                                                         <  1:2244503/149‑1 (MQ=255)
                                                                      gacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                                                         <  2:4484073/149‑1 (MQ=255)
                                                                       acgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGttt                                                                        >  2:2417040/1‑149 (MQ=255)
                                                                         gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                                      >  1:1821500/1‑149 (MQ=255)
                                                                           cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                                                    <  1:3864081/149‑1 (MQ=255)
                                                                           cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                                                    <  1:3709285/149‑1 (MQ=255)
                                                                             gTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGc                                                                  >  2:1483772/1‑149 (MQ=255)
                                                                               ttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                                <  1:3120992/149‑1 (MQ=255)
                                                                                tACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                               >  2:2325744/1‑149 (MQ=255)
                                                                                tACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                               >  1:403400/1‑149 (MQ=255)
                                                                                 aCGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAaca                                                              >  2:2687164/1‑149 (MQ=255)
                                                                                 aCGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAaca                                                              >  1:270310/1‑149 (MQ=255)
                                                                                 aCGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAaca                                                              >  1:2011210/1‑149 (MQ=255)
                                                                                   gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACAcc                                                            <  1:4034335/149‑1 (MQ=255)
                                                                                    cTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                           >  2:16815/1‑149 (MQ=255)
                                                                                          gTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                                                     >  2:1144217/1‑149 (MQ=255)
                                                                                            ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataaa                                                   <  1:3207421/149‑1 (MQ=255)
                                                                                              aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAcc                                                 <  1:4764337/149‑1 (MQ=255)
                                                                                              aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTAGCCTGGTTTAAACGCAACACCTtaaataaaac                                                 >  1:3487983/1‑147 (MQ=255)
                                                                                               gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                                >  1:4269807/1‑149 (MQ=255)
                                                                                                tGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                                               >  1:2177571/1‑149 (MQ=255)
                                                                                                 ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                                              >  2:1674853/1‑149 (MQ=255)
                                                                                                    gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                           >  2:3147599/1‑149 (MQ=255)
                                                                                                    gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                           >  2:3834036/1‑149 (MQ=255)
                                                                                                    gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                           >  1:196824/1‑149 (MQ=255)
                                                                                                    gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAACCGTGAt                                           >  2:2466997/1‑149 (MQ=255)
                                                                                                      cAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                         <  2:4216400/149‑1 (MQ=255)
                                                                                                        aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATg                                       >  2:3084904/1‑149 (MQ=255)
                                                                                                        aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATg                                       >  2:4765237/1‑149 (MQ=255)
                                                                                                        aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTGTAAACGCAACACCTTAAATAAACCCGTGATTATg                                       >  2:3242991/1‑149 (MQ=255)
                                                                                                         ttttttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                      >  1:4398725/1‑149 (MQ=255)
                                                                                                          tttttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                     <  1:4243645/149‑1 (MQ=255)
                                                                                                           ttttgtATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTGAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                    <  2:1250489/149‑1 (MQ=255)
                                                                                                            tttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAATCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                                   >  2:1602148/1‑149 (MQ=255)
                                                                                                               ttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                                <  2:788933/149‑1 (MQ=255)
                                                                                                                 aTCGGGAAATCTCAATGATCAGTCTGATTGCTGAGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                              <  1:4152656/149‑1 (MQ=255)
                                                                                                                     ggAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                          <  1:4695325/149‑1 (MQ=255)
                                                                                                                      gAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                         <  1:2163511/149‑1 (MQ=255)
                                                                                                                       aaaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTgg                        >  2:3402247/1‑149 (MQ=255)
                                                                                                                        aaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg                       >  1:2927217/1‑149 (MQ=255)
                                                                                                                           ctcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                    <  2:3856930/149‑1 (MQ=255)
                                                                                                                            tcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                   <  2:609230/149‑1 (MQ=255)
                                                                                                                            tcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                   >  1:2874524/1‑149 (MQ=255)
                                                                                                                            tcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                   >  2:1379617/1‑149 (MQ=255)
                                                                                                                            tcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                   >  1:3729826/1‑149 (MQ=255)
                                                                                                                               aTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                >  2:2995320/1‑149 (MQ=255)
                                                                                                                               aTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                <  2:1545617/149‑1 (MQ=255)
                                                                                                                               aTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                <  1:3883260/149‑1 (MQ=255)
                                                                                                                                   tCAGTCTGATTGCTGCGTTCGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAGTATGGGCCGCCATACCTGGGAATCAATCGgt            <  1:1734574/149‑1 (MQ=255)
                                                                                                                                   tCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt            <  2:1081270/149‑1 (MQ=255)
                                                                                                                                       tCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc        >  1:1004491/1‑149 (MQ=255)
                                                                                                                                          gATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCgc     >  1:1393174/1‑148 (MQ=255)
                                                                                                                                          gATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt     <  1:1002835/149‑1 (MQ=255)
                                                                                                                                           aTTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt    >  1:3466260/1‑149 (MQ=255)
                                                                                                                                            ttGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTg   >  2:1233566/1‑149 (MQ=255)
                                                                                                                                             tgctgcGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGc  <  2:2148988/149‑1 (MQ=255)
                                                                                                                                             tgctgcGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGc  <  2:3941289/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGC  >  CP009273/49696‑49985

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: