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,3850AT94.8% 147.5 / ‑0.6 58T127S (ACG→TCG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (28/27);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

CAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCG  >  NC_000913/4636138‑4636614
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
ccagcaggctcccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGCAGTTCTCGTCGc                                                                                                                                                                                                                                     >  1:155920/11‑251 (MQ=255)
   gcaggctcccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGt                                                                                                                                                                                                                                  >  1:46505/8‑251 (MQ=255)
             gACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGCCTATGTGGCTAAACCGATTTCACCCCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGCGCAGAAGGTCTCGTCGCCGTATCCCGACAt                                                                                                                                                                                                                        >  2:155610/1‑251 (MQ=255)
             gACGTCATGATTCTCGAGGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCAt                                                                                                                                                                                                                        <  2:394211/251‑1 (MQ=255)
                gTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCg                                                                                                                                                                                                                     >  1:172694/1‑251 (MQ=255)
                     gatgctcGATGTTGGTTTTCCTGATATCAGCGGCTTTTAATTATGCCGCCAGTTACTGGCGCTCCAGCCGGCGTTACCTGTACGGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTAtt                                                                                                                                                                                                                 <  2:291267/246‑1 (MQ=255)
                     gATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTAtt                                                                                                                                                                                                                 <  1:80121/250‑1 (MQ=255)
                     gATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTg                                                                                                                                                                                                                >  1:360972/1‑251 (MQ=255)
                           ctatgttgtgttgccttatatCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCACCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTggtc                                                                                                                                                                                                             <  2:326079/231‑3 (MQ=255)
                           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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
CAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCG  >  NC_000913/4636138‑4636614

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