Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,933,680 G→T 27.8% S98* (TCG→TAG)  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,6800GT27.8% 761.4 / 259.6 348S98* (TCG→TAG) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base G (127/124);  new base T (61/36);  total (188/160)
Fisher's exact test for biased strand distribution p-value = 4.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

GCCCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGC  >  CP009273/3933531‑3933829
                                                                                                                                                     |                                                                                                                                                     
gCCCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCg                                                                                                                                                       <  2:4361800/150‑1 (MQ=255)
      ccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAg                                                                                                                                                 <  2:4679155/150‑1 (MQ=255)
          aCCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCAACATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAg                                                                                                                                             >  2:1694379/1‑150 (MQ=255)
            ccGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACACCCAGCTGTGGTAGCGAATTAGAAAGTg                                                                                                                                           >  2:4504611/1‑150 (MQ=255)
             cGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACACCCAGCTGTGGTAGCGACTGAGAAAGTgc                                                                                                                                          >  2:4857500/1‑150 (MQ=255)
              ggCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcg                                                                                                                                         >  1:3334154/1‑150 (MQ=255)
              ggCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcg                                                                                                                                         <  2:3190356/150‑1 (MQ=255)
              ggCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcg                                                                                                                                         <  2:3256600/150‑1 (MQ=255)
               gCATGGCGACAAAATTCAATACTTGAATAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTgcgc                                                                                                                                        <  1:4878817/150‑1 (MQ=255)
                   ggCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGc                                                                                                                                    >  2:3420497/1‑150 (MQ=255)
                         aaaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAgaga                                                                                                                              >  2:1496320/1‑150 (MQ=255)
                         aaaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAgaga                                                                                                                              >  2:4929552/1‑150 (MQ=255)
                         aaaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAgaga                                                                                                                              >  2:2946840/1‑150 (MQ=255)
                          aaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGaa                                                                                                                             >  1:4931052/1‑150 (MQ=255)
                          aaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGaa                                                                                                                             >  2:1056510/1‑150 (MQ=255)
                           aaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAc                                                                                                                            <  2:2751921/150‑1 (MQ=255)
                             ttCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAAccc                                                                                                                          >  1:2814429/1‑150 (MQ=255)
                              tCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCa                                                                                                                         >  1:167429/1‑150 (MQ=255)
                                aaTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt                                                                                                                       <  2:2642317/150‑1 (MQ=255)
                                   aCTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgtg                                                                                                                    <  2:2039231/150‑1 (MQ=255)
                                       gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACa                                                                                                                <  1:2534900/150‑1 (MQ=255)
                                       gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACa                                                                                                                <  1:88667/150‑1 (MQ=255)
                                       gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACa                                                                                                                <  2:2742503/150‑1 (MQ=255)
                                        aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAg                                                                                                               <  1:2125493/150‑1 (MQ=255)
                                           aagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACt                                                                                                            >  2:3458104/1‑150 (MQ=255)
                                             gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCa                                                                                                          >  1:4240685/1‑150 (MQ=255)
                                               tCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACg                                                                                                        <  1:1004401/150‑1 (MQ=255)
                                                cATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACgg                                                                                                       >  1:3306285/1‑150 (MQ=255)
                                                 aTTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGc                                                                                                      >  1:4469732/1‑150 (MQ=255)
                                                   tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCaa                                                                                                    >  2:212722/1‑150 (MQ=255)
                                                    aCGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAg                                                                                                   >  2:1586620/1‑150 (MQ=255)
                                                     cGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAgg                                                                                                  >  1:3788896/1‑150 (MQ=255)
                                                     cGAGGAGACGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAgg                                                                                                  >  1:1188474/1‑150 (MQ=255)
                                                      gAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAggg                                                                                                 <  2:5098675/150‑1 (MQ=255)
                                                        ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTa                                                                                               <  2:2505758/150‑1 (MQ=255)
                                                          aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaaa                                                                                             <  1:528091/150‑1 (MQ=255)
                                                           tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACTGCAAGGGTaaaa                                                                                            >  1:5121674/1‑150 (MQ=255)
                                                           tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaaaa                                                                                            <  2:3620378/150‑1 (MQ=255)
                                                               cgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGc                                                                                        >  2:1952974/1‑150 (MQ=255)
                                                               cgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGc                                                                                        <  1:1130363/150‑1 (MQ=255)
                                                               cgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGACTTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCACGGGTAAAAATGc                                                                                        >  1:4487473/1‑150 (MQ=255)
                                                                gcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCg                                                                                       >  1:4739233/1‑150 (MQ=255)
                                                                 cgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGa                                                                                      >  2:3674536/1‑150 (MQ=255)
                                                                  gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAc                                                                                     <  1:1868921/150‑1 (MQ=255)
                                                                  gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAc                                                                                     <  1:2723400/150‑1 (MQ=255)
                                                                    aGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAcgc                                                                                   >  2:1691028/1‑150 (MQ=255)
                                                                      tAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGt                                                                                 >  1:1441820/1‑150 (MQ=255)
                                                                      tAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGt                                                                                 >  2:4561234/1‑150 (MQ=255)
                                                                       aaGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTa                                                                                >  1:3258893/1‑150 (MQ=255)
                                                                        aGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTAc                                                                               <  1:2381715/150‑1 (MQ=255)
                                                                         ggCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTAcc                                                                              >  1:2455527/1‑150 (MQ=255)
                                                                         ggCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTAcc                                                                              <  2:3306284/150‑1 (MQ=255)
                                                                           ccAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGa                                                                            <  2:2139330/150‑1 (MQ=255)
                                                                            cAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGaa                                                                           <  1:2269642/150‑1 (MQ=255)
                                                                            cAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGaa                                                                           <  1:1985246/150‑1 (MQ=255)
                                                                               ccGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCg                                                                        >  2:1009856/1‑150 (MQ=255)
                                                                                cGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGa                                                                       >  1:3036108/1‑150 (MQ=255)
                                                                                 gAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGAt                                                                      >  2:1812447/1‑150 (MQ=255)
                                                                                  aGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATc                                                                     <  2:1735239/150‑1 (MQ=255)
                                                                                     aaCAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGc                                                                  >  1:4708245/1‑150 (MQ=255)
                                                                                      aCAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCt                                                                 <  2:3856267/150‑1 (MQ=255)
                                                                                         ggCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGc                                                              >  1:80564/1‑150 (MQ=255)
                                                                                         ggCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGc                                                              >  2:663940/1‑150 (MQ=255)
                                                                                            atcatcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCAt                                                           <  1:5048158/150‑1 (MQ=255)
                                                                                             tcatcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATc                                                          <  2:3816843/150‑1 (MQ=255)
                                                                                              catcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATcc                                                         <  1:1867384/150‑1 (MQ=255)
                                                                                              catcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATcc                                                         <  2:1226900/150‑1 (MQ=255)
                                                                                                 catTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCt                                                      <  2:592794/150‑1 (MQ=255)
                                                                                                  atTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCtt                                                     >  2:412832/1‑150 (MQ=255)
                                                                                                       gccgccTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCgg                                                <  2:2960261/150‑1 (MQ=255)
                                                                                                        ccgccTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGa                                               <  1:1942380/150‑1 (MQ=255)
                                                                                                          gccTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAAt                                             >  2:3401799/1‑150 (MQ=255)
                                                                                                            cTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATaa                                           <  1:3452652/150‑1 (MQ=255)
                                                                                                            cTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATaa                                           >  1:460905/1‑150 (MQ=255)
                                                                                                              aTCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGc                                         >  1:4365026/1‑150 (MQ=255)
                                                                                                                ccccTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCAt                                       >  1:3437124/1‑150 (MQ=255)
                                                                                                                ccccTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCAt                                       <  2:405597/150‑1 (MQ=255)
                                                                                                                ccccTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCAt                                       <  2:405715/150‑1 (MQ=255)
                                                                                                                  ccTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCg                                     <  1:5076544/150‑1 (MQ=255)
                                                                                                                     gAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCa                                  >  1:857595/1‑150 (MQ=255)
                                                                                                                      aaTAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAc                                 >  2:4465297/1‑150 (MQ=255)
                                                                                                                       aTAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAcc                                <  2:2437412/150‑1 (MQ=255)
                                                                                                                        tAACCCGGAATACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCg                               <  1:4079599/150‑1 (MQ=255)
                                                                                                                        tAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCg                               >  1:4526985/1‑150 (MQ=255)
                                                                                                                           cccGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCa                            >  2:4450063/1‑150 (MQ=255)
                                                                                                                            ccGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAg                           >  1:1275970/1‑150 (MQ=255)
                                                                                                                            ccGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAg                           <  1:3616683/150‑1 (MQ=255)
                                                                                                                              ggAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCg                         <  1:1622790/150‑1 (MQ=255)
                                                                                                                               gAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGt                        <  1:2965763/150‑1 (MQ=255)
                                                                                                                                     aaCCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACt                  >  1:3482773/1‑150 (MQ=255)
                                                                                                                                     aaCCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACt                  <  1:4224360/150‑1 (MQ=255)
                                                                                                                                     aaCCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACt                  <  1:2669634/150‑1 (MQ=255)
                                                                                                                                       ccAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtga                >  2:502683/1‑150 (MQ=255)
                                                                                                                                           cTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTgatgat            >  1:388466/1‑150 (MQ=255)
                                                                                                                                            tgtgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGatgatg           <  2:4125401/150‑1 (MQ=255)
                                                                                                                                              tgGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGc         >  2:1072142/1‑150 (MQ=255)
                                                                                                                                               ggTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCt        >  1:1158183/1‑150 (MQ=255)
                                                                                                                                                    cGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTg   <  2:876102/150‑1 (MQ=255)
                                                                                                                                                     gAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGc  >  2:2312914/1‑150 (MQ=255)
                                                                                                                                                     gAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGc  <  1:3119798/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGC  >  CP009273/3933531‑3933829

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: