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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927349,8400GT94.8% 843.4 / 32.0 261A6A (GCG→GCTfolAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (9/5);  new base T (121/126);  total (130/131)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 271 total aligned reads displayed.

GCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCC  >  CP009273/49693‑49986
                                                                                                                                                   |                                                                                                                                                  
 cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgctgcg                                                                                                                                                 <  2:4307550/150‑1 (MQ=255)
 cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgctgcg                                                                                                                                                 <  2:2075060/150‑1 (MQ=255)
 cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgctgcg                                                                                                                                                 >  1:730160/1‑150 (MQ=255)
  gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGt                                                                                                                                                >  1:4534036/1‑150 (MQ=255)
    gTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTa                                                                                                                                              <  1:149489/150‑1 (MQ=255)
      ccAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGc                                                                                                                                            <  1:4858439/150‑1 (MQ=255)
      ccAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGc                                                                                                                                            >  1:346177/1‑150 (MQ=255)
         gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGGTAGCGGt                                                                                                                                         >  1:1082036/1‑150 (MQ=255)
          ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTa                                                                                                                                        <  2:3296853/150‑1 (MQ=255)
          ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTa                                                                                                                                        >  2:58445/1‑150 (MQ=255)
                 gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATcgcg                                                                                                                                 >  2:2679324/1‑150 (MQ=255)
                         tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGc                                                                                                                         >  2:2891545/1‑150 (MQ=255)
                         tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGc                                                                                                                         >  1:2155506/1‑150 (MQ=255)
                            ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTCGATCGCGTTATCGGCATg                                                                                                                      >  2:2641527/1‑150 (MQ=255)
                             cagcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATgg                                                                                                                     >  1:642200/1‑150 (MQ=255)
                             cagcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATgg                                                                                                                     >  1:4304826/1‑150 (MQ=255)
                                cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaa                                                                                                                  >  1:3110075/1‑150 (MQ=255)
                                cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaa                                                                                                                  <  1:2932576/150‑1 (MQ=255)
                                         aaaTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGc                                                                                                         >  1:1650784/1‑150 (MQ=255)
                                             tttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTg                                                                                                     >  2:4952217/1‑150 (MQ=255)
                                              ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTgg                                                                                                    <  1:3954476/150‑1 (MQ=255)
                                                  tCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcc                                                                                                >  1:4927408/1‑150 (MQ=255)
                                                   cAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcct                                                                                               <  2:1506655/150‑1 (MQ=255)
                                                    aaCCTCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg                                                                                              <  2:371454/150‑1 (MQ=255)
                                                       atcatcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcct                                                                                           <  1:4059159/150‑1 (MQ=255)
                                                        tcatcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcctg                                                                                          <  1:3552320/150‑1 (MQ=255)
                                                               cGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTGTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtct                                                                                   <  2:1187975/150‑1 (MQ=255)
                                                               cGCACCAGTCGACGACGGTTTACGCTTTACGTAGAGTGGCGACAATTTTTTTTATCGTGAAATCTCCATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtct                                                                                   >  2:4212518/1‑150 (MQ=255)
                                                                  aCCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGc                                                                                <  1:3136403/150‑1 (MQ=255)
                                                                    cAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCt                                                                              >  1:4333846/1‑150 (MQ=255)
                                                                     aGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTg                                                                             <  1:3884877/150‑1 (MQ=255)
                                                                      gTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTgg                                                                            >  2:422892/1‑150 (MQ=255)
                                                                      gTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTgg                                                                            <  2:1118303/150‑1 (MQ=255)
                                                                      gTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTgg                                                                            >  1:879170/1‑150 (MQ=255)
                                                                       tCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGt                                                                           >  2:1464029/1‑150 (MQ=255)
                                                                        cgacgaCGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                                                          <  2:873472/150‑1 (MQ=255)
                                                                         gacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGttt                                                                         >  1:3659538/1‑150 (MQ=255)
                                                                          acgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                                                        >  1:3137530/1‑150 (MQ=255)
                                                                           cgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                                       >  2:3503507/1‑150 (MQ=255)
                                                                            gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                                      <  2:4630391/150‑1 (MQ=255)
                                                                            gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                                      >  2:403093/1‑150 (MQ=255)
                                                                              cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACg                                                                    >  1:4345120/1‑150 (MQ=255)
                                                                              cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACg                                                                    <  1:4226341/150‑1 (MQ=255)
                                                                                gTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                                                  <  2:2827651/150‑1 (MQ=255)
                                                                                    aCGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAacac                                                              >  2:796853/1‑150 (MQ=255)
                                                                                    aCGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAacac                                                              <  1:2053021/150‑1 (MQ=255)
                                                                                     cGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACAcc                                                             <  1:1340235/150‑1 (MQ=255)
                                                                                      gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                            <  2:3022427/150‑1 (MQ=255)
                                                                                      gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                            >  1:1009824/1‑150 (MQ=255)
                                                                                      gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                            <  1:4123807/150‑1 (MQ=255)
                                                                                      gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                            <  2:4044565/150‑1 (MQ=255)
                                                                                      gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                                            <  2:165064/150‑1 (MQ=255)
                                                                                         ttACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaa                                                         <  2:29987/150‑1 (MQ=255)
                                                                                         tgaCGTATAGTGGTGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaa                                                         <  2:2632802/148‑1 (MQ=255)
                                                                                          tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaa                                                        >  1:4233652/1‑150 (MQ=255)
                                                                                            cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCATGGTTTAAACGCAACACCTtaaata                                                      >  1:2497943/1‑150 (MQ=255)
                                                                                               ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                                   <  1:816818/150‑1 (MQ=255)
                                                                                                taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAcc                                                  >  2:1916934/1‑150 (MQ=255)
                                                                                                taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAcc                                                  >  2:3815168/1‑150 (MQ=255)
                                                                                                 aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                                 <  2:3938456/150‑1 (MQ=255)
                                                                                                 aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                                 <  2:2410265/150‑1 (MQ=255)
                                                                                                       gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAtt                                           <  1:4410879/150‑1 (MQ=255)
                                                                                                       gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAtt                                           <  2:3282878/150‑1 (MQ=255)
                                                                                                         cAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                         <  1:5148726/150‑1 (MQ=255)
                                                                                                          acttttttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTCTg                                        >  2:341857/3‑150 (MQ=255)
                                                                                                            ttttttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                      <  2:4348419/150‑1 (MQ=255)
                                                                                                              ttttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                                    >  1:2380398/1‑150 (MQ=255)
                                                                                                              atttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                                    <  1:767161/149‑1 (MQ=255)
                                                                                                               tttttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                                   <  2:375385/150‑1 (MQ=255)
                                                                                                                  ttATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                >  1:3048713/1‑150 (MQ=255)
                                                                                                                   tATCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                               <  1:833479/150‑1 (MQ=255)
                                                                                                                    aTCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                              <  1:191823/150‑1 (MQ=255)
                                                                                                                    aTCGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                              <  1:4757416/150‑1 (MQ=255)
                                                                                                                      cGGGAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAcc                            >  2:1392028/1‑150 (MQ=255)
                                                                                                                       gggAAATCTCAATGATCAGTCTGATTGCTGCGTGAGCGGTAGCTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                           <  1:3974939/150‑1 (MQ=255)
                                                                                                                        ggAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                          <  2:4473551/150‑1 (MQ=255)
                                                                                                                         gAAATCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTgg                         >  1:3747292/1‑150 (MQ=255)
                                                                                                                          aaaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg                        >  2:1624664/1‑150 (MQ=255)
                                                                                                                          aaaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg                        <  1:3255739/150‑1 (MQ=255)
                                                                                                                          aaaTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACATggg                        <  2:4363305/150‑1 (MQ=255)
                                                                                                                            aTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGac                      >  1:110015/1‑149 (MQ=255)
                                                                                                                            aTCTCAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaa                      >  1:4547346/1‑150 (MQ=255)
                                                                                                                             tctcAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                     <  1:3271517/150‑1 (MQ=255)
                                                                                                                                cAATGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa                  <  2:215666/150‑1 (MQ=255)
                                                                                                                                 aaTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                 >  1:4072079/1‑150 (MQ=255)
                                                                                                                                 aaTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                 <  1:2551327/150‑1 (MQ=255)
                                                                                                                                 aaTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                 >  2:2058696/1‑150 (MQ=255)
                                                                                                                                  aTGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc                <  2:3110075/150‑1 (MQ=255)
                                                                                                                                   tGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCg               >  1:875746/1‑150 (MQ=255)
                                                                                                                                   tGATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAAGCg               >  2:4769616/1‑150 (MQ=255)
                                                                                                                                    gATCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg              >  2:1737262/1‑150 (MQ=255)
                                                                                                                                      tCAGTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc            >  1:2669142/1‑150 (MQ=255)
                                                                                                                                         gTCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc         <  2:2217304/150‑1 (MQ=255)
                                                                                                                                          tCTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTcc        >  2:917511/1‑150 (MQ=255)
                                                                                                                                           cTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg       >  2:2508349/1‑150 (MQ=255)
                                                                                                                                           cTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg       >  1:4271778/1‑150 (MQ=255)
                                                                                                                                           cTGATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg       >  1:3288952/1‑150 (MQ=255)
                                                                                                                                             gATTGCTGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt     <  2:4534036/150‑1 (MQ=255)
                                                                                                                                             gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt     <  1:2303144/150‑1 (MQ=255)
                                                                                                                                                tgctgcGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGcc  <  2:3635056/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCC  >  CP009273/49693‑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: