Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,612,191 G→T G98* (GGA→TGA)  hyfR → DNA‑binding transcriptional activator HyfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,612,1910GT97.1% 193.9 / ‑1.2 69G98* (GGA→TGA) hyfRDNA‑binding transcriptional activator HyfR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (31/36);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

CAATTGAAGCGGTAAACGGAATGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTC  >  NC_000913/2611946‑2612434
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cAATTGAAGCGGTAAACGGAATGCTCGCGGAGCGATTAGCACAGAAACACGGCAAGGCGACTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAACAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGAGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCCGTCGTTGAATg                                                                                                                                                                                                                                                 >  2:333087/1‑251 (MQ=255)
    tGAAGCGGTAAACGGAATGCTCGCGGAGCGGTTAGCACAGAAACACGGCCAGGCGTCTTTATTACGCGACTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGGTGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATggagg                                                                                                                                                                                                                                             >  1:21112/1‑248 (MQ=255)
             aaaCGGAATGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCt                                                                                                                                                                                                                                    <  1:36091/251‑1 (MQ=255)
                cGGAATGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCttt                                                                                                                                                                                                                                  >  2:48554/1‑250 (MQ=255)
                   caTGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGCTGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGACTGGCGGCGCTCCGCTTTATTc                                                                                                                                                                                                                              >  2:48132/2‑251 (MQ=255)
                   aaTGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTc                                                                                                                                                                                                                              <  1:201074/251‑1 (MQ=255)
                      gcttgcgGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCACCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGc                                                                                                                                                                                                                            <  2:83152/246‑1 (MQ=255)
                       cTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGGAGTTGAATGGCGGAGCTACGCTTTATTCGCAc                                                                                                                                                                                                                          >  2:149263/1‑251 (MQ=255)
                          gcgGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCAccc                                                                                                                                                                                                                        >  1:434443/1‑250 (MQ=255)
                          gcgGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCAccc                                                                                                                                                                                                                        >  1:434441/1‑250 (MQ=255)
                                gggTTAGCACAGAAACACGGCAAGTCGTCTTTATTAGCCGCCTTCACCCCGCTGCCGCCGCTGTTCATCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCGACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCTGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCTCACCCACGGCAg                                                                                                                                                                                                                 <  2:174732/250‑1 (MQ=255)
                                 ggTTAGCACAGAAACTCGGCAAGGCGTCTTTATTACGCGCCTTCTACCCGCTGCCGCCGCCGTTCAGCCCTGTACAACTTATTCAACTGCATGTTCTCAAAAGCATCTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCGTCTCTGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGt                                                                                                                                                                                                                <  1:57539/251‑1 (MQ=255)
                                 ggTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGt                                                                                                                                                                                                                <  2:513605/251‑1 (MQ=255)
                                   ttAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCCCCCACGGCAGTCg                                                                                                                                                                                                              >  1:73130/1‑251 (MQ=255)
                                   ttAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCAACGGCAGTCg                                                                                                                                                                                                              >  2:114649/1‑251 (MQ=255)
                                     agcccaGAACCACGGCAAGGCTTCTTTATTACGCGCCTTCATCCCGCTGCCACCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGATTTATTCGCACCCACGGCAGTCGt                                                                                                                                                                                                             <  1:15743/246‑1 (MQ=255)
                                       cacaGAAACACGGCATGGTGTCTTTATTCCTCGCTTTCACCCCGATGCCGCCGGCGGTCTGCCCTGTACAACTTATGGAACTGCATGTTCTCAAAAGCAATTTCTATTTCCGATACCATGATGATGGCAGCGATGTGACGGCAACACCAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTTAATGGCGGTGCTACGCTTTATTCGCACCCACGGCAGTCGtttt                                                                                                                                                                                                          <  2:503627/251‑1 (MQ=255)
                                          aGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACt                                                                                                                                                                                                      <  1:370507/251‑1 (MQ=255)
                                             aaCACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTcc                                                                                                                                                                                                   >  1:101561/1‑251 (MQ=255)
                                                 cGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAgg                                                                                                                                                                                               >  1:340250/1‑251 (MQ=255)
                                                       ggCGTCTTTATTGCGCGCTGTCATCCCGCTTCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGca                                                                                                                                                                                         <  2:269275/251‑1 (MQ=255)
                                                       ggCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTGT‑ACTTCCCAGGATTGca                                                                                                                                                                                         <  1:379951/251‑1 (MQ=255)
                                                         cGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCAAGGATTGcacc                                                                                                                                                                                       >  2:495647/1‑250 (MQ=255)
                                                                attaCGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGttt                                                                                                                                                                                 >  2:262913/1‑250 (MQ=255)
                                                                       gccggcATCACGCTGCCGCCGCCGTTCAGCCAGGTACAACTTATTGCCCTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGTAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGACGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCCCTGCAGGCGTTTTCATTCCCCAGG‑TGTCACACTGTTTAActggcc                                                                                                                                                                         >  2:523292/6‑250 (MQ=255)
                                                                       cccTTCATCCCCCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAActggct                                                                                                                                                                         <  2:73130/250‑1 (MQ=255)
                                                                            cATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGgcgc                                                                                                                                                                    >  1:279798/1‑251 (MQ=255)
                                                                                    tGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATc                                                                                                                                                            >  1:199916/1‑251 (MQ=255)
                                                                                     gccgccgccGTTCAGCCCGGTACAACTTATTGAACTGCAGGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCa                                                                                                                                                           <  1:221655/251‑1 (MQ=255)
                                                                                       cgccgccgTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCAcc                                                                                                                                                         >  1:144436/1‑251 (MQ=255)
                                                                                             cgTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCCACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTCACTGGCTGGCGCGGATAATCACCCCGGtt                                                                                                                                                   >  1:361561/1‑251 (MQ=255)
                                                                                                    cccGGTACAACTTATTGAACTGCATGTTCTCAAATGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACATTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAAt                                                                                                                                            <  1:238000/251‑1 (MQ=255)
                                                                                                          aCAACTTATTGAACTGCATGGTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTAGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGc                                                                                                                                      <  1:506376/251‑1 (MQ=255)
                                                                                                                    gAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAAtgatga                                                                                                                            <  1:103233/251‑1 (MQ=255)
                                                                                                                         gCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATgaagaa                                                                                                                        <  1:126915/250‑1 (MQ=255)
                                                                                                                         gCATGTTCTCAAAAGCAACTTCTATTACCGCTACAATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAAc                                                                                                                       <  1:149400/251‑1 (MQ=255)
                                                                                                                           aTGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGCTAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAAc                                                                                                                       >  1:191074/1‑249 (MQ=255)
                                                                                                                                  cAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTgg                                                                                                              <  1:336304/251‑1 (MQ=255)
                                                                                                                                             tCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGATACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGttt                                                                                                    <  1:85072/250‑1 (MQ=255)
                                                                                                                                               tATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAAAAGCAGGTGGCGCTGCGTTTGCt                                                                                                 >  1:75749/1‑251 (MQ=255)
                                                                                                                                                aTTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGt                                                                                                                                                                                                             <  1:190016/143‑1 (MQ=255)
                                                                                                                                                aTTACCGCTACCATGATGATGGCAGCGATGGGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCt                                                                                                 <  2:141338/250‑1 (MQ=255)
                                                                                                                                                aTTACCGCTACCAGGAGTATGGCGGTGATGTGGCGGCAACAACACAGTTTCAGGGCGAGATGGTCGATTATTGGCTGCACTCCGTCCGTCTCGGCAGTAGTTGAATGGCGGCGCATCGCTTTAGTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCGCCCCGGTTCTGCAATCATGGCTCAAGGATGAAGAACAGCGGGTGGCGCGGCGTctgctg                                                                                                <  2:104687/251‑1 (MQ=255)
                                                                                                                                                 ttACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTgg                                                                                               >  1:81737/1‑251 (MQ=255)
                                                                                                                                                   aCCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGgag                                                                                             <  2:199916/251‑1 (MQ=255)
                                                                                                                                                      gCTACCATGATGATGGCAGCGATGGGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGaaa                                                                                          <  2:257204/251‑1 (MQ=255)
                                                                                                                                                      gCTACCATGATGAAGGCAGAGAAGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATGCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTGAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGaaa                                                                                          <  2:509042/251‑1 (MQ=255)
                                                                                                                                                        tACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGGGaa                                                                                           >  2:224788/1‑248 (MQ=255)
                                                                                                                                                             ggttgatgGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGAGTATTCGCGTCACGCCGGCCGTCTCAGCAGTAGGTGAATGGCGGAGCTACGCTTTATGCGCACCCACGGCAGTCGGTTT‑ACTTCCCAGGATTGCACTCTGTTTAACTGGCTGGCGCGGATAATCACGCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCCTCGTTTGCTGGAGAAAGATcgcg                                                                                   <  2:309071/248‑1 (MQ=255)
                                                                                                                                                                    ggCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGGTTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGatcatc                                                                             <  1:328578/250‑1 (MQ=255)
                                                                                                                                                                     gCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTAGATTATTCGCGTCACGCCGTCCTGCTCGGCAGTAGTTGAATGGCGGAGCTAAGCTTTATTCGCACCCACGGCAGTCGTGGT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGATCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGatcatcag                                                                           <  2:361561/251‑2 (MQ=255)
                                                                                                                                                                         cGATGTGACGGCAACAACAGAGTATCAGGGCTAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGCATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCggg                                                                       >  1:254718/1‑251 (MQ=255)
                                                                                                                                                                          gATGTGACGACAACAACAAAGTATCAGGGCGAGCTGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTCGTTGAATGGCGGAGCTACGCTTTATTCGCCCCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATCATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGTACAGCAGGTGGCGCTGCGCTTGCTGGAGAAAGATCGCGCTCATCATCgggg                                                                      >  2:232040/1‑250 (MQ=255)
                                                                                                                                                                               gACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCTGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTg                                                                  <  2:326945/250‑1 (MQ=255)
                                                                                                                                                                                    caacaaCAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGat                                                            <  1:365169/251‑1 (MQ=255)
                                                                                                                                                                                     aacaacAGAGTATCAGGGCAAGATGGTCGATTATTCGCGTCACGCCGTCCTACTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTACTCGCACCCACGGCAGTCGTTTT‑ACTTCCCATGATTGCACACTGTTTAACTGGCTGGCGCGGGTAATCACCCCGGTTCTTCAATCATCGCTCCATGCTTAAGCACAGCATGTGGCGCTTCGTTTGCTGGAGAAAGATCGCGATCATAATCGTGTAATGgatgat                                                            >  2:73146/1‑250 (MQ=255)
                                                                                                                                                                                       caacaGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATAtt                                                         >  1:134141/1‑251 (MQ=255)
                                                                                                                                                                                       caacaGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCCTCTCGGCAGGAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGCTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGtggacagt                                                         <  2:21112/251‑8 (MQ=255)
                                                                                                                                                                                                   aGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTg                                             >  1:80769/1‑251 (MQ=255)
                                                                                                                                                                                                     ggCGCGATGGTCGATTATTCGCGGCACGCCGTCCTTCTCGTCAGTAGTTGAATGGCGGAGCTACGCTGTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTTGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTCCTAATGCAGTGCt                                           <  2:434443/251‑1 (MQ=255)
                                                                                                                                                                                                     ggCGAGATGTGCGATTATTCGCGTCCCGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCATGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCt                                           <  2:434441/251‑1 (MQ=255)
                                                                                                                                                                                                         cgaTGTTCGATTATTCGCGTCACGCCGTCCTTCTAGGCCGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCACGGATTGCACACTGTTTCACTAGCTGGCGCGGATAATCACCCCGGTTCTGCAATATTGGCTCAATGATGAAGAAAAGCAGGTGGCGCTGCCTGTGCTGGAGCAAGTTCGCGATCATACTCGGGGAATGGTTTTTATTAATAATGGAGTGCTGTc                                        >  2:354292/2‑250 (MQ=255)
                                                                                                                                                                                                             ggTCGATTATTCTCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATc                                   <  2:340250/251‑1 (MQ=255)
                                                                                                                                                                                                             ggTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATc                                   <  1:48132/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCTGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGGGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTg                  >  1:122088/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ccTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTg           >  1:277372/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ctcGGCAGTAGTTGAATGGCGGAGCTAAGCTTTATTCGCCCCCACGGAAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGAGCGGATAATCACACCGGTTCTGCAATCATGGCTCAATGCTGAAGAACAGCAGGTGGCGCTGCGTTTGCTGAATAAAGATCACGATCATCATCGGGTACTGGTTGCTATTACTAATGCAGTGCTGTCACTTCTTGATCTCGACGATCTGATCGCTGACGt       >  1:284250/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGGCAGTAGTTGAATGGCGGAGCTACGCTTTAGTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCGTGATCTCGACGATCTGATCGCTGACGTc      <  2:277372/250‑1 (MQ=255)
                                                                                                                                                                                                                                              cAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATATCGACGGTCTGATCGCTGGCGTCGCTc  >  2:16206/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CAATTGAAGCGGTAAACGGAATGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTT‑ACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTC  >  NC_000913/2611946‑2612434

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