Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,933,926 G→T 6.4% A16E (GCG→GAG)  hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,933,9260GT 6.4% 971.2 / 28.3 250A16E (GCG→GAG) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base G (107/127);  new base T (9/7);  total (116/134)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 262 total aligned reads displayed.

TACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAAT  >  CP009273/3933778‑3934075
                                                                                                                                                    |                                                                                                                                                     
tACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTCAGAATAGTGGCATCAAGTGCCTGCAGGAAGAAGGCCATcgccgc                                                                                                                                                      >  1:118239/1‑150 (MQ=255)
    ggAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATc                                                                                                                                                  <  1:1916209/150‑1 (MQ=255)
    ggAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATc                                                                                                                                                  >  1:3817112/1‑150 (MQ=255)
       aTAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCAc                                                                                                                                               <  2:2201393/150‑1 (MQ=255)
       aTAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCTAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCAc                                                                                                                                               <  1:587439/150‑1 (MQ=255)
        tAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCTCGATCCACg                                                                                                                                              >  1:1867574/1‑150 (MQ=255)
         aaGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACgg                                                                                                                                             <  1:2160485/150‑1 (MQ=255)
         aaGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACgg                                                                                                                                             <  1:2160466/150‑1 (MQ=255)
          aGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGc                                                                                                                                            >  2:3275424/1‑150 (MQ=255)
               cGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAAcc                                                                                                                                       >  1:2533139/1‑150 (MQ=255)
                gCCACCGTCAGCGTATAACTGATGATGGCTTATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCAATTGCGGGTAAGGCGGTATTAAGAATAGTGTCATCAAGTGCCTGCATGAAGAAGGCAATCGACGCGATCCACGGCAAAccc                                                                                                                                      >  2:676195/1‑150 (MQ=255)
                 ccACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCg                                                                                                                                     <  2:1351031/150‑1 (MQ=255)
                 ccACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCg                                                                                                                                     <  1:3369110/150‑1 (MQ=255)
                  cACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGc                                                                                                                                    >  1:3874304/1‑150 (MQ=255)
                    ccGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCa                                                                                                                                  <  2:256871/150‑1 (MQ=255)
                      gTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATa                                                                                                                                <  1:4146182/150‑1 (MQ=255)
                        cAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACt                                                                                                                              <  2:1913066/150‑1 (MQ=255)
                        cAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACt                                                                                                                              <  2:2214858/150‑1 (MQ=255)
                        cAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACt                                                                                                                              >  1:3699856/1‑150 (MQ=255)
                        cAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACt                                                                                                                              <  1:3145484/150‑1 (MQ=255)
                         aGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTg                                                                                                                             <  2:2388565/150‑1 (MQ=255)
                                aaCTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGcttct                                                                                                                      <  1:2891233/150‑1 (MQ=255)
                                aaCTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGcttct                                                                                                                      <  1:469653/150‑1 (MQ=255)
                                aaCTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGcttct                                                                                                                      <  2:3237476/150‑1 (MQ=255)
                                 aCTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGcttctt                                                                                                                     <  2:1556716/150‑1 (MQ=255)
                                 aCTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGcttctt                                                                                                                     <  2:837539/150‑1 (MQ=255)
                                 aCTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGcttctt                                                                                                                     <  2:837651/150‑1 (MQ=255)
                                    gatgatGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCttttt                                                                                                                  >  1:2798484/1‑150 (MQ=255)
                                    gatgatGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCttttt                                                                                                                  <  2:708996/150‑1 (MQ=255)
                                     atgatgGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTa                                                                                                                 >  2:183134/1‑150 (MQ=255)
                                       gatgGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATc                                                                                                               <  2:3705044/150‑1 (MQ=255)
                                        atgGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCTCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCg                                                                                                              <  1:2649348/150‑1 (MQ=255)
                                        atgGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCg                                                                                                              <  1:4001254/150‑1 (MQ=255)
                                         tgGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGc                                                                                                             >  2:3079441/1‑150 (MQ=255)
                                             tGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCAt                                                                                                         <  1:1068467/150‑1 (MQ=255)
                                              gATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCAtt                                                                                                        <  2:4050565/150‑1 (MQ=255)
                                              gATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCAtt                                                                                                        >  2:3597132/1‑150 (MQ=255)
                                              gATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCAtt                                                                                                        <  1:1382854/150‑1 (MQ=255)
                                                ttGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCTCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCa                                                                                                      <  2:2906658/150‑1 (MQ=255)
                                                 tGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCaa                                                                                                     >  2:946795/1‑150 (MQ=255)
                                                  gCATCGCGAGTGGTGTACGATTATTGCTATGAGCTATTGCTGGTAAGGCGGTATTCTGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCAGCGATCCACGGCATACCCGCCATTCTGCGGTTCTTTTTATCGCTCATTCAAt                                                                                                    <  1:272016/150‑1 (MQ=255)
                                                   cATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCTCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATg                                                                                                   >  1:4166017/1‑150 (MQ=255)
                                                      cgcgAGAGGAGAACGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTcc                                                                                                <  1:2778511/150‑1 (MQ=255)
                                                         gagaGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGt                                                                                             >  1:3340635/1‑150 (MQ=255)
                                                          agatGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAATTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACACGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGtt                                                                                            <  1:1108012/150‑1 (MQ=255)
                                                          agagGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGACTAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGtt                                                                                            <  1:3485058/150‑1 (MQ=255)
                                                            aggagAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTAt                                                                                          >  1:3030137/1‑150 (MQ=255)
                                                             ggagAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATc                                                                                         >  2:3172992/1‑150 (MQ=255)
                                                                gaACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCggg                                                                                      >  1:2309750/1‑150 (MQ=255)
                                                                   cGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTa                                                                                   <  1:3717243/150‑1 (MQ=255)
                                                                   cGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTa                                                                                   <  1:3717296/150‑1 (MQ=255)
                                                                   cGATTAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTa                                                                                   <  2:3961046/150‑1 (MQ=255)
                                                                    gATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTAt                                                                                  >  2:2937173/1‑150 (MQ=255)
                                                                     aTTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATc                                                                                 <  2:1683920/150‑1 (MQ=255)
                                                                      ttAAGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCa                                                                                <  2:1193936/150‑1 (MQ=255)
                                                                        aaGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACt                                                                              >  1:199682/1‑150 (MQ=255)
                                                                         aGGCTATGAGCTATTGCGGGTAAGGCGGTATTAAGTATAGTGGCATCAAGTGCCTTCATGTAGAAGGCCATCGCCGCGATCCACTGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACtt                                                                             <  1:1151496/150‑1 (MQ=255)
                                                                            cTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATc                                                                          >  2:1714940/1‑150 (MQ=255)
                                                                                gAGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGt                                                                      >  1:869011/1‑150 (MQ=255)
                                                                                 aGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTg                                                                     >  2:1469457/1‑150 (MQ=255)
                                                                                 aGCTATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTg                                                                     >  1:2364782/1‑150 (MQ=255)
                                                                                 aGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTg                                                                     <  1:3299422/150‑1 (MQ=255)
                                                                                    taTTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGc                                                                  <  2:1004004/149‑1 (MQ=255)
                                                                                    gATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGc                                                                  <  1:1493506/150‑1 (MQ=255)
                                                                                     aTTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCg                                                                 >  1:142977/1‑150 (MQ=255)
                                                                                       tGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTa                                                               <  2:3817112/150‑1 (MQ=255)
                                                                                       tGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTa                                                               >  1:834104/1‑150 (MQ=255)
                                                                                          gggTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTagca                                                            <  2:2319858/150‑1 (MQ=255)
                                                                                          gggTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTagca                                                            <  1:701948/150‑1 (MQ=255)
                                                                                           ggTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTagcag                                                           <  1:3286224/150‑1 (MQ=255)
                                                                                           ggTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTagcag                                                           <  1:2405903/150‑1 (MQ=255)
                                                                                              aaGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAgcgc                                                        >  2:2806948/1‑150 (MQ=255)
                                                                                              aaGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGAAGCAgcgc                                                        >  1:3532759/1‑150 (MQ=255)
                                                                                                ggCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCt                                                      >  2:3887580/1‑150 (MQ=255)
                                                                                                  cGGTATTAAGAAGAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGa                                                    <  2:2280414/150‑1 (MQ=255)
                                                                                                    gTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACa                                                  >  1:3607018/1‑150 (MQ=255)
                                                                                                               aGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATg                                       >  1:745687/1‑150 (MQ=255)
                                                                                                               aGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATg                                       >  2:3924264/1‑150 (MQ=255)
                                                                                                                  ggCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCcg                                    <  2:1077753/150‑1 (MQ=255)
                                                                                                                     aTCAAGTGCCTGCATGAAGAAGGCCATCGCCTCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGt                                 <  2:3459281/150‑1 (MQ=255)
                                                                                                                     aTCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGt                                 >  1:3214056/1‑150 (MQ=255)
                                                                                                                        aaGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGc                              >  1:1677910/1‑150 (MQ=255)
                                                                                                                           tGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCAt                           >  2:3671689/1‑150 (MQ=255)
                                                                                                                            gCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATc                          <  2:3996516/150‑1 (MQ=255)
                                                                                                                              cTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGc                        >  1:3682813/1‑150 (MQ=255)
                                                                                                                                gCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCtt                      >  2:2422505/1‑150 (MQ=255)
                                                                                                                                gCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCtt                      <  1:3099167/150‑1 (MQ=255)
                                                                                                                                 cATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCttt                     >  1:1775964/1‑150 (MQ=255)
                                                                                                                                      agaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATa                >  1:1000488/1‑150 (MQ=255)
                                                                                                                                       gaagGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATaa               <  2:748163/150‑1 (MQ=255)
                                                                                                                                         agGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATc             <  2:1392928/150‑1 (MQ=255)
                                                                                                                                             cATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCAt         >  2:3469806/1‑150 (MQ=255)
                                                                                                                                              aTCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATc        <  1:557355/150‑1 (MQ=255)
                                                                                                                                               tCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCg       <  1:3016358/150‑1 (MQ=255)
                                                                                                                                                 gccgcGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGAc     >  2:3759440/1‑150 (MQ=255)
                                                                                                                                                 gccgcGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGAc     >  2:1661968/1‑150 (MQ=255)
                                                                                                                                                  ccgcGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACa    >  1:2478091/1‑150 (MQ=255)
                                                                                                                                                   cgcgATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACaa   >  2:300926/1‑150 (MQ=255)
                                                                                                                                                   cgcgATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACaa   >  2:4086741/1‑150 (MQ=255)
                                                                                                                                                    gcgATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAAt  <  2:3708328/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGCATCGCGAGAGGAGAACGATTAAGGCTATGAGCGATTGCGGGTAAGGCGGTATTAAGAATAGTGGCATCAAGTGCCTGCATGAAGAAGGCCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAAT  >  CP009273/3933778‑3934075

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: