Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,933,988 C→A 6.8% intergenic (‑16/+7) hsrA ← / ← yieP putative multidrug or homocysteine efflux system/putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,933,9880CA 6.8% 1231.6 / 46.7 325intergenic (‑16/+7)hsrA/yiePputative multidrug or homocysteine efflux system/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (159/144);  new base A (15/7);  total (174/151)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 329 total aligned reads displayed.

GGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGAT  >  CP009273/3933839‑3934137
                                                                                                                                                     |                                                                                                                                                     
  agaACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCgg                                                                                                                                                     >  2:121527/1‑150 (MQ=255)
   gaACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCggg                                                                                                                                                    <  1:3954532/150‑1 (MQ=255)
   gaACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCggg                                                                                                                                                    >  2:3895928/1‑150 (MQ=255)
      cGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTa                                                                                                                                                 >  1:1105478/1‑150 (MQ=255)
        aTTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATc                                                                                                                                               >  2:2007075/1‑150 (MQ=255)
         ttAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCa                                                                                                                                              <  1:4168054/150‑1 (MQ=255)
             ggCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTa                                                                                                                                          <  1:838727/150‑1 (MQ=255)
                tATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCa                                                                                                                                       <  1:2369781/150‑1 (MQ=255)
                tATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCa                                                                                                                                       <  2:2629194/150‑1 (MQ=255)
                 aTGATCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAg                                                                                                                                      <  1:3660660/150‑1 (MQ=255)
                  tGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAgg                                                                                                                                     >  2:1571082/1‑150 (MQ=255)
                   gAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGt                                                                                                                                    <  2:3203091/150‑1 (MQ=255)
                     gCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGa                                                                                                                                  <  1:2833623/150‑1 (MQ=255)
                     gCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGa                                                                                                                                  <  1:2092857/150‑1 (MQ=255)
                     gCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGa                                                                                                                                  <  1:4509229/150‑1 (MQ=255)
                       gATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGc                                                                                                                                >  1:4720081/1‑150 (MQ=255)
                        aTTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCg                                                                                                                               <  1:2443938/150‑1 (MQ=255)
                               gTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTagcagc                                                                                                                        <  2:3635176/150‑1 (MQ=255)
                                 aaGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAgcgc                                                                                                                      >  2:107181/1‑150 (MQ=255)
                                 aaGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAgcgc                                                                                                                      >  2:2103733/1‑150 (MQ=255)
                                   ggCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCt                                                                                                                    <  1:884365/150‑1 (MQ=255)
                                     cGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGa                                                                                                                  >  2:438678/1‑150 (MQ=255)
                                     cGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATAGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGa                                                                                                                  <  2:1476934/150‑1 (MQ=255)
                                      ggTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGAc                                                                                                                 >  2:1764055/1‑150 (MQ=255)
                                      ggTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGAc                                                                                                                 >  1:552763/1‑150 (MQ=255)
                                           tAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCt                                                                                                            >  1:4175628/1‑150 (MQ=255)
                                           tAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCt                                                                                                            >  2:4610072/1‑150 (MQ=255)
                                           tAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATAGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCt                                                                                                            <  2:3641315/150‑1 (MQ=255)
                                             aGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCttt                                                                                                          <  2:3084157/150‑1 (MQ=255)
                                             aGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCttt                                                                                                          <  1:1021302/150‑1 (MQ=255)
                                                 tAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAAt                                                                                                      <  1:70874/150‑1 (MQ=255)
                                                       cATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCcgcg                                                                                                >  1:1926000/1‑150 (MQ=255)
                                                            aGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGcc                                                                                           <  2:1981967/150‑1 (MQ=255)
                                                             gTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCa                                                                                          >  2:4179156/1‑150 (MQ=255)
                                                             gTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCa                                                                                          >  1:2379644/1‑150 (MQ=255)
                                                                ccTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCg                                                                                       <  2:1005672/150‑1 (MQ=255)
                                                                 cTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACCAGCTTTAAATGCCGCGTCGCCATCGc                                                                                      >  2:257769/1‑150 (MQ=255)
                                                                    cATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCttt                                                                                   <  2:2878625/150‑1 (MQ=255)
                                                                    cATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCttt                                                                                   <  2:1920693/150‑1 (MQ=255)
                                                                      tGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTgg                                                                                 >  2:4236308/1‑150 (MQ=255)
                                                                        aagaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGAt                                                                               <  2:2664423/150‑1 (MQ=255)
                                                                        aagaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGAt                                                                               >  1:3295252/1‑150 (MQ=255)
                                                                        aagaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGAt                                                                               <  1:1651122/150‑1 (MQ=255)
                                                                          gaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATaa                                                                             >  1:2705919/1‑150 (MQ=255)
                                                                          gaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATaa                                                                             >  1:1665417/1‑150 (MQ=255)
                                                                          gaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATaa                                                                             <  2:4215086/150‑1 (MQ=255)
                                                                            agGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATc                                                                           >  2:692884/1‑150 (MQ=255)
                                                                            agGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATc                                                                           <  2:760975/150‑1 (MQ=255)
                                                                             ggCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCg                                                                          >  1:4566185/1‑150 (MQ=255)
                                                                                cATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCAt                                                                       <  2:2579721/150‑1 (MQ=255)
                                                                                cATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCAt                                                                       >  1:357547/1‑150 (MQ=255)
                                                                                 aTCGCCGCGATCCACGGCACACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATc                                                                      <  2:4684994/150‑1 (MQ=255)
                                                                                 aTCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATc                                                                      <  1:438678/150‑1 (MQ=255)
                                                                                       gcgATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATAGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAAt                                                                <  2:3470900/150‑1 (MQ=255)
                                                                                        cgATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATc                                                               <  1:3385238/150‑1 (MQ=255)
                                                                                         gATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCg                                                              >  1:3408864/1‑150 (MQ=255)
                                                                                          aTCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGc                                                             >  1:1972965/1‑150 (MQ=255)
                                                                                              aCGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGa                                                         <  2:488307/150‑1 (MQ=255)
                                                                                               cGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGAt                                                        >  1:2893334/1‑150 (MQ=255)
                                                                                               cGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGAt                                                        >  1:2277901/1‑150 (MQ=255)
                                                                                                  cAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCa                                                     <  1:4041009/150‑1 (MQ=255)
                                                                                                   aaaCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAg                                                    >  1:507469/1‑150 (MQ=255)
                                                                                                    aaCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGa                                                   >  1:267723/1‑150 (MQ=255)
                                                                                                      cccGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATc                                                 >  2:1897740/1‑150 (MQ=255)
                                                                                                       ccGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATcc                                                >  2:3250187/1‑150 (MQ=255)
                                                                                                        cGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCa                                               <  1:1326370/150‑1 (MQ=255)
                                                                                                        cGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCa                                               <  1:249488/150‑1 (MQ=255)
                                                                                                         gCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAg                                              <  1:4306797/150‑1 (MQ=255)
                                                                                                           cATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCt                                            <  1:3869704/150‑1 (MQ=255)
                                                                                                            aTACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCtt                                           >  2:2280631/1‑150 (MQ=255)
                                                                                                            aTACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCtt                                           <  1:2613441/150‑1 (MQ=255)
                                                                                                             tACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCttt                                          >  2:2429463/1‑150 (MQ=255)
                                                                                                              aCTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATAGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTa                                         >  2:2905854/1‑150 (MQ=255)
                                                                                                                  cgcTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCAc                                     <  1:4053871/150‑1 (MQ=255)
                                                                                                                  cgcTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCAc                                     >  1:3259760/1‑150 (MQ=255)
                                                                                                                   gcTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACt                                    >  1:3148647/1‑150 (MQ=255)
                                                                                                                     ttcttTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACtgt                                  <  1:609859/150‑1 (MQ=255)
                                                                                                                     ttcttTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACtgt                                  <  1:1506917/150‑1 (MQ=255)
                                                                                                                        tttttATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTc                               >  1:1621654/1‑150 (MQ=255)
                                                                                                                        tttttATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTc                               >  2:4236573/1‑150 (MQ=255)
                                                                                                                        tttttATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTGTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTc                               >  2:4649514/1‑150 (MQ=255)
                                                                                                                          tttATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATTAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGc                             >  2:4265452/1‑150 (MQ=255)
                                                                                                                              tCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGt                         >  2:3149467/1‑150 (MQ=255)
                                                                                                                              tCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGt                         >  2:3757882/1‑150 (MQ=255)
                                                                                                                                    aTTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGa                   <  2:2775593/150‑1 (MQ=255)
                                                                                                                                    aTTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGa                   <  1:2208221/150‑1 (MQ=255)
                                                                                                                                     ttCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGAc                  <  2:2343540/150‑1 (MQ=255)
                                                                                                                                      tCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACg                 <  1:3331560/150‑1 (MQ=255)
                                                                                                                                        aaTGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTg               <  2:62420/150‑1 (MQ=255)
                                                                                                                                        aaTGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTg               >  2:275209/1‑150 (MQ=255)
                                                                                                                                          tGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGaa             <  1:1055403/150‑1 (MQ=255)
                                                                                                                                          tGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGaa             <  2:1413187/150‑1 (MQ=255)
                                                                                                                                              cTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGaagtaa         <  2:2496315/150‑1 (MQ=255)
                                                                                                                                              cTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGaagtaa         >  2:700715/1‑150 (MQ=255)
                                                                                                                                                  tATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAgtgt     >  1:1372547/1‑150 (MQ=255)
                                                                                                                                                  tATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAgtgt     >  2:4017179/1‑150 (MQ=255)
                                                                                                                                                  tATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAgtgt     >  1:2546221/1‑150 (MQ=255)
                                                                                                                                                    tCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGa   <  1:308659/150‑1 (MQ=255)
                                                                                                                                                     cGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGAt  <  2:2947224/150‑1 (MQ=255)
                                                                                                                                                     cGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGAt  >  1:1207297/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGAT  >  CP009273/3933839‑3934137

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: