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,8400GT97.7% 729.8 / 7.5 216A6A (GCG→GCTfolAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (122/89);  total (125/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 218 total aligned reads displayed.

CGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCC  >  CP009273/49694‑49986
                                                                                                                                                  |                                                                                                                                                  
 gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGt                                                                                                                                                >  2:1269021/1‑150 (MQ=255)
  aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGtt                                                                                                                                               >  1:2961387/1‑150 (MQ=255)
   gTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTa                                                                                                                                              <  2:2061177/150‑1 (MQ=255)
     ccAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGc                                                                                                                                            <  1:4214193/150‑1 (MQ=255)
       aGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCgg                                                                                                                                          >  1:507946/1‑150 (MQ=255)
       aGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCgg                                                                                                                                          <  1:2418871/150‑1 (MQ=255)
        gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGt                                                                                                                                         >  1:2604507/1‑150 (MQ=255)
        gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGt                                                                                                                                         >  1:3810350/1‑150 (MQ=255)
          gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAg                                                                                                                                       >  1:194701/1‑150 (MQ=255)
          gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAg                                                                                                                                       >  1:217328/1‑150 (MQ=255)
            gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGAt                                                                                                                                     >  1:3265763/1‑150 (MQ=255)
             agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATc                                                                                                                                    >  1:1623786/1‑150 (MQ=255)
             agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATc                                                                                                                                    <  2:2065097/150‑1 (MQ=255)
                 cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGt                                                                                                                                >  2:2771933/1‑150 (MQ=255)
                      aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGATTACGCTTTACGTATAGTGGCGACAATTTTATTTATCGGGAAAACTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGtgagcg                                                                                                                           >  2:1226970/1‑145 (MQ=255)
                       cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCgg                                                                                                                          <  1:953631/150‑1 (MQ=255)
                         cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCa                                                                                                                        >  2:3357019/1‑150 (MQ=255)
                          gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCAt                                                                                                                       <  2:4328439/150‑1 (MQ=255)
                                agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaaa                                                                                                                 >  2:3062775/1‑150 (MQ=255)
                                agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaaa                                                                                                                 >  2:2317081/1‑150 (MQ=255)
                                   atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCACTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGc                                                                                                              >  1:1990242/1‑150 (MQ=255)
                                   atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGc                                                                                                              <  2:1536518/150‑1 (MQ=255)
                                         aaTTTTCCTCAACATAATCCTCGCACCAGTCGACGAAGGTTTACGCTTTACGTATAGTGGAGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGAGTTATCGGCATGGAAAACGCCATGcc                                                                                                        <  2:1390068/150‑1 (MQ=255)
                                          aTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCg                                                                                                       >  1:2594329/1‑150 (MQ=255)
                                               ccTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                                                  <  1:4251950/150‑1 (MQ=255)
                                               ccTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                                                  <  1:1082646/150‑1 (MQ=255)
                                                cTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAc                                                                                                 <  2:1059426/150‑1 (MQ=255)
                                                   aaCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg                                                                                              >  1:446834/1‑150 (MQ=255)
                                                    aCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgc                                                                                             >  1:3915412/1‑150 (MQ=255)
                                                         atcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctgcc                                                                                        >  2:937153/1‑150 (MQ=255)
                                                           ccTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGa                                                                                      <  1:4417802/150‑1 (MQ=255)
                                                           ccTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGa                                                                                      <  1:2157832/150‑1 (MQ=255)
                                                             tCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                                                                    >  2:334098/1‑150 (MQ=255)
                                                               gCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtctc                                                                                  >  2:4456088/1‑150 (MQ=255)
                                                                cACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCg                                                                                 >  1:583232/1‑150 (MQ=255)
                                                                    aGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTg                                                                             >  1:3864255/1‑150 (MQ=255)
                                                                    aGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTg                                                                             >  1:1421903/1‑150 (MQ=255)
                                                                     gTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTgg                                                                            >  1:799484/1‑150 (MQ=255)
                                                                      tCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGt                                                                           >  2:916407/1‑150 (MQ=255)
                                                                       cgacgaCGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                                                          >  2:3668839/1‑150 (MQ=255)
                                                                        gacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGttt                                                                         <  1:670393/150‑1 (MQ=255)
                                                                         acgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                                                        <  1:189561/150‑1 (MQ=255)
                                                                         acgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                                                        <  1:3303557/150‑1 (MQ=255)
                                                                            acgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                                                     >  1:1874074/1‑150 (MQ=255)
                                                                            acgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                                                     >  1:1844400/1‑150 (MQ=255)
                                                                              ggTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGTAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGc                                                                   >  2:2521686/1‑150 (MQ=255)
                                                                              ggTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGc                                                                   >  2:4444900/1‑150 (MQ=255)
                                                                                   aCGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAacac                                                              >  2:104059/1‑150 (MQ=255)
                                                                                        ttACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaa                                                         >  1:291094/1‑150 (MQ=255)
                                                                                         tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaa                                                        >  1:2570888/1‑150 (MQ=255)
                                                                                            gTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                                     >  2:4126123/1‑150 (MQ=255)
                                                                                            gTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                                     >  1:2930604/1‑150 (MQ=255)
                                                                                              ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                                   >  1:925212/1‑150 (MQ=255)
                                                                                              ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                                   >  1:3215232/1‑150 (MQ=255)
                                                                                              ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                                   >  1:3055235/1‑150 (MQ=255)
                                                                                                aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                                 <  1:2069287/150‑1 (MQ=255)
                                                                                                 gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                                                <  1:3854686/150‑1 (MQ=255)
                                                                                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                            <  2:264812/150‑1 (MQ=255)
                                                                                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                            <  1:4299117/150‑1 (MQ=255)
                                                                                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                            <  1:1825728/150‑1 (MQ=255)
                                                                                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                            <  2:1852311/150‑1 (MQ=255)
                                                                                                      gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAtt                                           >  2:2727858/1‑150 (MQ=255)
                                                                                                       aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                          >  2:1074895/1‑150 (MQ=255)
                                                                                                       aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                          <  1:3562793/150‑1 (MQ=255)
                                                                                                          aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                       >  1:3339848/1‑150 (MQ=255)
                                                                                                                tttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                                 >  2:578239/1‑150 (MQ=255)
                                                                                                                  tATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                               <  1:516799/150‑1 (MQ=255)
                                                                                                                   aTCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                              >  1:2747960/1‑150 (MQ=255)
                                                                                                                    tCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAc                             >  1:1676234/1‑150 (MQ=255)
                                                                                                                    tCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAc                             >  1:1498921/1‑150 (MQ=255)
                                                                                                                     cGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAcc                            <  1:3161411/150‑1 (MQ=255)
                                                                                                                      gggAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                           <  1:3737647/150‑1 (MQ=255)
                                                                                                                         aaaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg                        <  2:3737484/150‑1 (MQ=255)
                                                                                                                         aaaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg                        >  2:3088993/1‑150 (MQ=255)
                                                                                                                          aaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                       <  1:2474613/150‑1 (MQ=255)
                                                                                                                          aaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                       >  2:4350329/1‑150 (MQ=255)
                                                                                                                           aTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaa                      <  1:1651529/150‑1 (MQ=255)
                                                                                                                            tctcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                     <  1:2757556/150‑1 (MQ=255)
                                                                                                                            tctcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                     >  1:905790/1‑150 (MQ=255)
                                                                                                                             ctcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                    >  1:3580546/1‑150 (MQ=255)
                                                                                                                              tcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                   <  1:1368377/150‑1 (MQ=255)
                                                                                                                                aaTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                 >  2:726565/1‑150 (MQ=255)
                                                                                                                                aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                 >  2:1632948/1‑150 (MQ=255)
                                                                                                                                 aTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCTCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc                <  2:1561945/150‑1 (MQ=255)
                                                                                                                                 aTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGactcaatc                >  2:2226905/1‑150 (MQ=255)
                                                                                                                                  tGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCg               <  1:2231771/150‑1 (MQ=255)
                                                                                                                                     tCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc            <  1:3665014/150‑1 (MQ=255)
                                                                                                                                      cattCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcg           >  2:493305/4‑150 (MQ=255)
                                                                                                                                      cAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcg           >  1:4497997/1‑150 (MQ=255)
                                                                                                                                        gTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc         >  2:2409789/1‑150 (MQ=255)
                                                                                                                                        gTCTGATTGCTGCGTTAGCGGTAGATCGAGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc         <  2:3481699/150‑1 (MQ=255)
                                                                                                                                         tCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTcc        <  1:3420512/150‑1 (MQ=255)
                                                                                                                                           tGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt      >  1:1618733/1‑150 (MQ=255)
                                                                                                                                           tGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt      >  2:560022/1‑150 (MQ=255)
                                                                                                                                           tGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt      >  2:163367/1‑150 (MQ=255)
                                                                                                                                           tGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt      >  2:1735193/1‑150 (MQ=255)
                                                                                                                                            gATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt     <  2:1491665/150‑1 (MQ=255)
                                                                                                                                             aTTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTg    <  1:4549199/150‑1 (MQ=255)
                                                                                                                                             aTTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTg    <  2:1569654/150‑1 (MQ=255)
                                                                                                                                               tgctgcGTTAGCGGTAGATCGAGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGcc  >  1:4076985/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: