Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,943 T→G L416L (CTT→CTG)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,9430TG98.7% 217.3 / ‑5.2 74L416L (CTT→CTG) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (42/31);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑ACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC  >  NC_000913/4637699‑4638189
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGGGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGt                                                                                                                                                                                                                                                  >  1:386117/1‑251 (MQ=255)
     gTCAGCGAAGCGCGCTCCGTGCAGTGAGAAGAAAAAAACAGCACGCTACAGGTTAAGCCCACGGAGGTGAATGTTGCAGCTGAACCGGCGTTACCGTAGCAGGAGCTGGGGAATTTACTTGATAACGCCATCGAATTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGATGGATAGCGGTAGTGGTATTCCTGACATCGCGCTGTCACTTAtttt                                                                                                                                                                                                                                             <  2:69701/251‑1 (MQ=255)
      tCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGATATGCCCACTGAGGTTAATCTTGCTGCTGAACCGGCGTTACTGGAGCAGACGCTGCGGAATTTACTGGATAACGCAATAGATTTTACCCCCGAAAGCGGTCGCATCACGCTAAGCGCAGAAGGGGATCAGGAACAAGTC‑CCCCTTAAGGTGCAGGAAACCGGTAGGGGTTTTCATGACTACGCGATGTCACGAAttttt                                                                                                                                                                                                                                            >  1:668988/1‑251 (MQ=255)
             aGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCt                                                                                                                                                                                                                                     <  1:10495/251‑1 (MQ=255)
                 cgcACCGTGCAGTTGGCAGAAAAAAACAGCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGGCGCATAACGCTAAGCGCCGAAGTGGATCAGGATCACGTC‑GCCCTTAAGGTGCTGGATTCCGGTAGTGGTATTCCTGACTGCGCGCTGTCATGTATTTTTGAGCGTTTTTa                                                                                                                                                                                                                                 <  2:318241/251‑1 (MQ=255)
                    cccGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTActc                                                                                                                                                                                                                              <  1:384307/250‑1 (MQ=255)
                       gTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTTCTCttt                                                                                                                                                                                                                           <  2:165030/251‑1 (MQ=255)
                          cAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAAGCTAAGCGCCGAAGTGGATCAGGCACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTATTCTTTGcc                                                                                                                                                                                                                        >  2:7121/1‑251 (MQ=255)
                           aGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGccg                                                                                                                                                                                                                       <  1:60501/251‑2 (MQ=255)
                                gCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATACCGCCATCGATTTTACCCCCGAGAGCGGTCGCCTAACGCTAAGCGCCGAAGTAGATCAGGAACACGTC‑GCCCTTCAGGTGCCGGCTACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGccgcggtc                                                                                                                                                                                                                  >  2:86737/1‑245 (MQ=255)
                                 cAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGc                                                                                                                                                                                                                  >  1:535481/1‑250 (MQ=255)
                                 cAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGATGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGAC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGGATTCCTGCCTACGCGCTGTCACGTCTTGTTGAACGCTTTTACTCTTTGCCGCGTGc                                                                                                                                                                                                                  >  2:598509/1‑250 (MQ=255)
                                  aGAAAAAAACATCACTCTACATGTTATGCACACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGCTCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGAGGTATTCCTGACTACGCGCTGTCCCGTATTTTTGAACGCTTTTAATCTTTGCCGCGTGcc                                                                                                                                                                                                                 >  1:433202/1‑249 (MQ=255)
                                    aaaaaaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAAt                                                                                                                                                                                                              <  1:628074/251‑1 (MQ=255)
                                        aaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGc                                                                                                                                                                                                          >  1:281782/5‑251 (MQ=255)
                                             tCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaa                                                                                                                                                                                                      >  1:547763/1‑250 (MQ=255)
                                                ctctACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTCGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGAAAAAAagc                                                                                                                                                                                                  >  1:655592/1‑251 (MQ=255)
                                                  ctacATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAagc                                                                                                                                                                                                  >  2:161893/4‑249 (MQ=255)
                                                        gTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCt                                                                                                                                                                                          >  1:263633/1‑251 (MQ=255)
                                                             gCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGt                                                                                                                                                                                     >  1:562475/1‑251 (MQ=255)
                                                                   tgAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGCTAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGCAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGCTACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTATTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGTCGt                                                                                                                                                                               >  2:409127/2‑251 (MQ=255)
                                                                     aGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc                                                                                                                                                                             <  1:3515/251‑1 (MQ=255)
                                                                     aGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc                                                                                                                                                                             <  2:281782/251‑1 (MQ=255)
                                                                            tgTTGCTGCTGAACCCGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGCTTTTACCACCGAGAGCGGTGGCATAACGCTAAGCGACGAAGTGGCTCAGGGACACGTC‑GCCCTTAAGGTGCTGGCTACCGGTAGTGGAATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGGGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGTGTTAGTCAGTg                                                                                                                                                                      >  2:264863/2‑251 (MQ=255)
                                                                               tgctgcTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAgg                                                                                                                                                                   >  1:181873/1‑251 (MQ=255)
                                                                                gctgctGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATTCCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTTAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGt                                                                                                                                                                  <  1:264863/251‑1 (MQ=255)
                                                                                gctgctGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGGAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCAGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTccttggtggagg                                                                                                                                                                  >  2:587711/1‑240 (MQ=255)
                                                                                    ctgAACCGGCGTTACTGGAGCCGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTTCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGGTGGCGTTCGTCAGTGAGGTCGcc                                                                                                                                                              >  2:121530/1‑251 (MQ=255)
                                                                                      gAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCg                                                                                                                                                            <  1:428969/251‑1 (MQ=255)
                                                                                        aCCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCgtt                                                                                                                                                          <  1:161893/251‑1 (MQ=255)
                                                                                          cGGCGTTACTGGAGCAGACGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGCGCGGTCGCATAACGCTAAGCGCCGCAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGAGGTATTCCTGACTACGCGCTGTCCAGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAACTGGGCAAAATAGCAGCGGCCTGGGGTTGGCGTTCGTCAGTggtgtcgcacggctg                                                                                                                                                        >  2:644496/1‑237 (MQ=255)
                                                                                           ggCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCgtttgt                                                                                                                                                       <  1:7121/251‑1 (MQ=255)
                                                                                                    ggAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCAc                                                                                                                                                                                                                                                    <  1:490717/149‑1 (MQ=255)
                                                                                                    ggAGCAGGCACTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAAGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCAc                                                                                                                                                                                                                                                    >  2:490717/1‑149 (MQ=255)
                                                                                                       gCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGGCGCCCGTTTGTTTAACGGCGa                                                                                                                                             >  2:62712/1‑249 (MQ=255)
                                                                                                                    aaTTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCa                                                                                                                               >  1:524426/1‑250 (MQ=255)
                                                                                                                      tttACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAACTGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGACGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACg                                                                                                                            >  1:56344/1‑251 (MQ=255)
                                                                                                                       ttACTGGATATCGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGTACACGGC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGGCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGGCTGGGGTTGGCGTGCGTCAGTGAGGTCGCCCGTTTGTTTAACGTCGAAGTCACGCTGCGCAACg                                                                                                                            <  2:655592/250‑1 (MQ=255)
                                                                                                                        tACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGt                                                                                                                           <  1:184287/250‑1 (MQ=255)
                                                                                                                        tACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGCCGCATAACGCTAAGCGACGAAGTGGATCAGGAACACGTC‑GCCCTTATGGAGCTGGATACCGGTAGTGGTATTCCTGACTACGCGATGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGAGAGGCCGCCCGATTGTTTAACGGAGAAGTCACGCAGCGCaaaggg                                                                                                                          >  2:622637/1‑247 (MQ=255)
                                                                                                                            ggATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAgg                                                                                                                      <  1:459385/251‑1 (MQ=255)
                                                                                                                             gATAACGCCATCGATTTTACCCCCGAGAGCGGACGCATAACGCTAAGCGCCGAAGTGGATCAGGAACCCGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCCCGTATTTTTGAACGCTTTTACTCTTTGCCGCGGGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGGCAGTGAGGTCGCCCGTTTGTTTAACGGCGAATACACGCTGCGCTACGTGCAgg                                                                                                                      >  2:329964/1‑250 (MQ=255)
                                                                                                                               taactacATCGATTTTACACCCGAGAGCGGCCGCATAACGCTAAGCGCCGAAATTGATCAGGACCACGTCGAC‑CTTAAGGTGCTGGTTACCGGAAGTGGTATTCCTGACTAAGCGCTGACACGTATTTTTGAACGAATTTCCTCTTTGCCGCGTGCACATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGCCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTAACGCTGCGCAACGTGCAGGa                                                                                                                     >  2:586035/7‑249 (MQ=255)
                                                                                                                                      aTCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaaaa                                                                                                                                                                                                    <  1:423594/163‑1 (MQ=255)
                                                                                                                                      aTCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaaaa                                                                                                                                                                                                    >  2:423594/1‑163 (MQ=255)
                                                                                                                                               aCCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGCCTACGCGCTGTCACGTATTTTTGAAAGCGTTTACTCTTTGTCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGAGAGGGCGCCCGTTTTTTTAACGGCGAAGGCACGCGGCGCGCAGTGCAGGAAGGGGGCGTGCTGGGCTc                                                                                                   >  2:558172/1‑251 (MQ=255)
                                                                                                                                                cccccGATAGCGGTCGCATAAAGCTAAGCGCCGAAGTGAATCAGGCACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCATGACTACGCGCTGTCACGTATTTATGAACGCTTTTACTCTTTTCCGCGTGCAAATTGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAAAGGAGGATCCACGCTGCGAAACGTACAGGAAGGGGGCGTgatgtccgcg                                                                                                  >  2:182666/1‑242 (MQ=255)
                                                                                                                                                cccccGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCATGCTGGCCTCg                                                                                                  >  1:352120/1‑251 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCCGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAg                                                                                                                                                                        >  2:465871/1‑171 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAg                                                                                                                                                                        <  1:465871/171‑1 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTc                                                                                        >  1:29651/1‑251 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGAGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTACGTGAGGTCCCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGAGCAGGCCTCGCTTCGTCttt                                                                                        >  2:527712/1‑250 (MQ=255)
                                                                                                                                                                    aCGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCg                                                                                             <  1:405711/236‑1 (MQ=255)
                                                                                                                                                                    aCGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGGGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCg                                                                                             >  2:405711/1‑236 (MQ=255)
                                                                                                                                                                        tAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCACATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCACCGTGCAGGAGGATGGCGTGCTGGCCTCGCTTCGACTTCACCGTCCCTTccc                                                                           >  2:220653/1‑248 (MQ=255)
                                                                                                                                                                               cGAAGTGGATCAGGAACACGTC‑ACGCGTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCACGAAGGTGGCGTGCTTGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTc                                                                   <  1:69898/251‑1 (MQ=255)
                                                                                                                                                                                          aGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTccc                                                        <  1:14315/251‑1 (MQ=255)
                                                                                                                                                                                           gtgacacGTC‑GCCCTTAAGGTGCTGTATACCGGTAGTGGTATTCCTGATTACGCGCTGTCACGTATTTGTGAACGCTTTTTCGCTTTGCCGCGTGGAAATGGGCAAAAAAGCATCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGGGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCca                                                       <  2:545939/248‑1 (MQ=255)
                                                                                                                                                                                                 acgtcgcccTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGAGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGAGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCAACAt                                                    >  2:487096/9‑249 (MQ=255)
                                                                                                                                                                                                  cgtcgcccTTAAGGTGCTGGATCCCGGTAGTGGGATTCCTGACTACGCGCTGTCACGCATTTTTGACCGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGGGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCCCATCGCTTCAAATTCTTCCCACATAgc                                                 >  2:576003/8‑250 (MQ=255)
                                                                                                                                                                                                   gtcgcccTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGCGGCGGGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGATTCAAATTCTTCCCACATAgcc                                                >  2:571366/7‑249 (MQ=255)
                                                                                                                                                                                                     cgcccTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATaacctt                                              >  2:239199/5‑246 (MQ=255)
                                                                                                                                                                                                        ccTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGt                                           >  1:418704/2‑251 (MQ=255)
                                                                                                                                                                                                           tAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCTCGTATTTTTGAACGCTTTTACTCGTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTGCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCTGCACTTCACATAGCTGCAAATTCTTCCCACATAGCCGTCGTa                                          <  1:225508/249‑1 (MQ=255)
                                                                                                                                                                                                            aaGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATc                                        <  1:86737/250‑1 (MQ=255)
                                                                                                                                                                                                             aggttccggATACCGGTAGTGTGATTCCTGACTACGCGCTGTCACGTAGTTTTGAACGTTTTTACTCTTGGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTTCTGTCCTCGCTTCGACTTCACCTTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTTATCct                                      <  2:119926/244‑1 (MQ=255)
                                                                                                                                                                                                                    ggATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGcct                               >  1:183997/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCAtt                             <  1:239199/251‑1 (MQ=255)
                                                                                                                                                                                                                              aGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGga                     <  1:571366/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGCGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCcaca                                                     >  2:183393/1‑204 (MQ=255)
                                                                                                                                                                                                                                             tACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGGCACTTCACATAGCTTCAAATTCTTCCcaca                                                     <  1:183393/204‑1 (MQ=255)
                                                                                                                                                                                                                                              aCGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGaaa     <  1:69019/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgcTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTTGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATcc  <  1:31996/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgcTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATcc  <  1:121530/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑ACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC  >  NC_000913/4637699‑4638189

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑