Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 553,161 G→A intergenic (‑61/+57) purE ← / ← lpxH N(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913553,1610GA98.5% 217.9 / ‑4.5 66intergenic (‑61/+57)purE/lpxHN(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/35);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

ATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCAC  >  NC_000913/552912‑553399
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
aTAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGTTGCACAAATCTCGGCTTTAGAGGGCACAGAGACTAGCACgg                                                                                                                                                                                                                                                >  2:28044/1‑250 (MQ=255)
         ctctTCGGCGCTTTCGGCGAGGCTGAACAGTTTATCGGGGGCGCGGTGCGCAGAAACCGCTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGACGGCGACACGCGCCGGATTATTGCGGGACGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGCGAGc                                                                                                                                                                                                                                       >  2:98599/1‑250 (MQ=255)
         ctctTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGc                                                                                                                                                                                                                                       >  1:560717/1‑250 (MQ=255)
          tctTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCa                                                                                                                                                                                                                                      >  1:164916/1‑250 (MQ=255)
               ggCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGaaa                                                                                                                                                                                                                                <  2:436469/251‑1 (MQ=255)
                 cgcTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAAc                                                                                                                                                                                                                              >  1:349254/1‑251 (MQ=255)
                     ttCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCATCGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTTCATGGTAGCCCATTCGCTTTTGGACCCCATCACGATGTCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTTGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGt                                                                                                                                                                                                                           <  1:382369/250‑1 (MQ=255)
                       gggCGAAGCTGGACAGTTTAACGGGGGTGCGGTGAGTGGAAACCACTTCAACGTGGTTCGGGGGAGTGAGGATTTCGAAGATTTCGGCGGCAAACTGCATGGTAGCTCAGTCGCTTGCGGACTCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGGCATGCGTCTTAAAACTCCGGGGGTGCACAACTCGCGGCTTTAGAGGGCACAGATACAAGCACGAAAAGAGAGCAAGGAAAACGGTTGc                                                                                                                                                                                                                        <  2:439015/250‑1 (MQ=255)
                                  aaCAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGaa                                                                                                                                                                                                             <  1:311183/251‑1 (MQ=255)
                                   aCAGTGTATCGGGGGTGCGGTGCGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCGTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGaaa                                                                                                                                                                                                            <  2:268027/251‑1 (MQ=255)
                                         tATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGc                                                                                                                                                                                                       >  1:333947/1‑250 (MQ=255)
                                            cGGGGGTGCGGTGAGCAGAAACCACTGCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCaaa                                                                                                                                                                                                    <  2:402363/250‑1 (MQ=255)
                                             gggggTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGt                                                                                                                                                                                                  >  2:482384/1‑251 (MQ=255)
                                             gggggTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGt                                                                                                                                                                                                  <  2:315885/251‑1 (MQ=255)
                                              ggggTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGAAAAGtt                                                                                                                                                                                                 >  2:424278/1‑251 (MQ=255)
                                              ggggTGCGGGGAGCAGAAACCACTTCAACGTGGTTCGTGACATTCAGGATTTGGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGAACCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGtt                                                                                                                                                                                                 <  1:479103/251‑1 (MQ=255)
                                                 gTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCg                                                                                                                                                                                              >  1:510434/1‑251 (MQ=255)
                                                 gTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCg                                                                                                                                                                                              <  2:52594/251‑1 (MQ=255)
                                                     ggggAGCAGAAACCACCTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGtt                                                                                                                                                                                          <  2:360731/248‑1 (MQ=255)
                                                      gagAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTAGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGAGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAACAGGGAGCAAGGAAAACGGTTGCGTGGCTGGGAAATCAGCAAAGTTGGGGGttt                                                                                                                                                                                         >  2:579002/3‑251 (MQ=255)
                                                        gAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGtttt                                                                                                                                                                                        <  1:44925/250‑1 (MQ=255)
                                                         aGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGtttttt                                                                                                                                                                                      <  1:514241/251‑1 (MQ=255)
                                                          gCAGAAACCACTTCAACGTGGTGCGGGCCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACGCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTa                                                                                                                                                                                     <  2:384947/251‑1 (MQ=255)
                                                                      tCAACGTGGTGCTGTACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAAt                                                                                                                                                                          <  1:206902/250‑1 (MQ=255)
                                                                      tCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGaaaa                                                                                                                                                                           <  1:141636/249‑1 (MQ=255)
                                                                       cAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATg                                                                                                                                                                         <  1:531543/250‑1 (MQ=255)
                                                                       cAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGGTTTTTAAAACGGAAAATg                                                                                                                                                                         >  2:427056/1‑250 (MQ=255)
                                                                               tgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTc                                                                                                                                                                >  1:109691/1‑251 (MQ=255)
                                                                                    gACATTAAGGATATCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGAGCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAAGGACTAGCACGAAAAGAGAGCAAGGCAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGt                                                                                                                                                           >  2:278008/1‑251 (MQ=255)
                                                                                       agtCAGGATTTCGAAGATTGCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCATGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCAt                                                                                                                                                        <  1:506669/249‑1 (MQ=255)
                                                                                                ttCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTg                                                                                                                                                >  1:357115/1‑250 (MQ=255)
                                                                                                           cggcggCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCAtt                                                                                                                                    >  2:446432/1‑251 (MQ=255)
                                                                                                            gggggggAACTGCATGGTAGCCCAGACGTTTTTGGACCCCATCACGATGGCGACCCGGGCCGGATTTTTGGGGGAAAACAGGCGGCTTAAAACCCCGGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAACCGGTGGCGTGGCTGGGAAAGCAGCAAATTTGCGGGTTTTTTAAAACGGAAAATGAAACAGCTCAACGTCATCCGCCGTGACTTTCACCATTg                                                                                                                                   <  2:346846/245‑1 (MQ=255)
                                                                                                                    aCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGCACCTTcc                                                                                                                           >  1:340102/1‑251 (MQ=255)
                                                                                                                     cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCg                                                                                                                          >  1:473855/1‑251 (MQ=255)
                                                                                                                        cATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCAAAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTAt                                                                                                                       >  1:357888/1‑251 (MQ=255)
                                                                                                                           ggTAGCCCAGTCGCTTTTGGACCCCATCCCGAGGGCTACCGGCGACGAATTTTTGGGGGAAGACATGCGTCTAAAAACTCCTGGGGTGCACAACTCTCGGCTTGAGAGGGCAAAGAGAATAGAACGAAAAGAGAGCAAGGAAAACGGTGGCGGGGCTGTGAAATCAGCAAAGTGGCGGTTTTTTTAAAACGGAAAATGATCCATCTCACCGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGcc                                                                                                                    <  2:125829/251‑1 (MQ=255)
                                                                                                                           ggTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGCCATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGcc                                                                                                                    >  2:161337/1‑251 (MQ=255)
                                                                                                                            gTAGCCCAGGCGCTTTTGGACCCCAACAGGATGGCGACACGCGCCGGATTATTGCGGAAAGACATGCGCCTTAAAACTCCTGTGGGGCACAACTCTCGGCTTTAGAGGGCACAGGAAATAGCACGAAAAGAGAGCAAGAAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGATTCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCa                                                                                                                   <  1:446432/251‑1 (MQ=255)
                                                                                                                                      cGCTTTTGGACCCCATCTCGATGGCCACATGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACCTCTGTGGTGCAAAATTCGCGGTTTTAGAGGTCACAGAGAATAGCACGAAAAGAGAGAAAGGAAACCGTTTGCGTGGCGTTGAAATCAGTAATGTGGCGGGTTTTTTAAAACGGAAAATGAATCAGTTAAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAg                                                                                                          <  1:508635/250‑1 (MQ=255)
                                                                                                                                            tGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCAc                                                                                                    <  1:482384/250‑1 (MQ=255)
                                                                                                                                              aaCCCCATCCCGATGGCGACACGCGCCGTATTATTGCGGGAAGAAATGCGCCTTAAAACCCCTGTGGTGCACACCTCGCGGCTTTAGAGGGCACAGAGAATAGCAGGAAAAGAGAGCAAGGAAAACGGTTGCGTGGTTGTGAAAGCAGCAAAGTTGGGGGTTTTTTAAACCGGAAAATGAATCAGCTCAACGCCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCAcgcg                                                                                                 <  2:228914/250‑1 (MQ=255)
                                                                                                                                                             gCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGGGGTGCACAACTCTCGGTTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTg                                                                                  <  1:487997/251‑1 (MQ=255)
                                                                                                                                                               gACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCGCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGAt                                                                                <  1:490963/251‑1 (MQ=255)
                                                                                                                                                                acacGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAAt                                                                                                                                                                          >  1:403647/1‑160 (MQ=255)
                                                                                                                                                                acacGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAAt                                                                                                                                                                          <  2:403647/160‑1 (MQ=255)
                                                                                                                                                                       ccGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAAAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTAAAAACGGAAAATGAATCAGCTCACCGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATa                                                                        <  1:391688/251‑1 (MQ=255)
                                                                                                                                                                                 gCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAg                                                                                                                                                                                                   <  1:430051/118‑1 (MQ=255)
                                                                                                                                                                                 gCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAg                                                                                                                                                                                                   >  2:430051/1‑118 (MQ=255)
                                                                                                                                                                                   gggAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACAGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCAcc                                                            >  2:51071/1‑251 (MQ=255)
                                                                                                                                                                                                 cTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCa                                                                                                                                                               <  1:119330/138‑1 (MQ=255)
                                                                                                                                                                                                 cTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCa                                                                                                                                                               >  2:119330/1‑138 (MQ=255)
                                                                                                                                                                                                    aaaaCTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGAGGGGTATGc                                          >  2:78992/1‑251 (MQ=255)
                                                                                                                                                                                                      aaCTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTg                                                                              <  1:74996/214‑1 (MQ=255)
                                                                                                                                                                                                      aaCTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTg                                                                              >  2:74996/1‑214 (MQ=255)
                                                                                                                                                                                                         tCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGt                                      >  1:311285/1‑251 (MQ=255)
                                                                                                                                                                                                          ccgggggggcacaccgttcGGTTTTAGAGGGAACAGAGAATAGCAAGAAAAGAGAGCAAGGAAAACGGTGGCGTGGCTGTGAAATCAGAAAAGTTGCGGGTTTTTTAAAACGGAAAATGAACCAGCAAAACGTCATCCGCCGTGACTTTCACCATTGACCCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTg                                     <  2:420200/234‑1 (MQ=255)
                                                                                                                                                                                                             ggggtgCACAACTCCCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTGGCGTGGCTGTGAAATCAGCAAAGTTGGGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGAt                                  <  1:177973/249‑1 (MQ=255)
                                                                                                                                                                                                                 ggaccaattctcGGCTTTAGGGGGCACAAAAAATAGCACAAAAAGAGAGAAAGGAAAACGGTTGCGGGGCTTTGAAATCAGCAAAGTTTGGGGTTTTTTAAAACGGAAAAGGAATCACCTCAGCCTCTTCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGc                              <  1:425877/243‑1 (MQ=255)
                                                                                                                                                                                                                  gcccacctcgcGGCTTTAGAGGGCAAAGCGAATAGATAGAAAAGAAAGCAAGGAAAACGGTTGCGTGGTTGTAAAATCAAAAAAGTGGCGGGTTTTTTAAAACGGAAAATGAATCAGCTAAACGTCTTCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGcc                             <  1:557371/241‑1 (MQ=255)
                                                                                                                                                                                                                      aaCTCTCGGCTTTAGAGGGCACAGAGAATAGCACAAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAACGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTg                         >  1:67933/1‑251 (MQ=255)
                                                                                                                                                                                                                        ctctCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGc                        >  2:446179/1‑250 (MQ=255)
                                                                                                                                                                                                                                  tAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCACTTGCACCTGATGttt             >  2:428193/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cacaGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAAACTTCCGTATGCCAGGCACCCAGTACCACGCGTAAAGCAGGATGCTGATTGGCGAAAAGGTACTGCACAGCCGGCAGATGGGAATGCCAGTGGGTCAGCAATTGCACACGCGGTTTTTCCATc      >  2:144751/1‑251 (MQ=255)
                                                                                                                                                                                                                                           caGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGCAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGc    >  2:94203/1‑251 (MQ=255)
                                                                                                                                                                                                                                            agagAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCga    >  2:417236/1‑249 (MQ=255)
                                                                                                                                                                                                                                             aagaATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGTTTTTTTAAAACAGAAAATGAATCAGCTCACCGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCAc  <  1:51071/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
ATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCAC  >  NC_000913/552912‑553399

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: