Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,484,827 C→T L152F (CTT→TTT)  evgS → sensor histidine kinase EvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,484,8270CT97.6% 127.4 / NA 42L152F (CTT→TTT) evgSsensor histidine kinase EvgS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/24);  minor base A (1/0);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTAAATGT  >  NC_000913/2484592‑2485062
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           
cGATTCGCAGCAACGGGGTAGTGGTATTAAATCTGATTATTTAAAGCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATAAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGc                                                                                                                                                                                                                              <  2:201018/251‑1 (MQ=255)
            aggggatcctgatattcatgcaaattattTTAATCATTTAAAAAGAGCGTTAAATATCAAATTAACACTCTGTGAATACGCTGATCATCAACAAAAAATGTATCAC‑CTTGAAGATGCTCGAGTAGATATA‑GAGTTATCTCATTTACTACCTTCGC‑GACTCCTTACTAATGCCATTGCTCCAACCAAACCATTGATAATTCCATTTTTGTCGCTGGCACCCACCATTCACGCCTAAATGCGACCGCTGACCt                                                                                                                                                                                                                  >  2:525261/24‑251 (MQ=255)
            aCGGGTTCATGGTATTAATGCTGATTATTTAAATCTGTTAAAAATAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAATCAATGGA‑CGCGCCTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCGCCTTTCACGACTCAATTCGACCGCTTACCt                                                                                                                                                                                                                  <  2:118519/251‑1 (MQ=255)
                   cGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGAGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTAACCaa                                                                                                                                                                                                            >  2:688256/1‑250 (MQ=255)
                                        ttAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGt                                                                                                                                                                                      >  1:76526/1‑251 (MQ=255)
                                         tAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGt                                                                                                                                                                                      >  2:291036/1‑250 (MQ=255)
                                                 ttAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAAtt                                                                                                                                                                              >  1:332873/1‑250 (MQ=255)
                                                   aaaaaGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcc                                                                                                                                                                           >  2:14547/1‑251 (MQ=255)
                                                      aaGAGCGTTAAATATCAAATTAACACTCCGTGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATGTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccc                                                                                                                                                                        <  1:382262/251‑1 (MQ=255)
                                                           cGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGAc                                                                                                                                                                   <  1:14547/251‑1 (MQ=255)
                                                               aaaTATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAgg                                                                                                                                                               >  1:611160/1‑251 (MQ=255)
                                                                      aaaTTAACACTCCGGGAATACGAAAATCATCAACAAGAAATGGA‑CGCGCTTGCAGAAGGTGAAGCCGATATA‑GTGTTATCACATTTAGTTACTTCGACACATC‑TTAATAATGAAATTACTGAAACCAAACCATTGATAATTACCTTTCCGACGCTGGTAACCACCTTTCACGACTAAATGCGACCGCTTACCTCAACAACACCAGTAAATATTGCTCGGGTAGCAAATTACCACACAGACGAGGAAATtc                                                                                                                                                         >  2:779484/1‑250 (MQ=255)
                                                                       aaTTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATtcat                                                                                                                                                       >  1:487236/1‑251 (MQ=255)
                                                                         ttAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCATTTCTACCGCTTACATCACCAAAACCAGTAAATATTGCTCGGGAAGCAACTTACCCCCCAGACGAGGAAATtcatc                                                                                                                                                      >  2:786118/1‑250 (MQ=255)
                                                                                   gggAATACGCAGAACAGCAAAAAGCAATGGT‑CGCGCTTGCAGGAGGTGAAGTCTATATA‑GTGTTATCACATTTAGTTACGTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATGGCCCCCCAGACGAGGTAAGTCATCAACCATTTCCa                                                                                                                                           <  2:446347/251‑1 (MQ=255)
                                                                                        aaCGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTCCTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGTGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                                                                      <  2:487236/250‑1 (MQ=255)
                                                                                           gCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaac                                                                                                                                   <  2:407744/251‑1 (MQ=255)
                                                                                             agatctgcaAAAAGCAAGGGA‑CGCGCTTGCAGACGTTGAATTCGATATC‑ATTTTATCACAGTTAGTAACTTCTCCGCCGC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAAGTACCTTTCCGGCGCTGTTAACCCCCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAAACATTAACTATTTCTCTGGTAGCAAATTCCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaacaa                                                                                                                                 <  2:4306/244‑1 (MQ=255)
                                                                                               atcatcAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAtt                                                                                                                               <  1:768875/251‑1 (MQ=255)
                                                                                                 catcaAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATGGCTGCAACCAAACCATTGATAATTTCCTTTCCGGCGCTGGTAACCCCCTTTCACGACTCAATGCGACCGCTTACCTCCCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt                                                                                                                             <  1:550693/251‑1 (MQ=255)
                                                                                                           cAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACa                                                                                                                   <  1:816375/251‑1 (MQ=255)
                                                                                                                    gctctggCAGATGGTGAAGTCGATATAAGT‑TTATCACATTTATTTACTTCACCGCCTC‑TTAATAATGACATTTCTGCAACCAAACCATTGATAATTACCTTTCCGCCGCTCGTAACCACCTTTCACTACTCAATGCGACCGCTTACCTCTCCAAAACCAGTAAATATTGCTCGGGTAGCAAATAACCCCCCAGACGATGTAATTCATCAATCATGTCCAAAAGCAACAATTATCTCTTTTACAAATTtata                                                                                                           <  2:6093/245‑1 (MQ=255)
                                                                                                                          gCAGAAGGAGAAGTCGGTATA‑GTGTTATCTCATTTAGTTACTTCGCCGCCTA‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTAGTAACCACCTGTCACGACTCAGTGCTACCGATTACCTAAACAAAACCAATACATATTGCGCGGGTAGCAAATTACACCCCATACGAGGTAGTTCATCAATCATTACCAACAGCAACAATTATCTCTTTTCCAAATTTATATCAGGc                                                                                                     >  2:132973/1‑251 (MQ=255)
                                                                                                                              aaGGTGAAGTCGATATA‑TTTTTATCACATTTAGTTTCTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCattcatta                                                                                                 <  2:392007/251‑8 (MQ=255)
                                                                                                                                ggTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGc                                                                                               <  1:5929/251‑1 (MQ=255)
                                                                                                                                        cGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTa                                                                                                                                                                                     >  2:389502/1‑157 (MQ=255)
                                                                                                                                              a‑GTGTTATCACATTTAATTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCt                                                                                  <  2:872342/250‑1 (MQ=255)
                                                                                                                                                        acaTTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATg                                                                        >  1:344659/1‑251 (MQ=255)
                                                                                                                                                                           ccTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCg                                                                                          >  2:470417/1‑214 (MQ=255)
                                                                                                                                                                           ccTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCg                                                                                          <  1:470417/214‑1 (MQ=255)
                                                                                                                                                                                   ataatGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAcat                                                  <  1:115170/247‑1 (MQ=255)
                                                                                                                                                                                   ataatGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAcat                                                  >  2:115170/1‑247 (MQ=255)
                                                                                                                                                                                       tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTAcc                                           >  1:458229/1‑250 (MQ=255)
                                                                                                                                                                                            ttGCTGCAACCAACCCATTGATAATTACCTTTCCGGCGCTGGAAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaacaa                                                                                                                                 >  1:270129/1‑159 (MQ=255)
                                                                                                                                                                                            ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAAGTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaacaa                                                                                                                                 <  2:270129/159‑1 (MQ=255)
                                                                                                                                                                                             tgctgcAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGt                                     <  1:231136/250‑1 (MQ=255)
                                                                                                                                                                                                      caaaccaTTGATAATTTCCTTTCCGGCGCTGGTTACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTccc                           <  2:714121/251‑1 (MQ=255)
                                                                                                                                                                                                           caTTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTAt                      <  2:668972/251‑1 (MQ=255)
                                                                                                                                                                                                                           tCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAACATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTa       <  2:540701/250‑1 (MQ=255)
                                                                                                                                                                                                                               gcgcTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTAAATGt  <  1:688256/251‑1 (MQ=255)
                                                                                                                                                                                                                                          caccTTTCACGGCTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCAATTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTa       <  2:599921/235‑1 (MQ=255)
                                                                                                                                                                                                                                          caccTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTa       >  1:599921/1‑235 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           
CGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGA‑CGCGCTTGCAGAAGGTGAAGTCGATATA‑GTGTTATCACATTTAGTTACTTCGCCGCCTC‑TTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTAAATGT  >  NC_000913/2484592‑2485062

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: ‑