Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,407,008 G→A 27.9% intergenic (‑184/‑215) 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,0080GA27.9% 300.5 / 134.8 165intergenic (‑184/‑215)acrS/acrEacrAB operon transcriptional repressor/cytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (63/56);  new base A (32/14);  total (95/70)
Fisher's exact test for biased strand distribution p-value = 5.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 168 total aligned reads displayed.

AGTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTATTTTT  >  CP009273/3406861‑3407155
                                                                                                                                                       |                                                                                                                                                   
agTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACCTCAATATATATTTCAATTTACGAGGtt                                                                                                                                                   <  2:264249/150‑1 (MQ=255)
 gTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGttt                                                                                                                                                  >  2:420848/1‑150 (MQ=255)
 gTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTAAATCAATATATATTTCAATTTACGAGGttt                                                                                                                                                  >  2:4324185/1‑150 (MQ=255)
   gagaAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTCTTTACCTAAGATACAGTCACTACATCAATATATATTTCAAGTTACGAAGTTTTa                                                                                                                                                >  2:1599478/1‑150 (MQ=255)
    agaAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTaa                                                                                                                                               >  1:11452/1‑150 (MQ=255)
    agaAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTaa                                                                                                                                               <  1:3607682/150‑1 (MQ=255)
    agaAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTaa                                                                                                                                               <  1:3607689/150‑1 (MQ=255)
       acaTAGCGAAGGTTAATCTATCACCTAATGAGTCTTTAGACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTc                                                                                                                                            <  2:1141058/148‑1 (MQ=255)
       aaaTAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTc                                                                                                                                            <  1:666706/150‑1 (MQ=255)
       aaaTAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTc                                                                                                                                            <  2:1141221/150‑1 (MQ=255)
         aTAGCGAAGGTTAATCTATCACCTACTGTGTATTTCTACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTg                                                                                                                                          <  2:4313954/150‑1 (MQ=255)
           aGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGcc                                                                                                                                        >  1:1857019/1‑150 (MQ=255)
            gCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCct                                                                                                                                       >  2:2008345/1‑150 (MQ=255)
             cGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCctc                                                                                                                                      >  2:1849947/1‑150 (MQ=255)
             cGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCctc                                                                                                                                      >  1:2889275/1‑150 (MQ=255)
              gaattttAATCTATCACATAATTTGTATTTATACGAGAGTATAATATTGAGTTGCTATAAATATT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCctct                                                                                                                                     <  2:3728907/145‑1 (MQ=255)
              gAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCctct                                                                                                                                     >  2:3902184/1‑150 (MQ=255)
              gAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCctct                                                                                                                                     <  2:4096026/150‑1 (MQ=255)
                 ggTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTc                                                                                                                                  <  1:3438510/150‑1 (MQ=255)
                     aaTCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAAcc                                                                                                                              >  1:2034649/1‑150 (MQ=255)
                       tCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCcg                                                                                                                            <  1:1689573/150‑1 (MQ=255)
                           tCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCTAACCGCGTc                                                                                                                        >  2:4024626/1‑150 (MQ=255)
                           tCACCTAATGTGTATTAATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCcg                                                                                                                            <  2:2311680/150‑1 (MQ=255)
                            cACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCa                                                                                                                       <  2:1352806/150‑1 (MQ=255)
                            cACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCa                                                                                                                       <  1:2353090/150‑1 (MQ=255)
                             aCCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCaa                                                                                                                      >  2:2370313/1‑150 (MQ=255)
                             aCCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCaa                                                                                                                      >  1:1341987/1‑150 (MQ=255)
                                tAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAAt                                                                                                                   <  2:697810/150‑1 (MQ=255)
                                  aTGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTACTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATaa                                                                                                                 >  1:2393941/1‑150 (MQ=255)
                                   tgtgTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATaaa                                                                                                                >  2:898370/1‑150 (MQ=255)
                                         tttATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTa                                                                                                          <  1:1303352/150‑1 (MQ=255)
                                         tttATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTa                                                                                                          <  2:2872555/150‑1 (MQ=255)
                                         tttATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAAGATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTa                                                                                                          <  2:2550181/150‑1 (MQ=255)
                                         tttAGACGAGAGGCTAATATTGAGGTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTa                                                                                                          <  2:2550138/150‑1 (MQ=255)
                                          ttATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAg                                                                                                         <  2:4892/150‑1 (MQ=255)
                                          ttATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAg                                                                                                         <  2:3969546/150‑1 (MQ=255)
                                            ataCGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGaa                                                                                                       >  1:2639305/1‑150 (MQ=255)
                                               cGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAatat                                                                                                    <  1:1668729/150‑1 (MQ=255)
                                                gagaGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATAtt                                                                                                   <  2:2950823/150‑1 (MQ=255)
                                                gagaGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATAtt                                                                                                   <  1:914133/150‑1 (MQ=255)
                                                 agagGCTAATATTGAGTTTCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTa                                                                                                  <  1:2912218/150‑1 (MQ=255)
                                                 agagGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTa                                                                                                  <  1:2349103/150‑1 (MQ=255)
                                                 agagGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTa                                                                                                  <  1:2349075/150‑1 (MQ=255)
                                                  gagGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTaa                                                                                                 >  1:1036778/1‑150 (MQ=255)
                                                  gagGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTaa                                                                                                 >  1:1036624/1‑150 (MQ=255)
                                                   agGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAAt                                                                                                <  2:1568267/150‑1 (MQ=255)
                                                      cTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCtt                                                                                             <  2:3259074/150‑1 (MQ=255)
                                                      cTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATAGATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCtt                                                                                             >  1:1795358/1‑150 (MQ=255)
                                                         atatTGAGTTGCTATAAATAGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCttttt                                                                                          <  2:4210240/150‑1 (MQ=255)
                                                           atTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCttttttt                                                                                        >  1:459093/1‑150 (MQ=255)
                                                           atTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCttttttt                                                                                        >  1:2632283/1‑150 (MQ=255)
                                                             tGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTtgt                                                                                      >  1:2251701/1‑150 (MQ=255)
                                                              gAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTtgtg                                                                                     >  2:1074068/1‑150 (MQ=255)
                                                               aGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTgtgt                                                                                    <  2:4251667/150‑1 (MQ=255)
                                                                gTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGtt                                                                                   >  1:2103351/1‑150 (MQ=255)
                                                                   gCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAt                                                                                >  1:3195183/1‑150 (MQ=255)
                                                                    cTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtg                                                                               >  1:2225438/1‑150 (MQ=255)
                                                                    cTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtg                                                                               >  1:2225312/1‑150 (MQ=255)
                                                                    cTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTAtg                                                                               >  2:219013/1‑150 (MQ=255)
                                                                         aaTCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCt                                                                          <  2:125058/150‑1 (MQ=255)
                                                                         aaTCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCt                                                                          >  2:4329077/1‑150 (MQ=255)
                                                                         aaTCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCt                                                                          >  2:972120/1‑150 (MQ=255)
                                                                             gT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTgag                                                                      <  1:266547/150‑1 (MQ=255)
                                                                                    aaataaatAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATg                                                                   <  2:1312294/150‑1 (MQ=255)
                                                                                           taataTATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTAt                                                            >  1:3428219/1‑150 (MQ=255)
                                                                                           taataTATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTAt                                                            >  1:2566272/1‑150 (MQ=255)
                                                                                                tatatTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATa                                                       <  2:2034649/150‑1 (MQ=255)
                                                                                                tatatTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATa                                                       <  2:3703248/150‑1 (MQ=255)
                                                                                                  tattatTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGATTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAAc                                                     <  1:193770/150‑1 (MQ=255)
                                                                                                   attattTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACt                                                    >  2:531134/1‑150 (MQ=255)
                                                                                                   attattTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACt                                                    >  2:530278/1‑150 (MQ=255)
                                                                                                     tattTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTAt                                                  >  1:4189769/1‑150 (MQ=255)
                                                                                                     tattTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTAt                                                  >  1:4310603/1‑150 (MQ=255)
                                                                                                         tACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTCGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAAATATTCCt                                              >  1:2404145/1‑150 (MQ=255)
                                                                                                         tACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCt                                              >  1:3101793/1‑150 (MQ=255)
                                                                                                                gATACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcg                                       <  1:2947546/150‑1 (MQ=255)
                                                                                                                  tACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcgac                                     >  1:3696777/1‑150 (MQ=255)
                                                                                                                  tACATTCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcgac                                     <  1:467981/150‑1 (MQ=255)
                                                                                                                   aCATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATcgacg                                    >  1:3508895/1‑150 (MQ=255)
                                                                                                                   aCATTCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCATTACATcgacg                                    >  2:2936654/1‑150 (MQ=255)
                                                                                                                     aTTCACTACATCAATATATATTTCAATTTAAGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGaa                                  <  2:1479171/150‑1 (MQ=255)
                                                                                                                       tCACTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATg                                >  1:811891/1‑150 (MQ=255)
                                                                                                                       tCACTACATCAATATATATTTCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATg                                >  1:2524409/1‑150 (MQ=255)
                                                                                                                          cTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATa                             <  2:3425115/150‑1 (MQ=255)
                                                                                                                          cTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATa                             <  2:3425246/150‑1 (MQ=255)
                                                                                                                          cTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATa                             >  1:2811200/1‑150 (MQ=255)
                                                                                                                          cTACATCAATATATATTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATa                             <  2:1856488/150‑1 (MQ=255)
                                                                                                                              aTCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAAttt                         >  1:4252408/1‑150 (MQ=255)
                                                                                                                              aTCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAAttt                         >  1:87053/1‑150 (MQ=255)
                                                                                                                                cAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGt                       <  2:2307098/150‑1 (MQ=255)
                                                                                                                                       tatTTCAATTTACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAg                >  1:3383285/1‑150 (MQ=255)
                                                                                                                                         tttCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCg              <  2:1686404/150‑1 (MQ=255)
                                                                                                                                          ttCAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGa             >  2:3064657/1‑150 (MQ=255)
                                                                                                                                            cAATTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAAc           >  2:684731/1‑150 (MQ=255)
                                                                                                                                              aTTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTg         >  2:1957190/1‑150 (MQ=255)
                                                                                                                                              aTTTACGAAGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTACTTACATCGACGAATGATAATTTGTAGGATAGCGAACTg         >  1:1517699/1‑150 (MQ=255)
                                                                                                                                                ttACGAGATTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATGTGTAGGATAGCGAACTGTa       >  1:2689968/1‑150 (MQ=255)
                                                                                                                                                  aCGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAtt     <  1:184312/150‑1 (MQ=255)
                                                                                                                                                     aagTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAttttt  <  1:298787/148‑1 (MQ=255)
                                                                                                                                                     aagTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTAttttt  <  1:1174138/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
AGTGAGAAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGT‑‑‑‑TAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATATATTTCAATTTACGAGGTTTTAATTCTGCCTCTTTCAACCCGCGTCAAAATAAAACAGTAGAATATTAATCTTTTTTTGTGTTTATGTGCCTTGAGATGCCTGTATTCATAACTATTCCTTACATCGACGAATGATAATTTGTAGGATAGCGAACTGTATTTTT  >  CP009273/3406861‑3407155

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: