Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,857 T→C C388R (TGC→CGC)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,8570TC98.6% 207.3 / ‑6.4 73C388R (TGC→CGC) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (44/28);  minor base G (1/0);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑ACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTC  >  NC_000913/4637616‑4638103
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
ggTAGAAACGTTACGACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACg                                                                                                                                                                                                                                                 <  1:77041/250‑1 (MQ=255)
   aGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGGCGTTCTGACTGATGTTGATGTGGCGGAGTTATTTCGCCGCGTCAGCGAAACGCGCACCGTGCAGTTGGCAGAAAAAACCATCACTCTACATGTTATGCCCCCTGAGGTTACTGTGGCTGCTGAACCGGCGTTACTGGAGGAGGCGCTTGGGAATTTACTTGATAACGACATCGATTTTACCCCCGAGATCGGTCGCATAACGcgat                                                                                                                                                                                                                                             >  2:314724/1‑248 (MQ=255)
           tactacGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGa                                                                                                                                                                                                                                      <  1:178703/250‑1 (MQ=255)
           tactacGCCAGGCAAGACTGGAGAATCGTCAGGAAGCCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCCCCGGGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCCTCGATTTTACCCCCGCGAGCGGTCGCATAACCCTAAGCGCCGa                                                                                                                                                                                                                                      >  1:587465/1‑250 (MQ=255)
            actacGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATACCGCTAAGCGCCGa                                                                                                                                                                                                                                      >  1:123975/1‑249 (MQ=255)
                 gCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGAt                                                                                                                                                                                                                               <  1:352773/251‑1 (MQ=255)
                       cAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTGTTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCCCTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGCACCGGCGGTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCTATTTTACCCCCGAGGGGGGTCGCATTATGCTAAGCGCCGAAGTggggcaggg                                                                                                                                                                                                                          >  2:459385/1‑243 (MQ=255)
                            cTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGCTCAGGAACACGt                                                                                                                                                                                                                    >  1:186844/1‑251 (MQ=255)
                              ggAGAATCGTCATGAAGTCGTTCTGACTGTTGTTGATGTGGCTGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTcg                                                                                                                                                                                                                  <  2:526041/251‑2 (MQ=255)
                              ggAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGt                                                                                                                                                                                                                    <  1:306957/249‑1 (MQ=255)
                                            aaGTCGTTCTGACTGTTGTTGATGTGTCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTg                                                                                                                                                                                                   <  1:306417/251‑1 (MQ=255)
                                              gTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTgg                                                                                                                                                                                                  >  1:61088/1‑250 (MQ=255)
                                                 gTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCCCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGCGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGCGAGCGGTCGCATAACGCTAAGCGCCGCAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATAc                                                                                                                                                                                              >  1:69701/1‑251 (MQ=255)
                                                    cTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCTCCGTGCAGTTGGCAGAAAAAAACATCACTCTACCTGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCCAGAGCGGTCGCCTAACGCTAAGCGCCGAAATGGAGCAGGAACACGGC‑GCCCTTAAGGTGCTGGATAcccg                                                                                                                                                                                           >  1:430984/1‑249 (MQ=255)
                                                      gactgtggTTGATGTGGCGGCATTATTTCTCCGCGTCAGCGAAGCGCGCACCGTGCCGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTACTGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGAAACCCGTC‑GCCCTTAAGGTGCTGGATACCGGt                                                                                                                                                                                          >  1:465482/8‑250 (MQ=255)
                                                           ttgTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGGGCTGGATACCGGTAGTggg                                                                                                                                                                                    >  1:137586/2‑250 (MQ=255)
                                                                  ggtgGCTGCATTATTTCGCCGCGTCAGCGAAGCGCGAACCGTGCAGTTGGAAGAAAAAAACATCACTCTACATGTTATTCCCACTGAGTTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAGCGCCATCGATTTTCCCCCCGAGAGCGGTCGCATATCGCTAAGCGCCGAAGTGGATCATGAACACGTC‑GCCCTTAAGGGGCTGGATACCGGTAGTCGTATTCCt                                                                                                                                                                              <  2:604728/249‑1 (MQ=255)
                                                                       cggcATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGACAGAAAAAAACATAACTCTACATGTTATGCCCACTGCGGATAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTAc                                                                                                                                                                        >  1:309899/1‑251 (MQ=255)
                                                                        ggcATTATTTCGCCGCGTCAGCGGAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTAcg                                                                                                                                                                       <  1:489582/251‑1 (MQ=255)
                                                                                   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)
                                                                                                                                                  ccactgagtgtAACGTTGCTGCTGAACCGGCGTGACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATGTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGCC‑TCCCTTAAGGTGCTGGATACCGGGAGTGGTATTCCTGACTACGCGCtgtc                                                                                                                                                                <  2:95481/174‑4 (MQ=255)
                                                                                                                                                  ccACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑GCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCtgtc                                                                                                                                                                >  1:95481/1‑181 (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)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
GGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTC‑ACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTC  >  NC_000913/4637616‑4638103

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