Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,139 A→C K45Q (AAG→CAG)  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,1390AC85.7% 98.6 / 1.6 42K45Q (AAG→CAG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base C (13/23);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTC  >  NC_000913/4635896‑4636383
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
cgcaacgcGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAGGCTCGCCAGCAggc                                                                                                                                                                                                                                                >  1:85261/1‑250 (MQ=255)
    acgcGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAggctccc                                                                                                                                                                                                                                            <  1:309708/251‑6 (MQ=255)
    acgcGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAggctccc                                                                                                                                                                                                                                            <  1:152241/251‑6 (MQ=255)
                    ggCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATtctc                                                                                                                                                                                                                            <  1:165746/250‑1 (MQ=255)
                          cgccGATAAAGTGGGAGAAGGATCGCCGATAAACTAGATTAGCGCGGTTCTTGTCATGCCGTGGCGTCAATAACAGAGGCGATTTATGCAACGGGAACCGGTCTGGTTATTGGAAGATGAGGAAGGGATAGCCGACACGTTGGCCTACATGTTGCAGCTGGAAGGTTTTGCCGTCGAGGTTTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTg                                                                                                                                                                                                                     <  2:219543/249‑1 (MQ=255)
                            ccgataatgtgtgaGAAGGTTCGCCGAAAAACGCGATGAGCGCGGTTCCTGTCATGCCGTTGCGGCCATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCTAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCGCGCCAGCATGCTCCCGACGTCATGATTCTCGATGtt                                                                                                                                                                                                                      <  2:203734/235‑1 (MQ=255)
                                   aGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATAATTCTCGATGTTGGTTTGcc                                                                                                                                                                                                              >  1:306760/1‑249 (MQ=255)
                                   aGTTGGAGAAGTTGCGACGAAAAACGCTATGAGCCCTGTTCCTGTCATGCCGTGGCGGCAATAACATAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTGTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATAATTCTCGATGTTGGTTTGcc                                                                                                                                                                                                              <  2:306760/249‑1 (MQ=255)
                                      tGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGat                                                                                                                                                                                                          <  1:62396/250‑1 (MQ=255)
                                       gTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGGCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGata                                                                                                                                                                                                         <  1:199036/250‑1 (MQ=255)
                                                 cGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCttt                                                                                                                                                                                              <  1:274046/251‑1 (MQ=255)
                                                  gCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGAAATAAAAGAGGCGATTTATGCCACGGGAAACGGTCTGGGTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCAACATGTTGCAGCAGGAAGGTTTTGCAGTCGGAGGCTTTGAGCGCGGATTGCCGGTTCTGGGTAAAGCTCGCACGCCGGTCCCCGGCGTTCCGATTCTCGATGGTGGTTTGCCTGatataacaggctttg                                                                                                                                                                                             >  2:380193/1‑240 (MQ=255)
                                                        aaaCGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGCGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGCCGTCATGATTCTCGATGTTGGTTTGCATGATATCAGCGGCTTTGAAttat                                                                                                                                                                                       >  1:169072/1‑249 (MQ=255)
                                                             cgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTTCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTgccgc                                                                                                                                                                                  >  1:183358/1‑251 (MQ=255)
                                                             cgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgccgc                                                                                                                                                                                  <  1:128164/251‑1 (MQ=255)
                                                                  aGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCCTGATTCTCGATGTTGGTTTGCCTGCTATCAGCGGCTTTGAATTATGCCGCCGGtt                                                                                                                                                                             >  2:181600/1‑251 (MQ=255)
                                                                                 aTGCCGTGGCGGCATAAACAGAGGCGATTTAGGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTTCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATcc                                                                                                                                                              <  2:183358/251‑1 (MQ=255)
                                                                                 aTGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATcc                                                                                                                                                              >  1:108331/1‑251 (MQ=255)
                                                                                       tggcgcaaaTAACAGAGGCGATTTATGCAACGGGACAAGGCATGGTTAGAGGAAGATGAGCAAGGGATAGACGACACGCTGGTCTCCATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGAGAGGATTGCCGGTGCTGGATAAAGCTCGCAATCAGGCACCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGCATTATGCCCACTGTTCCTGGCGCTCCATCCGGCGtt                                                                                                                                                        >  2:163376/8‑251 (MQ=255)
                                                                                                      agGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCt                                                                                                                                         <  1:20254/251‑1 (MQ=255)
                                                                                                      agGCGATTTATGCAACGGGAAACGGTCTGGATAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCt                                                                                                                                         >  1:99509/1‑251 (MQ=255)
                                                                                                                   aacggggaccgttCTGGTTAGTGGAAGATGAGGAAGGGATAGCCGACACGCTGGGCTACATTTTGCAGCAGGAAGGTGTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGt                                                                                                                            <  2:181463/239‑1 (MQ=255)
                                                                                                                            tggTCTGGTTAGTGGAAGATGAGAAAGGGATAGCCGACACGCTGGTCTACATGGTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGCGCTTGATAAAGCTCGCCAGCAGGCTCCCGACGGCATGATTGTCGATGTTGGTATGCCTGACATCAGCGGCTTTGAATACTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGAACTGTTCCTGACGGCCCGactcgaaggagtc                                                                                                                   >  2:376864/2‑239 (MQ=255)
                                                                                                                               gcTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCGGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGAt                                                                                                                <  1:135023/250‑1 (MQ=255)
                                                                                                                                   gTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGcc                                                                                                            >  1:127271/1‑251 (MQ=255)
                                                                                                                                            aGATGAGCTAGGAATAGGGGACACGCTAGGATACATGTTGCAGCAGGAAGGTTTTGACGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGATACAGACGTCATGAATCGCGATGTTGGTTTGCTAGATATCAGCGGCTTTACATTATGCAGCCAGTTACTGGAGATCCATACGGCGTTACCTGAACTGATCCTGACGGCCCGAAGGGGAGAGGTCGATCGCCTGCCTGGGc                                                                                                   >  2:58798/1‑251 (MQ=255)
                                                                                                                                              agtaGCAAGGGATAGCCGACACGCTGGTCTACATGTGGGTGTAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCGTGATAAAGCTCGCCAGCAGGCTCCCGGCGTCATGATTCTCGATGATGGTTTGCCTGATACCAGCGGCTTTGAATTATGCCGCTAGTTGCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGTCTCGAAGTGAAGAGGTCCATCGCCTGCTTGGGcgg                                                                                                 <  2:33071/248‑3 (MQ=255)
                                                                                                                                                 ggCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTTCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGaa                                                                                              <  1:279421/250‑1 (MQ=255)
                                                                                                                                                    aaGGGATAGCCGCCACGATGGTCTACATGTTGCAGCAGGAAGGATTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGTCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGAAGGCCCGAAGTGAAGAGGTAGATCGCCTGCTTGGGCTGGAAAtt                                                                                           >  2:289076/1‑251 (MQ=255)
                                                                                                                                                        gATAGCCGACACGCTGGGCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGATTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTg                                                                                       <  1:240390/251‑1 (MQ=255)
                                                                                                                                                         aTAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGc                                                                                      <  1:67820/251‑1 (MQ=255)
                                                                                                                                                           aGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTg                                                                                    <  2:150012/251‑1 (MQ=255)
                                                                                                                                                                            aCATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGACTGATATCAGCGGCTTTGAATTATGCCGCAAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGCGACGAGGCCGATCGCCTGCTTGGGCTGGAAATTTGTGCGGACGACTATGTGGCTaa                                                                    >  2:105500/1‑250 (MQ=255)
                                                                                                                                                                                         ggAAGGCTTCGCAGTCGAGGTCTTTGAGCGCGGCTTGCAGGTCCTGGATAAAGCTCGCAAGCAGGCTCCCGACGTCCTGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGATTTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGATACTGACGGCACGAAGTGAAGCGGTCGATCGCCTGCTTGGGATGGCACTTGGTGCTGACGACTATGTGGCAAAACCGTTTTCAcccc                                                      >  2:351702/1‑251 (MQ=255)
                                                                                                                                                                                                     cGTCGAGGTCTTTGCGCTAGGCTTGCCGGTGCTGGATAAAGCACGCCAGCAGGCTCCCGACATGATGATTCTCGATGTTGTTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTTCTGGTGCTCCATCCGGTGTTACCTGCACTTTTCCTGACGGCCCGAAGTGAAGAGGTCGATCTCCTGCTTGGGCTGGAAATTGTGGCTGATGTCTATGTGGCTAAAGCGTTTTCatctcgattagtgtgc                                           <  2:332323/250‑16 (MQ=255)
                                                                                                                                                                                                      gTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                                         <  1:380193/251‑1 (MQ=255)
                                                                                                                                                                                                            ggTCATTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGAGGCTTTGAATTATGCCGCCAGTTACTGGCGCCCCATACGGCGTTACCTGTACTTTTCCTTACGCTCCAATGTGTAGATGGCGATCGCCTGCTTGGGCTGGAAATAGGTGCTAACCAGTATGTATCTACACCGTTTTCAcccaaagaactgctccccggtg                                    >  2:118215/1‑232 (MQ=255)
                                                                                                                                                                                                                   gAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCa                              <  1:169337/250‑1 (MQ=255)
                                                                                                                                                                                                                          gCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTg                      >  1:279904/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTgaagaa        <  1:143261/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTGCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTtc    <  1:337170/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTtctc  <  1:373323/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
CGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTC  >  NC_000913/4635896‑4636383

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