Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,385 A→T T127S (ACG→TCG)  creB → DNA‑binding transcriptional regulator CreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,3850AT97.1% 200.1 / ‑4.2 68T127S (ACG→TCG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (38/28);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑ACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCT  >  NC_000913/4636140‑4636611
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
agcaggctcccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCg                                                                                                                                                                                                                               >  1:111256/9‑251 (MQ=255)
          cGACGCCATGATTCGCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTGACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCa                                                                                                                                                                                                                     <  2:228956/251‑1 (MQ=255)
           gaggttaTGATTCTCGATGTTGGTTTGCCTGATTTCAGCGGCTTTGAATTATGCCGCCAGTTACAGGCGCTCCATCCGGCGTTACCGGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACGATGTGGCTAAACCGTTTGCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCAt                                                                                                                                                                                                                    <  2:559795/245‑1 (MQ=255)
             cGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATcc                                                                                                                                                                                                                  <  1:541801/251‑1 (MQ=255)
              ggctttaTTCTCGATGGTGGTTTGCTTTATATCAGCGTCTTTGAAGTATGCAGCGAGTTCCTGGCGCTCCATCCGGCGTGACCTGTAATGTGCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGTGCTGGAAATTGGTGCTGAGGACTATGTGGCTAAACCGTTTTCTCCCTGCGAAGTGTGCGCCATGGTGCGCACCTTACTTCGTCGGGTGAAGAAGTTCTCT‑TCGCCGTCTCCCGTCATCCg                                                                                                                                                                                                                 <  2:380993/245‑1 (MQ=255)
               tCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCCTCCGt                                                                                                                                                                                                                >  2:140312/1‑251 (MQ=255)
                   gATTCTCGATGTTGGTTTGCTGGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCGGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTAtt                                                                                                                                                                                                             <  2:526141/250‑1 (MQ=255)
                       ctcGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGGGGTCGATCGCCTGCTTGGGCTGGCAGTTGGTGCTGACGACTATGTGGCTAATCCGTTTTCACCCCGCGACGTGTGCGCCAGGGTTCGCCCCTTACTGCGTCGGGGGAAGCAGTTCTCG‑TCGCAGTCTCCCATCATACGTGTTggtcc                                                                                                                                                                                                        >  2:401724/1‑248 (MQ=255)
                        tcGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCAt                                                                                                                                                                                                       >  1:535767/1‑251 (MQ=255)
                           aTGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCAtttt                                                                                                                                                                                                    <  1:98412/251‑1 (MQ=255)
                                 gtttGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaattg                                                                                                                                                                                              >  1:490375/4‑251 (MQ=255)
                                     gCCGGATTTCAGCGGATGTGAATTATGCCTCCAGTTACTGGCGCTCCATCCTGCTTTACCTGTACTGTTCCTGACGGCCCTAAGTGAAGAGGTCGATCGCCTTCTTGGGCTGGAAATTGGTGCGGACGACTATGTGGCTAATCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCCTCTTACTTCGTCGGGTGAATAAGTTCTCG‑TCGCCGTCTCCCGTCATCCTTATTGGTCATTTTTAATttaatg                                                                                                                                                                                          <  2:432344/251‑1 (MQ=255)
                                      cctgATATCAGCGTCTTGGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGGGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATtgaatga                                                                                                                                                                                         <  2:506981/248‑1 (MQ=255)
                                        tgATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGCCGAGGTCGATCGCCTGCTTGGGCTGGAACTTGGTGCTGCCGCCTATGTGGCTCAACCGTTTTCACCCCGCGAAGAGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGGGAAGAAGTTCTCGTAC‑CCGGCTCCCGTCATCCGTTTTGGTCCTTTTGACTtgaatgacc                                                                                                                                                                                       >  2:380739/2‑249 (MQ=255)
                                        tgATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGAACTGTTCCTGACGGCCCGAAGTGCAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGCGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTCCGCACCTTACTGCGTCGGGTGAAGACGATCTCG‑TCGCCGACTCCCGTCATCCGTATTGGTCATTTTGAATTGCATGAAc                                                                                                                                                                                       >  2:36553/2‑251 (MQ=255)
                                           tatcaGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCcg                                                                                                                                                                                    >  1:81994/5‑251 (MQ=255)
                                            accaGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCcgc                                                                                                                                                                                   <  2:318059/248‑1 (MQ=255)
                                               aGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCgcggc                                                                                                                                                                                >  1:36595/1‑251 (MQ=255)
                                                gCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCgcggcg                                                                                                                                                                               >  1:563450/1‑251 (MQ=255)
                                                   gCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAg                                                                                                                                                                            <  1:551445/251‑1 (MQ=255)
                                                           ttatgCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTAATCCGTATTGGTCATTTTGAATTGAATGAACCCGAGGCGCCGATCAGCTg                                                                                                                                                                    >  2:445279/4‑251 (MQ=255)
                                                            tatgCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTgg                                                                                                                                                                   >  1:139903/3‑251 (MQ=255)
                                                             atgCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCTAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGt                                                                                                                                                                  <  1:277938/250‑1 (MQ=255)
                                                             atgCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCgcggc                                                                                                                                                                                <  1:347231/236‑1 (MQ=255)
                                                             atgCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCgcggc                                                                                                                                                                                >  2:347231/2‑237 (MQ=255)
                                                               gccgccAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGGTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGtgt                                                                                                                                                                <  2:36595/251‑3 (MQ=255)
                                                                  gccAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGac                                                                                                                                                             <  1:60199/251‑1 (MQ=255)
                                                                    cAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGGGCGCCAGGGTGCGCACCTTACTGCGTCGGGGGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACTCGCGGCGCAGATCAGCGGGTTTGaaac                                                                                                                                                           >  2:468226/1‑248 (MQ=255)
                                                                     aGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGGAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACg                                                                                                                                                          >  1:615415/1‑251 (MQ=255)
                                                                     aGTTACTGGCGCCCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACg                                                                                                                                                          >  1:432833/1‑251 (MQ=255)
                                                                         aCTGGGGCTCGAGCCGTCGTTACCTGTGCTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTGGGGCTGGAAATTGGTGCTGACGACTAGGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGTAGATCAGCTGGTTTGACACGCCAt                                                                                                                                                      <  2:490375/251‑1 (MQ=255)
                                                                         aCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCGCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCAt                                                                                                                                                      <  1:91497/251‑1 (MQ=255)
                                                                               gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAacgc                                                                                                                                                <  1:445279/251‑1 (MQ=255)
                                                                                    aTCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGTAATTGGAGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTAGGGAGAAGCAGTTCTCG‑TCGCCGTCTCCCGTAATCCGTATTGGTCATTTTGAATTGAATGAACCCGCAGCGCATATCAGCTGGTTTGACACGCCagtaaagctgccg                                                                                                                                           >  2:665260/1‑239 (MQ=255)
                                                                                      ccggggttccctgtacttgcccTGACGGCCCGAAGTGAAGAGGTCGATCGCCGGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGACCCCGCGGCGCAGATCAGCGGGTTTGACACGCCATTAACGCTGAcgcg                                                                                                                                         <  1:149506/231‑4 (MQ=255)
                                                                                          cGTTACCTGTACTGTTCCTGACGGCCAGAAGTGAAGAGGTCGATCGCCTGATTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAACCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGCAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGACCAGCTGGATTGACACGCCTTTAACGCTGACGCGGTAt                                                                                                                                     >  2:71973/1‑251 (MQ=255)
                                                                                           gTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTtgaattgaat                                                                                                                                                                                           >  1:208515/1‑196 (MQ=255)
                                                                                           gTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTtgaattgaat                                                                                                                                                                                           <  2:208515/196‑1 (MQ=255)
                                                                                           gTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATg                                                                                                                                    >  1:482054/1‑251 (MQ=255)
                                                                                            ttACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATgg                                                                                                                                   >  1:321587/1‑250 (MQ=255)
                                                                                               ccTGTGCTGTTCCTGACGGCCCGGAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGGGCTGTCGACTAGGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCTTCTCCCGTCATCCGTTTTGGTCATTTTGAATTGAAGGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGtt                                                                                                                                <  2:583153/251‑1 (MQ=255)
                                                                                                      tGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAttg                                                                                                                         >  1:397432/1‑251 (MQ=255)
                                                                                                            tGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGCAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCAATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGa                                                                                                                     >  2:311061/1‑249 (MQ=255)
                                                                                                                     gAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACCATGTGGCTAAACCGTTTTACCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTc                                                                                                          <  1:334850/251‑1 (MQ=255)
                                                                                                                       ggTGAAGAGGTCGATCGCCTGCTTGGGCGGGAAATTGGTGCTGACGACGATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCTCTTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCaa                                                                                                        <  2:36546/250‑1 (MQ=255)
                                                                                                                                  cGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGCGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGcc                                                                                             >  1:327344/1‑251 (MQ=255)
                                                                                                                                   gATCGCCTGCTTGGGCTGGAAATTGGTGCTGAAGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAAACGTTCTCG‑ACGCCATCTCCCGTCATCCGTATTGGACATTTTGCATTGACTGAACCCGCGGCGCAGATCACCTGGTTTGCCCCACCATTAGCGATGACGCGGGATGAGTTTTTATTGTTGAAGTCGTTACTCAAGTCCCCGGGCcg                                                                                            >  2:482146/1‑251 (MQ=255)
                                                                                                                                      cGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTTGGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                          <  1:344495/250‑1 (MQ=255)
                                                                                                                                      cGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCGGAAGACGTTTCTCCAGTCGCCAGGCCGCGt                                                                                         >  1:596986/1‑251 (MQ=255)
                                                                                                                                         cTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCCTTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTg                                                                                      >  2:624645/1‑251 (MQ=255)
                                                                                                                                           gCTTGGGCTGGAAATGGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAATTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTgg                                                                                     <  1:457204/250‑1 (MQ=255)
                                                                                                                                               gggCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                 <  1:311061/250‑1 (MQ=255)
                                                                                                                                                        aTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCCACt                                                                       >  2:100626/1‑251 (MQ=255)
                                                                                                                                                            gTGCTGACGACTATGTGGCTAAACCGTTTTAACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTACGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATg                                                                   >  1:214944/1‑251 (MQ=255)
                                                                                                                                                                        atgTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTgg                                                       >  2:466021/3‑251 (MQ=255)
                                                                                                                                                                              cTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGa                                                     >  1:379891/1‑247 (MQ=255)
                                                                                                                                                                                aaaCCGTTTTCACCTCGCGAAGTGAGCGACAGGGTGCGCACCTTACTGCGTCGGGTAAAGAAGTTCTCG‑TCGCCGTCTCACGTCATACGTATGGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGGCAGGCACTTAACGCTGACGCGGAATGAGTTTTTATTTATTAAGACGTTACTCATGTCGCCAGTCCGCGACCGGTCCCGCCAGCCACTACTGGCAAGCGCGTGGGAAGATgc                                               >  2:60120/1‑251 (MQ=255)
                                                                                                                                                                                  aCCGTTTTCACCCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGCGAAACAGTTCTCG‑TCGCCGTCCCCCGTCCTCCGTATTGGTCCTTTTGTATTGAATGAACCCGCAGCGCAGATCCGCTGGTTTGACACGCAATTAACGCTGACGCGGTATGAGTTTTTCTTGCTGCAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGACGTTgcgc                                             >  1:443676/1‑251 (MQ=255)
                                                                                                                                                                                      ttttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAgg                                          >  1:542961/1‑250 (MQ=255)
                                                                                                                                                                                      ttttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAgg                                          >  1:159866/1‑250 (MQ=255)
                                                                                                                                                                                                 cgaagttagcgcCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTa                                                                                                                                      >  2:503482/9‑147 (MQ=255)
                                                                                                                                                                                                 cgAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCTTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCGGGTTTGACACGCCATTAACGCTGACGCGGTa                                                                                                                                      <  1:503482/147‑1 (MQ=255)
                                                                                                                                                                                                      tgtgCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                         >  1:557482/1‑251 (MQ=255)
                                                                                                                                                                                                      ggtgCGCCAGGGTGCTCCCCTTACTGCGTGGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCACCCGTATTGGTCATTTTGAATTGAATGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                         <  1:4903/250‑1 (MQ=255)
                                                                                                                                                                                                       gtgCGCCAGGGGGCGCTCCTTCCTGCGTCGGGTAAAGAATTTCTCG‑TCGCCGTCTCCCTTCACCCGTATTGGTCATTTTAAATGGAATGACCCCGCGCCGCAGATCAGCTGTTTTGCAACGCCATTACCGCGGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGTCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCg                        <  1:70950/251‑1 (MQ=255)
                                                                                                                                                                                                               gggTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGGCGCGGTATGAGTTTTTATTGCTGAAGCCGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCCACTAATGGATAGCGTGTGGAAAGAGGCGGAGGGCACCTCCGATCGCACCGttgctacc                >  2:21529/1‑244 (MQ=255)
                                                                                                                                                                                                                       aCCTTACTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTa         >  1:336388/1‑250 (MQ=255)
                                                                                                                                                                                                                             cTGCGTCGGGTGAAGAAGTTCTCG‑TCGCCGTCTACCGTCATCCGGATTGGTCATTTTGGATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTACTGGATAGCGTGTGGGAAGGTGCGCAGGAAACCTACGATCGCACCGTTGATCCCCACCTTAAAACgtt  >  2:221253/1‑249 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
AGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCG‑ACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCT  >  NC_000913/4636140‑4636611

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