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

CAATTGAAGCGGTAAACGGAATGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTC  >  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)
                                          aGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACt                                                                                                                                                                                                      <  1:370507/251‑1 (MQ=255)
                                             aaCACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTcc                                                                                                                                                                                                   >  1:101561/1‑251 (MQ=255)
                                                 cGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAgg                                                                                                                                                                                               >  1:340250/1‑251 (MQ=255)
                                                       ggCGTCTTTATTGCGCGCTGTCATCCCGCTTCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGca                                                                                                                                                                                         <  2:269275/251‑1 (MQ=255)
                                                       ggCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTGTACTTCCCAGGATTGca                                                                                                                                                                                         <  1:379951/251‑1 (MQ=255)
                                                         cGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCAAGGATTGcacc                                                                                                                                                                                       >  2:495647/1‑250 (MQ=255)
                                                                attaCGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGttt                                                                                                                                                                                 >  2:262913/1‑250 (MQ=255)
                                                                       gccggcATCACGCTGCCGCCGCCGTTCAGCCAGGTACAACTTATTGCCCTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGTAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGACGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCCCTGCAGGCGTTTTCATTCCCCAGGTGTCACACTGTTTAActggcc                                                                                                                                                                         >  2:523292/6‑250 (MQ=255)
                                                                       cccTTCATCCCCCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAActggct                                                                                                                                                                         <  2:73130/250‑1 (MQ=255)
                                                                            cATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGgcgc                                                                                                                                                                    >  1:279798/1‑251 (MQ=255)
                                                                                    tGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATc                                                                                                                                                            >  1:199916/1‑251 (MQ=255)
                                                                                     gccgccgccGTTCAGCCCGGTACAACTTATTGAACTGCAGGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCa                                                                                                                                                           <  1:221655/251‑1 (MQ=255)
                                                                                       cgccgccgTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCAcc                                                                                                                                                         >  1:144436/1‑251 (MQ=255)
                                                                                             cgTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCCACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTCACTGGCTGGCGCGGATAATCACCCCGGtt                                                                                                                                                   >  1:361561/1‑251 (MQ=255)
                                                                                                    cccGGTACAACTTATTGAACTGCATGTTCTCAAATGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACATTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAAt                                                                                                                                            <  1:238000/251‑1 (MQ=255)
                                                                                                          aCAACTTATTGAACTGCATGGTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTAGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGc                                                                                                                                      <  1:506376/251‑1 (MQ=255)
                                                                                                                    gAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAAtgatga                                                                                                                            <  1:103233/251‑1 (MQ=255)
                                                                                                                         gCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATgaagaa                                                                                                                        <  1:126915/250‑1 (MQ=255)
                                                                                                                         gCATGTTCTCAAAAGCAACTTCTATTACCGCTACAATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAAc                                                                                                                       <  1:149400/251‑1 (MQ=255)
                                                                                                                           aTGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGCTAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAAc                                                                                                                       >  1:191074/1‑249 (MQ=255)
                                                                                                                                  cAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTgg                                                                                                              <  1:336304/251‑1 (MQ=255)
                                                                                                                                             tCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGATACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGttt                                                                                                    <  1:85072/250‑1 (MQ=255)
                                                                                                                                               tATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAAAAGCAGGTGGCGCTGCGTTTGCt                                                                                                 >  1:75749/1‑251 (MQ=255)
                                                                                                                                                aTTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGt                                                                                                                                                                                                             <  1:190016/143‑1 (MQ=255)
                                                                                                                                                aTTACCGCTACCATGATGATGGCAGCGATGGGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCt                                                                                                 <  2:141338/250‑1 (MQ=255)
                                                                                                                                                aTTACCGCTACCAGGAGTATGGCGGTGATGTGGCGGCAACAACACAGTTTCAGGGCGAGATGGTCGATTATTGGCTGCACTCCGTCCGTCTCGGCAGTAGTTGAATGGCGGCGCATCGCTTTAGTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCGCCCCGGTTCTGCAATCATGGCTCAAGGATGAAGAACAGCGGGTGGCGCGGCGTctgctg                                                                                                <  2:104687/251‑1 (MQ=255)
                                                                                                                                                 ttACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTgg                                                                                               >  1:81737/1‑251 (MQ=255)
                                                                                                                                                   aCCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGgag                                                                                             <  2:199916/251‑1 (MQ=255)
                                                                                                                                                      gCTACCATGATGATGGCAGCGATGGGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGaaa                                                                                          <  2:257204/251‑1 (MQ=255)
                                                                                                                                                      gCTACCATGATGAAGGCAGAGAAGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATGCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTGAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGaaa                                                                                          <  2:509042/251‑1 (MQ=255)
                                                                                                                                                        tACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGGGaa                                                                                           >  2:224788/1‑248 (MQ=255)
                                                                                                                                                             ggttgatgGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGAGTATTCGCGTCACGCCGGCCGTCTCAGCAGTAGGTGAATGGCGGAGCTACGCTTTATGCGCACCCACGGCAGTCGGTTTACTTCCCAGGATTGCACTCTGTTTAACTGGCTGGCGCGGATAATCACGCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCCTCGTTTGCTGGAGAAAGATcgcg                                                                                   <  2:309071/248‑1 (MQ=255)
                                                                                                                                                                    ggCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGGTTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGatcatc                                                                             <  1:328578/250‑1 (MQ=255)
                                                                                                                                                                     gCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTAGATTATTCGCGTCACGCCGTCCTGCTCGGCAGTAGTTGAATGGCGGAGCTAAGCTTTATTCGCACCCACGGCAGTCGTGGTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGATCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGatcatcag                                                                           <  2:361561/251‑2 (MQ=255)
                                                                                                                                                                         cGATGTGACGGCAACAACAGAGTATCAGGGCTAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGCATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCggg                                                                       >  1:254718/1‑251 (MQ=255)
                                                                                                                                                                          gATGTGACGACAACAACAAAGTATCAGGGCGAGCTGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTCGTTGAATGGCGGAGCTACGCTTTATTCGCCCCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATCATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGTACAGCAGGTGGCGCTGCGCTTGCTGGAGAAAGATCGCGCTCATCATCgggg                                                                      >  2:232040/1‑250 (MQ=255)
                                                                                                                                                                               gACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCTGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTg                                                                  <  2:326945/250‑1 (MQ=255)
                                                                                                                                                                                    caacaaCAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGat                                                            <  1:365169/251‑1 (MQ=255)
                                                                                                                                                                                     aacaacAGAGTATCAGGGCAAGATGGTCGATTATTCGCGTCACGCCGTCCTACTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTACTCGCACCCACGGCAGTCGTTTTACTTCCCATGATTGCACACTGTTTAACTGGCTGGCGCGGGTAATCACCCCGGTTCTTCAATCATCGCTCCATGCTTAAGCACAGCATGTGGCGCTTCGTTTGCTGGAGAAAGATCGCGATCATAATCGTGTAATGgatgat                                                            >  2:73146/1‑250 (MQ=255)
                                                                                                                                                                                       caacaGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATAtt                                                         >  1:134141/1‑251 (MQ=255)
                                                                                                                                                                                       caacaGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCCTCTCGGCAGGAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGCTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGtggacagt                                                         <  2:21112/251‑8 (MQ=255)
                                                                                                                                                                                                   aGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTg                                             >  1:80769/1‑251 (MQ=255)
                                                                                                                                                                                                     ggCGCGATGGTCGATTATTCGCGGCACGCCGTCCTTCTCGTCAGTAGTTGAATGGCGGAGCTACGCTGTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTTGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTCCTAATGCAGTGCt                                           <  2:434443/251‑1 (MQ=255)
                                                                                                                                                                                                     ggCGAGATGTGCGATTATTCGCGTCCCGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCATGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCt                                           <  2:434441/251‑1 (MQ=255)
                                                                                                                                                                                                         cgaTGTTCGATTATTCGCGTCACGCCGTCCTTCTAGGCCGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCACGGATTGCACACTGTTTCACTAGCTGGCGCGGATAATCACCCCGGTTCTGCAATATTGGCTCAATGATGAAGAAAAGCAGGTGGCGCTGCCTGTGCTGGAGCAAGTTCGCGATCATACTCGGGGAATGGTTTTTATTAATAATGGAGTGCTGTc                                        >  2:354292/2‑250 (MQ=255)
                                                                                                                                                                                                             ggTCGATTATTCTCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATc                                   <  2:340250/251‑1 (MQ=255)
                                                                                                                                                                                                             ggTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATc                                   <  1:48132/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCTGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGGGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTg                  >  1:122088/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ccTTCTCGGCAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTg           >  1:277372/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ctcGGCAGTAGTTGAATGGCGGAGCTAAGCTTTATTCGCCCCCACGGAAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGAGCGGATAATCACACCGGTTCTGCAATCATGGCTCAATGCTGAAGAACAGCAGGTGGCGCTGCGTTTGCTGAATAAAGATCACGATCATCATCGGGTACTGGTTGCTATTACTAATGCAGTGCTGTCACTTCTTGATCTCGACGATCTGATCGCTGACGt       >  1:284250/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGGCAGTAGTTGAATGGCGGAGCTACGCTTTAGTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCGTGATCTCGACGATCTGATCGCTGACGTc      <  2:277372/250‑1 (MQ=255)
                                                                                                                                                                                                                                              cAGTAGTTGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATATCGACGGTCTGATCGCTGGCGTCGCTc  >  2:16206/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CAATTGAAGCGGTAAACGGAATGCTCGCGGAGCGGTTAGCACAGAAACACGGCAAGGCGTCTTTATTACGCGCCTTCATCCCGCTGCCGCCGCCGTTCAGCCCGGTACAACTTATTGAACTGCATGTTCTCAAAAGCAACTTCTATTACCGCTACCATGATGATGGCAGCGATGTGACGGCAACAACAGAGTATCAGGGCGAGATGGTCGATTATTCGCGTCACGCCGTCCTTCTCGGCAGTAGTGGAATGGCGGAGCTACGCTTTATTCGCACCCACGGCAGTCGTTTTACTTCCCAGGATTGCACACTGTTTAACTGGCTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTC  >  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: