Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,653,713 T→A 18.7% V463E (GTG→GAG)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,653,7130TA18.7% 515.7 / 115.9 198V463E (GTG→GAG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (76/85);  new base A (14/23);  total (90/108)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 201 total aligned reads displayed.

ACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATAT  >  CP009273/1653565‑1653861
                                                                                                                                                    |                                                                                                                                                    
aCCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGAtgtg                                                                                                                                                     <  1:4200340/150‑1 (MQ=255)
    ataCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCa                                                                                                                                                 <  2:1759039/150‑1 (MQ=255)
     taccattaccatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCAt                                                                                                                                                >  1:302760/1‑150 (MQ=255)
       ccattaccatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTa                                                                                                                                              >  1:3604130/1‑150 (MQ=255)
       ccattaccatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTa                                                                                                                                              <  2:1897182/150‑1 (MQ=255)
       ccattaccatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGACAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTa                                                                                                                                              >  1:548707/1‑150 (MQ=255)
         attaccatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAg                                                                                                                                            <  2:1746055/150‑1 (MQ=255)
           taccatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGtt                                                                                                                                          >  2:3334037/1‑150 (MQ=255)
           taccatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGAGAAACTGGATGTGCCCATTAAGtt                                                                                                                                          >  2:3897321/1‑150 (MQ=255)
            accatAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTc                                                                                                                                         <  1:3149791/150‑1 (MQ=255)
               atAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATc                                                                                                                                      >  1:4426177/1‑150 (MQ=255)
                tAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCt                                                                                                                                     <  1:868117/150‑1 (MQ=255)
                  gAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTgg                                                                                                                                   >  1:2213087/1‑150 (MQ=255)
                   aaCGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGa                                                                                                                                  >  2:1786019/1‑150 (MQ=255)
                    aCGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGaa                                                                                                                                 <  1:3359816/150‑1 (MQ=255)
                       ccTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTa                                                                                                                              <  1:4379548/150‑1 (MQ=255)
                        cTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTAc                                                                                                                             >  1:494635/1‑150 (MQ=255)
                         tGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTAcg                                                                                                                            <  2:2965685/150‑1 (MQ=255)
                                gCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGATGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGaaa                                                                                                                     <  2:3319453/150‑1 (MQ=255)
                                 cTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAAt                                                                                                                    <  2:1072690/150‑1 (MQ=255)
                                  tGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATa                                                                                                                   <  1:994319/150‑1 (MQ=255)
                                   ggATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATAc                                                                                                                  <  2:4264755/150‑1 (MQ=255)
                                    gATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATAcc                                                                                                                 >  1:427726/1‑150 (MQ=255)
                                      tAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCt                                                                                                               >  1:4317982/1‑150 (MQ=255)
                                      tAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCt                                                                                                               <  1:384966/150‑1 (MQ=255)
                                       aaTCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTc                                                                                                              <  1:3572594/150‑1 (MQ=255)
                                        aTCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCg                                                                                                             <  2:4297561/150‑1 (MQ=255)
                                         tCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGt                                                                                                            <  2:3638066/150‑1 (MQ=255)
                                          ccGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGtt                                                                                                           >  1:4318646/1‑150 (MQ=255)
                                          ccGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGtt                                                                                                           <  2:3540838/150‑1 (MQ=255)
                                            ggTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTaa                                                                                                         >  2:4554476/1‑150 (MQ=255)
                                            ggTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTaa                                                                                                         >  2:1374709/1‑150 (MQ=255)
                                            ggTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTaa                                                                                                         >  1:2800068/1‑150 (MQ=255)
                                             gTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAAt                                                                                                        >  2:598472/1‑150 (MQ=255)
                                             gTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAAt                                                                                                        <  2:3431398/150‑1 (MQ=255)
                                              tCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATc                                                                                                       <  1:3821736/150‑1 (MQ=255)
                                               cAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCa                                                                                                      >  2:2018905/1‑150 (MQ=255)
                                               cAAAACGTCAATATCCAGCTTCAGCTGGACAGATGCTAACGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGATCATCTGGAACTACGCGGGAAATACCCTCGATAATCa                                                                                                      >  2:305801/1‑150 (MQ=255)
                                                aaaaCGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAg                                                                                                     >  2:3249059/1‑150 (MQ=255)
                                                 aaaCGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGc                                                                                                    <  2:3713321/150‑1 (MQ=255)
                                                  aaCGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCa                                                                                                   <  1:4257556/150‑1 (MQ=255)
                                                   aCGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCAt                                                                                                  <  1:753348/150‑1 (MQ=255)
                                                    cGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCAtt                                                                                                 <  2:2198813/150‑1 (MQ=255)
                                                     gTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTc                                                                                                <  2:3270188/150‑1 (MQ=255)
                                                     gTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTc                                                                                                <  2:896928/150‑1 (MQ=255)
                                                     gTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTc                                                                                                >  1:3246441/1‑150 (MQ=255)
                                                      tCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCt                                                                                               <  2:2373106/150‑1 (MQ=255)
                                                      tCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCt                                                                                               <  1:1599420/150‑1 (MQ=255)
                                                       cAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTg                                                                                              >  2:2767098/1‑150 (MQ=255)
                                                           ctcCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACAt                                                                                          <  1:3118569/150‑1 (MQ=255)
                                                            tcCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATc                                                                                         >  2:4008264/1‑150 (MQ=255)
                                                            tcCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATc                                                                                         >  2:1621856/1‑150 (MQ=255)
                                                               tGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATcaac                                                                                      <  2:4211976/150‑1 (MQ=255)
                                                                   tCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACaaaa                                                                                  >  2:437061/1‑150 (MQ=255)
                                                                    cAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAAc                                                                                 >  2:2163810/1‑150 (MQ=255)
                                                                      gCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGc                                                                               >  1:3216041/1‑150 (MQ=255)
                                                                       cTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCa                                                                              <  1:1330661/150‑1 (MQ=255)
                                                                       cTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCa                                                                              <  1:1330517/150‑1 (MQ=255)
                                                                        tGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCAt                                                                             >  2:1565229/1‑150 (MQ=255)
                                                                         ggACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATg                                                                            <  1:1197572/150‑1 (MQ=255)
                                                                           aCAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGaa                                                                          <  2:4038015/150‑1 (MQ=255)
                                                                           aCAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGaa                                                                          <  2:1847419/150‑1 (MQ=255)
                                                                             aGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAAt                                                                        <  2:735606/150‑1 (MQ=255)
                                                                                  cTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGc                                                                   <  1:4298742/150‑1 (MQ=255)
                                                                                  cTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGc                                                                   >  1:4393944/1‑150 (MQ=255)
                                                                                   tATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCa                                                                  >  2:3299511/1‑150 (MQ=255)
                                                                                    aTCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAg                                                                 <  1:1647285/150‑1 (MQ=255)
                                                                                          cATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGaa                                                           <  2:463885/150‑1 (MQ=255)
                                                                                           aTGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAAt                                                          <  1:2152416/150‑1 (MQ=255)
                                                                                                ccGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaa                                                     >  1:1400526/1‑150 (MQ=255)
                                                                                                ccGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaa                                                     >  1:1400334/1‑150 (MQ=255)
                                                                                                 cGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaat                                                    <  1:2430611/150‑1 (MQ=255)
                                                                                                   cAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgc                                                  <  2:302760/150‑1 (MQ=255)
                                                                                                     aaTGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgcga                                                >  2:4739908/1‑150 (MQ=255)
                                                                                                      aTGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgcgaa                                               >  1:258791/1‑150 (MQ=255)
                                                                                                      aTGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgcgaa                                               >  1:757665/1‑150 (MQ=255)
                                                                                                       tGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgcgaaa                                              <  2:494453/150‑1 (MQ=255)
                                                                                                        gACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGGAATTCTGCAGGACGAATcgaaatgcgaaat                                             >  2:2322589/1‑150 (MQ=255)
                                                                                                           ggCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGAt                                          >  2:2788981/1‑150 (MQ=255)
                                                                                                               aTCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATtgtg                                      <  1:1821723/150‑1 (MQ=255)
                                                                                                                tCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTgg                                     <  1:2381844/150‑1 (MQ=255)
                                                                                                                tCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTgg                                     >  1:3190352/1‑150 (MQ=255)
                                                                                                                  cgcgACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGAGTGTGGTc                                   >  1:4621545/1‑150 (MQ=255)
                                                                                                                  cgcgACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCAGGAAATTCTGCAGGACGAATCGAAATGCGAAATGAGTGTGGTc                                   >  2:723371/1‑150 (MQ=255)
                                                                                                                     gACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATc                                >  1:532821/1‑150 (MQ=255)
                                                                                                                           gTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAAc                          >  2:2631195/1‑150 (MQ=255)
                                                                                                                           gTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAAc                          <  2:1589631/150‑1 (MQ=255)
                                                                                                                           gTCCGCGGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAAc                          <  2:1589580/150‑1 (MQ=255)
                                                                                                                            tCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACt                         <  1:3762785/150‑1 (MQ=255)
                                                                                                                              cgcgGCAAAGATAAACTGGATGAGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACttt                       <  2:794506/150‑1 (MQ=255)
                                                                                                                                cgGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTAt                     >  2:1257173/1‑150 (MQ=255)
                                                                                                                                  gCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGa                   <  2:150218/150‑1 (MQ=255)
                                                                                                                                   cAAAGATACACTGGATGTGCCCATTCACTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAACCGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGAc                  <  1:1172724/150‑1 (MQ=255)
                                                                                                                                     aaGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACGACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCt                <  2:2923939/150‑1 (MQ=255)
                                                                                                                                        aTAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCATGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCAt             >  1:3547075/1‑150 (MQ=255)
                                                                                                                                         tAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATc            >  2:4518938/1‑150 (MQ=255)
                                                                                                                                              tGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCaa       >  2:4486209/1‑150 (MQ=255)
                                                                                                                                                 aTGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAat    <  1:1492082/150‑1 (MQ=255)
                                                                                                                                                  tgtgCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAata   >  1:4595428/1‑150 (MQ=255)
                                                                                                                                                   gtgCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAatat  <  1:1744120/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATAT  >  CP009273/1653565‑1653861

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: