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→CTGcreCsensory 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

GCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC  >  NC_000913/4637699‑4638189
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGGGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGt                                                                                                                                                                                                                                                  >  1:386117/1‑251 (MQ=255)
     gTCAGCGAAGCGCGCTCCGTGCAGTGAGAAGAAAAAAACAGCACGCTACAGGTTAAGCCCACGGAGGTGAATGTTGCAGCTGAACCGGCGTTACCGTAGCAGGAGCTGGGGAATTTACTTGATAACGCCATCGAATTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGATGGATAGCGGTAGTGGTATTCCTGACATCGCGCTGTCACTTAtttt                                                                                                                                                                                                                                             <  2:69701/251‑1 (MQ=255)
      tCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGATATGCCCACTGAGGTTAATCTTGCTGCTGAACCGGCGTTACTGGAGCAGACGCTGCGGAATTTACTGGATAACGCAATAGATTTTACCCCCGAAAGCGGTCGCATCACGCTAAGCGCAGAAGGGGATCAGGAACAAGTCCCCCTTAAGGTGCAGGAAACCGGTAGGGGTTTTCATGACTACGCGATGTCACGAAttttt                                                                                                                                                                                                                                            >  1:668988/1‑251 (MQ=255)
             aGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCt                                                                                                                                                                                                                                     <  1:10495/251‑1 (MQ=255)
                 cgcACCGTGCAGTTGGCAGAAAAAAACAGCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGGCGCATAACGCTAAGCGCCGAAGTGGATCAGGATCACGTCGCCCTTAAGGTGCTGGATTCCGGTAGTGGTATTCCTGACTGCGCGCTGTCATGTATTTTTGAGCGTTTTTa                                                                                                                                                                                                                                 <  2:318241/251‑1 (MQ=255)
                    cccGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTActc                                                                                                                                                                                                                              <  1:384307/250‑1 (MQ=255)
                       gTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTTCTCttt                                                                                                                                                                                                                           <  2:165030/251‑1 (MQ=255)
                          cAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAAGCTAAGCGCCGAAGTGGATCAGGCACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTATTCTTTGcc                                                                                                                                                                                                                        >  2:7121/1‑251 (MQ=255)
                           aGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGccg                                                                                                                                                                                                                       <  1:60501/251‑2 (MQ=255)
                                gCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATACCGCCATCGATTTTACCCCCGAGAGCGGTCGCCTAACGCTAAGCGCCGAAGTAGATCAGGAACACGTCGCCCTTCAGGTGCCGGCTACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGccgcggtc                                                                                                                                                                                                                  >  2:86737/1‑245 (MQ=255)
                                 cAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGc                                                                                                                                                                                                                  >  1:535481/1‑250 (MQ=255)
                                 cAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGATGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGACGCCCTTAAGGTGCTGGATACCGGTAGTGGGATTCCTGCCTACGCGCTGTCACGTCTTGTTGAACGCTTTTACTCTTTGCCGCGTGc                                                                                                                                                                                                                  >  2:598509/1‑250 (MQ=255)
                                  aGAAAAAAACATCACTCTACATGTTATGCACACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGCTCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGAGGTATTCCTGACTACGCGCTGTCCCGTATTTTTGAACGCTTTTAATCTTTGCCGCGTGcc                                                                                                                                                                                                                 >  1:433202/1‑249 (MQ=255)
                                    aaaaaaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAAt                                                                                                                                                                                                              <  1:628074/251‑1 (MQ=255)
                                        aaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGc                                                                                                                                                                                                          >  1:281782/5‑251 (MQ=255)
                                             tCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaa                                                                                                                                                                                                      >  1:547763/1‑250 (MQ=255)
                                                ctctACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTCGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGAAAAAAagc                                                                                                                                                                                                  >  1:655592/1‑251 (MQ=255)
                                                  ctacATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAagc                                                                                                                                                                                                  >  2:161893/4‑249 (MQ=255)
                                                        gTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCt                                                                                                                                                                                          >  1:263633/1‑251 (MQ=255)
                                                             gCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGt                                                                                                                                                                                     >  1:562475/1‑251 (MQ=255)
                                                                   tgAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGCTAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGCAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGCTACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTATTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGTCGt                                                                                                                                                                               >  2:409127/2‑251 (MQ=255)
                                                                     aGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc                                                                                                                                                                             <  1:3515/251‑1 (MQ=255)
                                                                     aGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc                                                                                                                                                                             <  2:281782/251‑1 (MQ=255)
                                                                            tgTTGCTGCTGAACCCGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGCTTTTACCACCGAGAGCGGTGGCATAACGCTAAGCGACGAAGTGGCTCAGGGACACGTCGCCCTTAAGGTGCTGGCTACCGGTAGTGGAATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGGGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGTGTTAGTCAGTg                                                                                                                                                                      >  2:264863/2‑251 (MQ=255)
                                                                               tgctgcTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAgg                                                                                                                                                                   >  1:181873/1‑251 (MQ=255)
                                                                                gctgctGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATTCCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTTAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGt                                                                                                                                                                  <  1:264863/251‑1 (MQ=255)
                                                                                gctgctGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGGAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCAGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTccttggtggagg                                                                                                                                                                  >  2:587711/1‑240 (MQ=255)
                                                                                    ctgAACCGGCGTTACTGGAGCCGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTTCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGGTGGCGTTCGTCAGTGAGGTCGcc                                                                                                                                                              >  2:121530/1‑251 (MQ=255)
                                                                                      gAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCg                                                                                                                                                            <  1:428969/251‑1 (MQ=255)
                                                                                        aCCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCgtt                                                                                                                                                          <  1:161893/251‑1 (MQ=255)
                                                                                          cGGCGTTACTGGAGCAGACGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGCGCGGTCGCATAACGCTAAGCGCCGCAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGAGGTATTCCTGACTACGCGCTGTCCAGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAACTGGGCAAAATAGCAGCGGCCTGGGGTTGGCGTTCGTCAGTggtgtcgcacggctg                                                                                                                                                        >  2:644496/1‑237 (MQ=255)
                                                                                           ggCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCgtttgt                                                                                                                                                       <  1:7121/251‑1 (MQ=255)
                                                                                                    ggAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCAc                                                                                                                                                                                                                                                    <  1:490717/149‑1 (MQ=255)
                                                                                                    ggAGCAGGCACTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAAGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCAc                                                                                                                                                                                                                                                    >  2:490717/1‑149 (MQ=255)
                                                                                                       gCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGGCGCCCGTTTGTTTAACGGCGa                                                                                                                                             >  2:62712/1‑249 (MQ=255)
                                                                                                                    aaTTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCa                                                                                                                               >  1:524426/1‑250 (MQ=255)
                                                                                                                      tttACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAACTGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGACGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACg                                                                                                                            >  1:56344/1‑251 (MQ=255)
                                                                                                                       ttACTGGATATCGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGTACACGGCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGGCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGGCTGGGGTTGGCGTGCGTCAGTGAGGTCGCCCGTTTGTTTAACGTCGAAGTCACGCTGCGCAACg                                                                                                                            <  2:655592/250‑1 (MQ=255)
                                                                                                                        tACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGt                                                                                                                           <  1:184287/250‑1 (MQ=255)
                                                                                                                        tACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGCCGCATAACGCTAAGCGACGAAGTGGATCAGGAACACGTCGCCCTTATGGAGCTGGATACCGGTAGTGGTATTCCTGACTACGCGATGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGAGAGGCCGCCCGATTGTTTAACGGAGAAGTCACGCAGCGCaaaggg                                                                                                                          >  2:622637/1‑247 (MQ=255)
                                                                                                                            ggATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAgg                                                                                                                      <  1:459385/251‑1 (MQ=255)
                                                                                                                             gATAACGCCATCGATTTTACCCCCGAGAGCGGACGCATAACGCTAAGCGCCGAAGTGGATCAGGAACCCGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCCCGTATTTTTGAACGCTTTTACTCTTTGCCGCGGGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGGCAGTGAGGTCGCCCGTTTGTTTAACGGCGAATACACGCTGCGCTACGTGCAgg                                                                                                                      >  2:329964/1‑250 (MQ=255)
                                                                                                                               taactacATCGATTTTACACCCGAGAGCGGCCGCATAACGCTAAGCGCCGAAATTGATCAGGACCACGTCGACCTTAAGGTGCTGGTTACCGGAAGTGGTATTCCTGACTAAGCGCTGACACGTATTTTTGAACGAATTTCCTCTTTGCCGCGTGCACATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGCCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTAACGCTGCGCAACGTGCAGGa                                                                                                                     >  2:586035/7‑249 (MQ=255)
                                                                                                                                      aTCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaaaa                                                                                                                                                                                                    <  1:423594/163‑1 (MQ=255)
                                                                                                                                      aTCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaaaa                                                                                                                                                                                                    >  2:423594/1‑163 (MQ=255)
                                                                                                                                               aCCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGCCTACGCGCTGTCACGTATTTTTGAAAGCGTTTACTCTTTGTCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGAGAGGGCGCCCGTTTTTTTAACGGCGAAGGCACGCGGCGCGCAGTGCAGGAAGGGGGCGTGCTGGGCTc                                                                                                   >  2:558172/1‑251 (MQ=255)
                                                                                                                                                cccccGATAGCGGTCGCATAAAGCTAAGCGCCGAAGTGAATCAGGCACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCATGACTACGCGCTGTCACGTATTTATGAACGCTTTTACTCTTTTCCGCGTGCAAATTGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAAAGGAGGATCCACGCTGCGAAACGTACAGGAAGGGGGCGTgatgtccgcg                                                                                                  >  2:182666/1‑242 (MQ=255)
                                                                                                                                                cccccGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCATGCTGGCCTCg                                                                                                  >  1:352120/1‑251 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCCGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAg                                                                                                                                                                        >  2:465871/1‑171 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAg                                                                                                                                                                        <  1:465871/171‑1 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTc                                                                                        >  1:29651/1‑251 (MQ=255)
                                                                                                                                                          cGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGAGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTACGTGAGGTCCCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGAGCAGGCCTCGCTTCGTCttt                                                                                        >  2:527712/1‑250 (MQ=255)
                                                                                                                                                                    aCGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCg                                                                                             <  1:405711/236‑1 (MQ=255)
                                                                                                                                                                    aCGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGGGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCg                                                                                             >  2:405711/1‑236 (MQ=255)
                                                                                                                                                                        tAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCACATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCACCGTGCAGGAGGATGGCGTGCTGGCCTCGCTTCGACTTCACCGTCCCTTccc                                                                           >  2:220653/1‑248 (MQ=255)
                                                                                                                                                                               cGAAGTGGATCAGGAACACGTCACGCGTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCACGAAGGTGGCGTGCTTGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTc                                                                   <  1:69898/251‑1 (MQ=255)
                                                                                                                                                                                          aGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTccc                                                        <  1:14315/251‑1 (MQ=255)
                                                                                                                                                                                           gtgacacGTCGCCCTTAAGGTGCTGTATACCGGTAGTGGTATTCCTGATTACGCGCTGTCACGTATTTGTGAACGCTTTTTCGCTTTGCCGCGTGGAAATGGGCAAAAAAGCATCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGGGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCca                                                       <  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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC  >  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: