Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,784,952 A→G 19.0% F184L (TTT→CTT)  ydiU ← UPF0061 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,784,9520AG19.0% 705.4 / 162.2 296F184L (TTT→CTT) ydiUUPF0061 family protein
Reads supporting (aligned to +/- strand):  ref base A (110/129);  new base G (21/35);  total (131/165)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 309 total aligned reads displayed.

CATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCAC  >  CP009273/1784804‑1785101
                                                                                                                                                    |                                                                                                                                                     
 aTTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAAt                                                                                                                                                     <  2:3004587/150‑1 (MQ=255)
  ttGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATg                                                                                                                                                    <  1:4697142/150‑1 (MQ=255)
   tGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGt                                                                                                                                                   <  2:4200886/150‑1 (MQ=255)
    gggCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGtt                                                                                                                                                  <  1:4441468/150‑1 (MQ=255)
       cAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGAAGAACCTTTTCCGACTCGCGGAGGTAGTAAAAATGTTCGa                                                                                                                                               >  1:2412745/1‑150 (MQ=255)
        aaTTAACGACGCGGTACGTGCGACTACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTCACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGaa                                                                                                                                              >  2:1345018/1‑150 (MQ=255)
         aTTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGaaa                                                                                                                                             >  2:842656/1‑150 (MQ=255)
           tAACGACGCGGTACGTGCGACTACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATg                                                                                                                                           >  2:948865/1‑150 (MQ=255)
           tAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTCAAAATGTTCGAAATg                                                                                                                                           >  2:4459089/1‑150 (MQ=255)
           tAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATg                                                                                                                                           <  2:5006015/150‑1 (MQ=255)
            atcgacgCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGa                                                                                                                                          <  2:5225862/148‑1 (MQ=255)
            aaCGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACATTTTACGACTCGCGGCGGTAGTAAAGATGTTAGAAATGa                                                                                                                                          >  1:4399664/1‑150 (MQ=255)
             acgacgCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGAc                                                                                                                                         >  2:1920637/1‑150 (MQ=255)
                  gcgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAg                                                                                                                                    <  2:3902933/150‑1 (MQ=255)
                   cgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAgc                                                                                                                                   >  1:2581425/1‑150 (MQ=255)
                    ggTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAgcg                                                                                                                                  >  2:3637743/1‑150 (MQ=255)
                    ggTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAgcg                                                                                                                                  <  2:3550611/150‑1 (MQ=255)
                      tACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCa                                                                                                                                <  1:2317320/150‑1 (MQ=255)
                      tACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCa                                                                                                                                >  1:3563332/1‑150 (MQ=255)
                       aCGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAg                                                                                                                               <  1:3541933/150‑1 (MQ=255)
                         gTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGAt                                                                                                                             <  2:3147970/150‑1 (MQ=255)
                           gCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGa                                                                                                                           <  2:2319596/150‑1 (MQ=255)
                             gACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACt                                                                                                                         >  1:2657668/1‑150 (MQ=255)
                              acaacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACtt                                                                                                                        <  1:4362406/150‑1 (MQ=255)
                               caacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTCGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTg                                                                                                                       >  2:2671960/1‑150 (MQ=255)
                               caacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTg                                                                                                                       >  1:1966156/1‑150 (MQ=255)
                                aacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTgg                                                                                                                      >  1:3155477/1‑150 (MQ=255)
                                  caTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTg                                                                                                                    >  1:2882473/1‑150 (MQ=255)
                                       cTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCcaca                                                                                                               <  1:3073385/150‑1 (MQ=255)
                                         aaaCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACg                                                                                                             >  1:482800/1‑150 (MQ=255)
                                         aaaCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACg                                                                                                             <  2:2847519/150‑1 (MQ=255)
                                         aaaCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACg                                                                                                             <  1:919503/150‑1 (MQ=255)
                                              agagACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCa                                                                                                        <  1:4229052/150‑1 (MQ=255)
                                                 gaCGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTCACTGACGAACCTATTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCa                                                                                                     >  2:4953187/1‑150 (MQ=255)
                                                   cGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATc                                                                                                   <  1:4033937/150‑1 (MQ=255)
                                                      tATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATcgcg                                                                                                >  2:3077467/1‑150 (MQ=255)
                                                      tATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATcgcg                                                                                                >  2:1189763/1‑150 (MQ=255)
                                                         ttGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                             <  1:1985605/150‑1 (MQ=255)
                                                         ttGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                             <  1:1020712/150‑1 (MQ=255)
                                                         ttGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                             <  1:1472268/150‑1 (MQ=255)
                                                          tGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTg                                                                                            <  1:3676364/150‑1 (MQ=255)
                                                          tGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTg                                                                                            <  1:139057/150‑1 (MQ=255)
                                                             ccTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCt                                                                                         <  2:4438526/150‑1 (MQ=255)
                                                             ccTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCt                                                                                         >  1:1867452/1‑150 (MQ=255)
                                                              cTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTc                                                                                        <  2:4158179/150‑1 (MQ=255)
                                                              cTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTc                                                                                        >  2:548936/1‑150 (MQ=255)
                                                               tcatcaTCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTcc                                                                                       <  2:2423863/150‑1 (MQ=255)
                                                               tcatcaTCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTcc                                                                                       <  1:1646912/150‑1 (MQ=255)
                                                               tcatcaTCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTcc                                                                                       <  1:3851433/150‑1 (MQ=255)
                                                                catcatCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccg                                                                                      <  1:3775993/150‑1 (MQ=255)
                                                                      cTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGttt                                                                                >  2:4295275/1‑150 (MQ=255)
                                                                         cAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                                             <  1:3231470/150‑1 (MQ=255)
                                                                         cAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                                             <  2:234137/150‑1 (MQ=255)
                                                                         cAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                                             >  1:773559/1‑150 (MQ=255)
                                                                         cAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTAACCGACACGCGCCGGTAGTAAAAATGTTCTAAATGACCAAAGCGCAGATGACTTGGTGCAACACGCATCAGCATCGCGCATGGCTCCTCCGTTTccc                                                                             >  2:2167618/1‑150 (MQ=255)
                                                                         cAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACAGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                                             <  2:2216677/150‑1 (MQ=255)
                                                                             aTGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATa                                                                         >  1:701817/1‑150 (MQ=255)
                                                                             aTGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATa                                                                         >  2:1248132/1‑150 (MQ=255)
                                                                              tgtgAACAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATac                                                                        <  2:3591756/150‑1 (MQ=255)
                                                                               gtgACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATaca                                                                       >  1:3401977/1‑150 (MQ=255)
                                                                                 gACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACt                                                                     <  2:3767200/150‑1 (MQ=255)
                                                                                 gACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACt                                                                     >  2:4883456/1‑150 (MQ=255)
                                                                                        aaTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAAt                                                              <  1:1840862/150‑1 (MQ=255)
                                                                                         aTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATc                                                             <  2:721258/150‑1 (MQ=255)
                                                                                         aTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATc                                                             <  2:3466613/150‑1 (MQ=255)
                                                                                               aaTAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                                       <  1:4557725/150‑1 (MQ=255)
                                                                                               aaTAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                                       <  1:3659962/150‑1 (MQ=255)
                                                                                               aaTAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                                       <  1:4756624/150‑1 (MQ=255)
                                                                                                 tAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                                     <  2:2464160/150‑1 (MQ=255)
                                                                                                 tAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                                     >  1:4150331/1‑150 (MQ=255)
                                                                                                 tAGCAAAGTCAGCTAACTGACGAACCTGTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                                     <  2:159326/150‑1 (MQ=255)
                                                                                                  aGCAAAGTCAGCGAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGAc                                                    >  1:1293133/1‑150 (MQ=255)
                                                                                                        gTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACt                                              >  1:153067/1‑150 (MQ=255)
                                                                                                         tCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACtt                                             >  1:3081622/1‑150 (MQ=255)
                                                                                                          cAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTTCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTa                                            >  2:1934926/1‑150 (MQ=255)
                                                                                                          cAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTa                                            >  2:1821429/1‑150 (MQ=255)
                                                                                                           aGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                           <  1:1372763/150‑1 (MQ=255)
                                                                                                           aGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                           <  1:1509479/150‑1 (MQ=255)
                                                                                                           aGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                           >  2:1048927/1‑150 (MQ=255)
                                                                                                               aaCTGACGTACCTTTTCCGACTCGCGGCGGTAGTAAATATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCTCCTGGCTCCGCCGTTTCCCGATAGACTGGCGAATCGCTGGTGACGATACTTAAcgcg                                       <  1:2251001/150‑1 (MQ=255)
                                                                                                               aaCTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                                       >  1:1497640/1‑150 (MQ=255)
                                                                                                                    aCGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgt                                  <  1:4998573/150‑1 (MQ=255)
                                                                                                                      gAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcg                                >  1:1710921/1‑150 (MQ=255)
                                                                                                                        aCCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgtc                              <  2:4760160/150‑1 (MQ=255)
                                                                                                                         ccTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGTTGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGtcgtcg                             <  1:2011830/150‑1 (MQ=255)
                                                                                                                         ccTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGtcgtcg                             >  2:1717011/1‑150 (MQ=255)
                                                                                                                         ccTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGtcgtcg                             >  1:2523801/1‑150 (MQ=255)
                                                                                                                                cGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCATAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc                      <  1:3466134/150‑1 (MQ=255)
                                                                                                                                    tCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGa                  >  1:2979961/1‑150 (MQ=255)
                                                                                                                                      gcggcgGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATa                <  1:4512311/150‑1 (MQ=255)
                                                                                                                                       cggcggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATaa               <  1:4894395/150‑1 (MQ=255)
                                                                                                                                           ggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCa           >  2:1883193/1‑150 (MQ=255)
                                                                                                                                           ggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCa           >  2:1431836/1‑150 (MQ=255)
                                                                                                                                             tagtaAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATc         >  1:3504161/1‑150 (MQ=255)
                                                                                                                                               gtaAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGACTGCCCAGATAATGCATCGc       >  2:1255589/1‑150 (MQ=255)
                                                                                                                                                taAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGcc      >  1:2105595/1‑150 (MQ=255)
                                                                                                                                                 aaaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCt     <  1:594081/150‑1 (MQ=255)
                                                                                                                                                 aaaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCt     <  1:3147078/150‑1 (MQ=255)
                                                                                                                                                   aaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCa   <  2:4556623/150‑1 (MQ=255)
                                                                                                                                                    aaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCAc  >  2:5136505/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCAC  >  CP009273/1784804‑1785101

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: