Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,896 A→G T401A (ACG→GCC)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,8960AG97.0% 187.8 / ‑4.2 69T401A (ACG→GCC) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (40/25);  total (42/27)
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.43e-01

AAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCA  >  NC_000913/4637660‑4638136
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
aaGTCGTTCTGACTGTTGTTGATGTGTCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTg                                                                                                                                                                                                                                    <  1:306417/251‑1 (MQ=255)
  gTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTgg                                                                                                                                                                                                                                   >  1:61088/1‑250 (MQ=255)
     gTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCCCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGCGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGCGAGCGGTCGCATAACGCTAAGCGCCGCAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATAc                                                                                                                                                                                                                               >  1:69701/1‑251 (MQ=255)
        cTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCTCCGTGCAGTTGGCAGAAAAAAACATCACTCTACCTGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCCAGAGCGGTCGCCTAACGCTAAGCGCCGAAATGGAGCAGGAACACGGCGCCCTTAAGGTGCTGGATAcccg                                                                                                                                                                                                                            >  1:430984/1‑249 (MQ=255)
          gactgtggTTGATGTGGCGGCATTATTTCTCCGCGTCAGCGAAGCGCGCACCGTGCCGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTACTGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGAAACCCGTCGCCCTTAAGGTGCTGGATACCGGt                                                                                                                                                                                                                           >  1:465482/8‑250 (MQ=255)
               ttgTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGGGCTGGATACCGGTAGTggg                                                                                                                                                                                                                     >  1:137586/2‑250 (MQ=255)
                      ggtgGCTGCATTATTTCGCCGCGTCAGCGAAGCGCGAACCGTGCAGTTGGAAGAAAAAAACATCACTCTACATGTTATTCCCACTGAGTTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAGCGCCATCGATTTTCCCCCCGAGAGCGGTCGCATATCGCTAAGCGCCGAAGTGGATCATGAACACGTCGCCCTTAAGGGGCTGGATACCGGTAGTCGTATTCCt                                                                                                                                                                                                               <  2:604728/249‑1 (MQ=255)
                           cggcATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGACAGAAAAAAACATAACTCTACATGTTATGCCCACTGCGGATAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTAc                                                                                                                                                                                                         >  1:309899/1‑251 (MQ=255)
                            ggcATTATTTCGCCGCGTCAGCGGAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTAcg                                                                                                                                                                                                        <  1:489582/251‑1 (MQ=255)
                                       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)
                                                                                                      ccactgagtgtAACGTTGCTGCTGAACCGGCGTGACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATGTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGCCTCCCTTAAGGTGCTGGATACCGGGAGTGGTATTCCTGACTACGCGCtgtc                                                                                                                                                                                                 <  2:95481/174‑4 (MQ=255)
                                                                                                      ccACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCtgtc                                                                                                                                                                                                 >  1:95481/1‑181 (MQ=255)
                                                                                                          tgAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGCTAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGCAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGCTACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTATTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGTCGt                                                                                                                          >  2:409127/2‑251 (MQ=255)
                                                                                                            aGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc                                                                                                                        <  2:281782/251‑1 (MQ=255)
                                                                                                            aGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc                                                                                                                        <  1:3515/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)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
AAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCA  >  NC_000913/4637660‑4638136

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: