Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,411 A→T G135G (GGA→GGT 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,4110AT92.3% 138.0 / 1.6 54G135G (GGA→GGTcreBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (24/24);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACC  >  NC_000913/4636165‑4636650
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
cGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGAc                                                                                                                                                                                                                                                   <  1:16259/248‑1 (MQ=255)
 gATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCAttt                                                                                                                                                                                                                                               >  1:48596/1‑251 (MQ=255)
 gATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGTTCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACAttt                                                                                                                                                                                                                                               >  1:73296/1‑251 (MQ=255)
    gTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttg                                                                                                                                                                                                                                             <  1:12032/250‑1 (MQ=255)
          ttGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGGCATTttgaattga                                                                                                                                                                                                                                      >  2:11497/2‑250 (MQ=255)
               tgATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGACGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGACGGGTGAAGAAGTTCTCGTCGCCGTATCCCGTCCTCCGTATTGGTCATTTTGAATTGAATGAAc                                                                                                                                                                                                                                >  2:403328/2‑251 (MQ=255)
                             ttGAATTAGGCCGCCAGTTGCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGa                                                                                                                                                                                                                    <  2:447001/249‑1 (MQ=255)
                               gAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCGCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCa                                                                                                                                                                                                                 <  2:479954/250‑1 (MQ=255)
                                aaTTATGCCGCCTGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGc                                                                                                                                                                                                               <  2:158594/251‑1 (MQ=255)
                                 attatgcggccAGTTACTGGCGCTCCATCATGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGCAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGAGGCGCACCTTACTGCGTCGGGTGTCGAAGTTCTCGTCGCCGCCTCCCGTCCTCCGTATTGGTCATTTTTAATTGAATGAACCAGCGgcgccgtacagct                                                                                                                                                                                                              >  2:354579/9‑242 (MQ=255)
                                   tatgCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTAATGTTCCTGACGGTCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGTAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGAGCGCACCTTACTGAGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGGCATTTTGAGTTGAATGAACCCGCGGCGCAGTTCAGCTgg                                                                                                                                                                                                           >  2:90798/3‑251 (MQ=255)
                                           cAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGGGAAGAAGTTCTCGTCGCCGGCTCCCGTCATCCGGATTGGTCATTTTGAATTGAGTGAACCCGCGGCGCCGCTCAGCTGGTTTGacac                                                                                                                                                                                                   >  2:309595/1‑251 (MQ=255)
                                                aCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTGGACACGCCa                                                                                                                                                                                               <  2:284138/250‑1 (MQ=255)
                                                aCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCGTACGGCCCGAAGTGAAGAGGTCGATCGCCGGCTTGGGCTGGAAATTGGTGCTGACGATTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCAt                                                                                                                                                                                              <  2:264949/251‑1 (MQ=255)
                                                 cTGGCGATCCATCCGGCGTTTCCTGTACTGTTCCTGACGGCCCGAAGTGACGAGGTCGATCGCATGCTTGGGCTGTAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCAACCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCCTTTTGAATTGAATGAACCCGGGGCGACGATCAGCTTGTTTGACACGCCAtt                                                                                                                                                                                             >  2:311463/1‑251 (MQ=255)
                                                   ggCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTa                                                                                                                                                                                            >  2:329550/1‑250 (MQ=255)
                                                     cgcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATtgaatga                                                                                                                                                                                                                                  >  2:82883/1‑211 (MQ=255)
                                                     cgcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATtgaatga                                                                                                                                                                                                                                  <  1:82883/211‑1 (MQ=255)
                                                     cgcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTaacg                                                                                                                                                                                         <  1:336927/251‑4 (MQ=255)
                                                         ccttCCGGCTTTACCTGTACTTTTCCTGACGGCCCGATGTGAAGAGGTCGACCGCCTGCTTGGGCTGGAAATTGGGGCGGACGTCTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGTTGCGCCCCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTAATTTTGAATTGAATGACCCCGCGGCGAAGATCAGCTGGTTTGACACGCCATTAACGCTg                                                                                                                                                                                      <  1:265946/247‑1 (MQ=255)
                                                             ccgtcgttactggttcTGTACCTGACGCCCTGATGGGAAGAGGTCGAGCGCCTGCTTGGGCTGGAAATGGTGGCTGACGACTTTGTGGCTAACCCGTTTTCGCCCGGCGAAGTGTGCGCCGGGGTTCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTAGTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGACCAGCGGGTTTGACACGCCATTAACGCTGAcgcg                                                                                                                                                                                 <  2:350709/236‑4 (MQ=255)
                                                                     ccttGTACTGTTCCTGACGGCCCGAAGGGAGGAGGTCGATCGCCTGCTTGGGCTGGAAAGTGGTGCTCACGCATATGTGGCTAAACCTTTTTCCCCCCGCTAAGTTTGCCCCAGGTGGCGCACTTTACTGCGTCGGTTGAAGAATTTCTCGTCGCCGTCTCCCGTCATCCTTTTTGTTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGt                                                                                                                                                                        <  2:342039/248‑1 (MQ=255)
                                                                      cctgtatggTTCCTGTCGGCCCGAAGGGAAGAGGCCGAGCGCCTGCTTGGGTTGGAAATGGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGAGAAGTGTGCGCCAGGGTGCGCCCCTTACTGCGGCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCTTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGtt                                                                                                                                                                       <  2:429059/243‑1 (MQ=255)
                                                                          tACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCCTTAACGCTGACGCGGGATGAGttttt                                                                                                                                                                    >  2:494703/1‑250 (MQ=255)
                                                                              gtgccggACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAttg                                                                                                                                                                <  1:103301/244‑1 (MQ=255)
                                                                              gTTCCTGACGGCCAGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGAAGACTATGTGGCTCAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCCGACCAGATGGTTTGACACGCCCTTAACGCTGAACCGGTATGAGTTTTTAttgc                                                                                                                                                               >  2:259027/1‑250 (MQ=255)
                                                                                    aaCGGCCCGAAGGGACGGGGTCGAGTGCCTGCTTGGGTGGGAAATTGGTGCTGACGACTATGTGGCTAACCCGTTTCCACCCCGCGAAGTGGGCTCGAGGTTGCGCACCTTACGTCGTCGGGTGAAGAAGTTCTAGTAGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGTGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGa                                                                                                                                                         <  2:239618/250‑1 (MQ=255)
                                                                                         cccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTa                                                                                                                                                    <  1:159992/251‑1 (MQ=255)
                                                                                          ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGCAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTAc                                                                                                                                                   >  1:200455/1‑251 (MQ=255)
                                                                                                 gAAGAGGGCGATCGCTTGCTGGGGCAGTAAATTGGTGCTGACGACTATGTGGCTAAACCTTTTTCACCCCGCGAAGTGTGCTCCAGGGTGCGTATCTTACGGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTTCTGAAGACGGTACTCAAGTc                                                                                                                                            <  1:436268/251‑1 (MQ=255)
                                                                                                  aaGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcg                                                                                                                                           <  1:224148/251‑2 (MQ=255)
                                                                                                      ggTCGATCGCCTGCTTGGGCTGGAAATTGGTGATGAAGACTATGTGGATAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGtggcc                                                                                                                                         >  1:435995/1‑245 (MQ=255)
                                                                                                        tCGATCGCCTGCTTGGGCTGGAAATTGGGGCTGACGACTATGTGGCGACCCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGACCGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTAGTCGTTTTGAATTGAATGAACCCTCGGCGCAGATACGCTGGTTTGACACGCCATTATAGCTGCCGCGGTATGAGTTTTTATTGCTGAAGCCGTTACTCacatcgccggcc                                                                                                                                     >  2:170/1‑240 (MQ=255)
                                                                                                             cGCCTGCTTGGGCTGGAAATTGGTGCTGCCGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCCTTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTAATCAAGTCGCCAGGCCGCGt                                                                                                                                >  2:492715/1‑251 (MQ=255)
                                                                                                               ccTGCTTGGGCTGGAAATTGGTGCTGACGTCTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCt                                                                                                                              <  1:346616/251‑1 (MQ=255)
                                                                                                                      gggCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGGCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                                                        <  1:492715/250‑1 (MQ=255)
                                                                                                                       ggCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGc                                                                                                                      <  1:374008/251‑1 (MQ=255)
                                                                                                                           ggAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                                                                 >  1:138445/1‑236 (MQ=255)
                                                                                                                           ggAAATTGGGGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                                                                 <  2:138445/236‑1 (MQ=255)
                                                                                                                                                    ggCTAAACCGCTTTCACCCCGCGAAGTGTGCGCAAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGACTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTATGGTCCCGCCAGCAACTAATGGATAGCATGTGGGa                                                                                            >  1:181492/1‑249 (MQ=255)
                                                                                                                                                              ttttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGACTGGTCCCGCCCGCAACTAATGGATAGCGGGTGGGAAGATgcgccgg                                                                                 >  2:413571/1‑247 (MQ=255)
                                                                                                                                                                  cACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACAcc                                                                            >  1:167717/1‑251 (MQ=255)
                                                                                                                                                                              ggggcgcCAGGGTGCGCACTTTACTGCGTCGGGTGAAGAATTTCTCGTCGCCTTCCCCCGTCATCCGTATTGGTCATTTTGAATTAAATGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                                                                <  1:262215/248‑1 (MQ=255)
                                                                                                                                                                                      aGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGACGTCTCCCGTCATACGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAc                                                        >  1:197167/1‑251 (MQ=255)
                                                                                                                                                                                          tGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCTTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCACGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGAGCAGGACACCCACGAACGCACCGTTGATCCCcac                                                    >  2:419354/1‑251 (MQ=255)
                                                                                                                                                                                          tGCGCACCTTACTGCGTAGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACCCCTACGAACGCACCGTTGATACCcac                                                    >  2:326222/1‑251 (MQ=255)
                                                                                                                                                                                                    aCTGCGTCGGGTGAAGAAATTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAAcct                                          >  2:167476/1‑249 (MQ=255)
                                                                                                                                                                                                       gCGTGGGGGGAAGAAGTTCTCTTCGCCGCCCCCCTTCCTCCGTTTTGGTCATTTGAATTTGAAGGACCCCGCGGCGCAGACCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgttg                                        <  1:394652/250‑4 (MQ=255)
                                                                                                                                                                                                        cGTCGGGTGAAGAAGTTCTCGACGCAGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGATGGATTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTAGCCAGTCCGCGTCTGGTCCCGCCAGCAACTAATGGGTACCGGGGGGGAAGATGCGAAGGACCCCTAAGATAGCACAGTTGATACCCACATTAAAACGTTGCg                                      >  2:387618/1‑251 (MQ=255)
                                                                                                                                                                                                         gtggggtaaaaaagttccggaggcggTCTCCCGTCACCCGTGTGGGAAATTTTAAATGAAATGACCCCGGGGCGCAGACAACTTGGTTGGACACGCCATTAGCGCTGACTCGTGTTGAGTTTTTATTGGGGAAGCCGTTCCTCAAGTCCCCGGGCCGCGTTGGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGt                                     <  2:73296/226‑1 (MQ=255)
                                                                                                                                                                                                               caccccccttttccgatctgtcccccGTCATCCGTTTTGGTAATTTGAAATTAAATAAACCCGGGGCGCATATCAGCTGTTTTGCCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                               <  2:202321/228‑1 (MQ=255)
                                                                                                                                                                                                                               gCCGTCTCCCGTCCTCCGTATTGGTCATTTTGAATTGAATGAACCAGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGAATGAGTTTTTATTGCTGAAGACGTTACTAAAGTCGCCAGGCCGCGTATGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCCGGACACCTACGATCGACCCGTTGATCCCCACATTAAAACGTTGCGTGCCAAGCGGCGCGCCATCAAcc               >  2:448636/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                               >  1:202321/1‑232 (MQ=255)
                                                                                                                                                                                                                                            aTCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCAACATTAAAACGCTGCGTGCCAAGCTGCGCGCCCTCAACCCCGCTCTTTCAcc  >  2:450372/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACC  >  NC_000913/4636165‑4636650

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