Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,407,007 G→A 41.6% intergenic (‑183/‑216) acrS ← / → acrE acrAB operon transcriptional repressor/cytoplasmic membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,407,0070GA41.6% 169.8 / inf 247intergenic (‑183/‑216)acrS/acrEacrAB operon transcriptional repressor/cytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (71/73);  new base A (63/39);  total (134/113)
Fisher's exact test for biased strand distribution p-value = 6.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 254 total aligned reads displayed.

AGAGTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTATTTTTC  >  CP009273/3406859‑3407156
                                                                                                                                                    |                                                                                                                                                     
  agTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGtt                                                                                                                                                    >  1:3942459/1‑150 (MQ=255)
  agTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGtt                                                                                                                                                    <  1:1475451/150‑1 (MQ=255)
   gTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGttt                                                                                                                                                   <  2:74448/150‑1 (MQ=255)
    tGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGtttt                                                                                                                                                  <  2:462173/150‑1 (MQ=255)
       gaAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAAt                                                                                                                                               <  1:4161259/150‑1 (MQ=255)
        aaaaTAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAAtt                                                                                                                                              <  1:3029716/150‑1 (MQ=255)
          aaTAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCt                                                                                                                                            <  1:3616321/150‑1 (MQ=255)
            tAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGc                                                                                                                                          >  1:548536/1‑150 (MQ=255)
             aGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGcc                                                                                                                                         <  1:1360830/150‑1 (MQ=255)
                gAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCctct                                                                                                                                      <  2:3943908/150‑1 (MQ=255)
                 aaGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCtt                                                                                                                                     <  2:3495416/150‑1 (MQ=255)
                   ggTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTc                                                                                                                                   <  1:2438268/150‑1 (MQ=255)
                   ggTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTc                                                                                                                                   >  2:194666/1‑150 (MQ=255)
                      tAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAAc                                                                                                                                <  2:4811127/150‑1 (MQ=255)
                      tAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAAc                                                                                                                                <  2:2652517/150‑1 (MQ=255)
                         tCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCcg                                                                                                                             >  1:2280476/1‑150 (MQ=255)
                             tCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTc                                                                                                                         >  1:4232930/1‑150 (MQ=255)
                               aCCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCaa                                                                                                                       >  2:320686/1‑150 (MQ=255)
                               aCCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCaa                                                                                                                       >  2:1901542/1‑150 (MQ=255)
                                ccTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCaaa                                                                                                                      <  2:286993/150‑1 (MQ=255)
                                         tATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTAAACCCGCGTCAAAATAAAACAg                                                                                                             >  1:1945259/1‑150 (MQ=255)
                                         tATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAg                                                                                                             >  1:3238407/1‑150 (MQ=255)
                                              ataCGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACCGTAGaa                                                                                                        >  1:3041897/1‑150 (MQ=255)
                                              ataCGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGaa                                                                                                        >  1:1147847/1‑150 (MQ=255)
                                                   agagGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTa                                                                                                   >  1:4026607/1‑150 (MQ=255)
                                                       gCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCt                                                                                               >  1:714652/1‑150 (MQ=255)
                                                       gCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCt                                                                                               <  1:1029326/150‑1 (MQ=255)
                                                        cTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCtt                                                                                              >  2:3118018/1‑150 (MQ=255)
                                                         tAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCttt                                                                                             <  2:2228689/150‑1 (MQ=255)
                                                            tatTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCtttttt                                                                                          <  2:2753684/150‑1 (MQ=255)
                                                            tatTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCtttttt                                                                                          >  1:3827993/1‑150 (MQ=255)
                                                               tGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATGTTTTTTtgt                                                                                       >  2:1256800/1‑150 (MQ=255)
                                                                gAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTtgtg                                                                                      <  1:2722362/150‑1 (MQ=255)
                                                                  gTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGtt                                                                                    <  2:1394927/150‑1 (MQ=255)
                                                                   ttGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGttt                                                                                   <  2:193340/150‑1 (MQ=255)
                                                                   ttGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGttt                                                                                   >  1:3525815/1‑150 (MQ=255)
                                                                    tGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTa                                                                                  <  1:2808837/150‑1 (MQ=255)
                                                                     gCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAt                                                                                 <  2:3939459/150‑1 (MQ=255)
                                                                      cTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtg                                                                                >  2:4981098/1‑150 (MQ=255)
                                                                      cTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtg                                                                                <  2:4481321/150‑1 (MQ=255)
                                                                      cTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtg                                                                                >  2:448583/1‑150 (MQ=255)
                                                                       tataAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtgt                                                                               >  1:3661723/1‑150 (MQ=255)
                                                                       tataAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtgt                                                                               >  1:1738253/1‑150 (MQ=255)
                                                                        ataAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtgtg                                                                              >  1:4429803/1‑150 (MQ=255)
                                                                        ataAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtgtg                                                                              >  1:4430098/1‑150 (MQ=255)
                                                                        ataAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtgtg                                                                              >  2:4547770/1‑150 (MQ=255)
                                                                         taAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGc                                                                             <  1:2052540/150‑1 (MQ=255)
                                                                         taAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGc                                                                             <  1:5001624/150‑1 (MQ=255)
                                                                          aaaTCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGcc                                                                            >  2:522507/1‑150 (MQ=255)
                                                                           aaTCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCt                                                                           >  2:3904574/1‑150 (MQ=255)
                                                                           aaTCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCt                                                                           >  2:2087222/1‑150 (MQ=255)
                                                                            aTCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCtt                                                                          <  1:2462744/150‑1 (MQ=255)
                                                                            aTCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCtt                                                                          >  1:2079275/1‑150 (MQ=255)
                                                                             tCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTg                                                                         >  2:2034947/1‑150 (MQ=255)
                                                                               gTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTgag                                                                       >  1:223860/1‑150 (MQ=255)
                                                                               gTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTgag                                                                       >  1:4187219/1‑150 (MQ=255)
                                                                                 taaataaaTAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGAt                                                                     <  1:1683755/150‑1 (MQ=255)
                                                                                            tatatatTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCa                                                          >  1:3457029/1‑150 (MQ=255)
                                                                                               atatTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATaa                                                       <  1:1452552/150‑1 (MQ=255)
                                                                                                 attattTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACt                                                     <  1:1990524/150‑1 (MQ=255)
                                                                                                   tattTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCGACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTAt                                                   >  2:2174446/1‑150 (MQ=255)
                                                                                                    attTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTAtt                                                  <  1:1536014/150‑1 (MQ=255)
                                                                                                    attTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTAtt                                                  <  1:951861/150‑1 (MQ=255)
                                                                                                       tACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCt                                               >  1:4053922/1‑150 (MQ=255)
                                                                                                        aCCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCtt                                              <  2:1487047/150‑1 (MQ=255)
                                                                                                            aaGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACAt                                          <  1:159017/150‑1 (MQ=255)
                                                                                                             aGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATc                                         <  2:713296/150‑1 (MQ=255)
                                                                                                               aTACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcga                                       >  2:105848/1‑150 (MQ=255)
                                                                                                               aTACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcga                                       >  2:63298/1‑150 (MQ=255)
                                                                                                                 aCATTCACTACATCAATATATATTTCTATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcgacg                                     <  1:3407076/150‑1 (MQ=255)
                                                                                                                 aCATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcgacg                                     >  1:1006066/1‑150 (MQ=255)
                                                                                                                  cATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcgacga                                    >  1:3930957/1‑150 (MQ=255)
                                                                                                                    ttCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAAt                                  >  2:3963120/1‑150 (MQ=255)
                                                                                                                      cACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTGATAACTATTCCTTACATCGACGAATGa                                <  2:5001561/150‑1 (MQ=255)
                                                                                                                       aCTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGAt                               >  2:2059732/1‑150 (MQ=255)
                                                                                                                        cTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATa                              >  1:885692/1‑150 (MQ=255)
                                                                                                                        cTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATa                              <  2:3518932/150‑1 (MQ=255)
                                                                                                                           cATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAATCAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAAtt                           <  1:5144960/150‑1 (MQ=255)
                                                                                                                            aTCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAAttt                          >  2:913835/1‑150 (MQ=255)
                                                                                                                            aTCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAAttt                          >  2:3018970/1‑150 (MQ=255)
                                                                                                                             tCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTg                         >  2:4381999/1‑150 (MQ=255)
                                                                                                                             gcAATATATATTTCAATTTACGAGGATTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTg                         <  2:4793311/149‑1 (MQ=255)
                                                                                                                              cAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGt                        <  1:68588/150‑1 (MQ=255)
                                                                                                                                 tatatatTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAgg                     >  1:563413/1‑150 (MQ=255)
                                                                                                                                   tatatTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGAt                   >  1:2997950/1‑150 (MQ=255)
                                                                                                                                    atatTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATTATTTGTAGGATa                  <  2:3336626/150‑1 (MQ=255)
                                                                                                                                    atatTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATa                  >  1:181347/1‑150 (MQ=255)
                                                                                                                                           aaTTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACt           >  2:2135973/1‑150 (MQ=255)
                                                                                                                                           aaTTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACt           <  1:5091518/150‑1 (MQ=255)
                                                                                                                                           aaTTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACt           <  2:2301146/150‑1 (MQ=255)
                                                                                                                                            ntttACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTg          >  2:2935075/2‑150 (MQ=255)
                                                                                                                                            aTTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTg          >  2:4533861/1‑150 (MQ=255)
                                                                                                                                             tttACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGt         <  1:3118018/150‑1 (MQ=255)
                                                                                                                                               tACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAt       >  2:3057565/1‑150 (MQ=255)
                                                                                                                                                 cGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAttt     >  1:2686155/1‑150 (MQ=255)
                                                                                                                                                 cGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAttt     >  2:1747192/1‑150 (MQ=255)
                                                                                                                                                  gagaTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAtttt    >  2:1213096/1‑150 (MQ=255)
                                                                                                                                                  gagaTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAtttt    >  2:1033400/1‑150 (MQ=255)
                                                                                                                                                    ggTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAttttta  >  2:4307655/1‑149 (MQ=255)
                                                                                                                                                    ggTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTATTtttc  <  1:2795926/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGAGTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTATTTTTC  >  CP009273/3406859‑3407156

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: