Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,146 2 bp→CT coding (140‑141/690 nt) 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,1460TC87.2% 77.9 / 2.7 39V47A (GTA→GCT) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (13/21);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
*NC_0009134,636,1470AT86.8% 74.0 / 1.8 39V47A (GTA→GCT) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (12/21);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTC  >  NC_000913/4635916‑4636383
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                            
ggCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATtctc                                                                                                                                                                                                                            <  1:165746/250‑1 (MQ=255)
      cgccGATAAAGTGGGAGAAGGATCGCCGATAAACTAGATTAGCGCGGTTCTTGTCATGCCGTGGCGTCAATAACAGAGGCGATTTATGCAACGGGAACCGGTCTGGTTATTGGAAGATGAGGAAGGGATAGCCGACACGTTGGCCTACATGTTGCAGCTGGAAGGTTTTGCCGTCG‑AGGTTTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTg                                                                                                                                                                                                                     <  2:219543/249‑1 (MQ=255)
        ccgataatgtgtgaGAAGGTTCGCCGAAAAACGCGATGAGCGCGGTTCCTGTCATGCCGTTGCGGCCATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCTAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCGCGCCAGCATG‑CTCCCGACGTCATGATTCTCGATGtt                                                                                                                                                                                                                      <  2:203734/235‑1 (MQ=255)
               aGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATAATTCTCGATGTTGGTTTGcc                                                                                                                                                                                                              >  1:306760/1‑249 (MQ=255)
               aGTTGGAGAAGTTGCGACGAAAAACGCTATGAGCCCTGTTCCTGTCATGCCGTGGCGGCAATAACATAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTGTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATAATTCTCGATGTTGGTTTGcc                                                                                                                                                                                                              <  2:306760/249‑1 (MQ=255)
                  tGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGat                                                                                                                                                                                                          <  1:62396/250‑1 (MQ=255)
                   gTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGGCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGata                                                                                                                                                                                                         <  1:199036/250‑1 (MQ=255)
                             cGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCttt                                                                                                                                                                                              <  1:274046/251‑1 (MQ=255)
                              gCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGAAATAAAAGAGGCGATTTATGCCACGGGAAACGGTCTGGGTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCAACATGTTGCAGCAGGAAGGTTTTGCAGTCGGAGG‑CTTTGAGCGCGGATTGCCGGTTCTGGGTAAAGCTCGCACGCCGG‑TCCCCGGCGTTCCGATTCTCGATGGTGGTTTGCCTGatataacaggctttg                                                                                                                                                                                             >  2:380193/1‑240 (MQ=255)
                                    aaaCGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑CGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGCCGTCATGATTCTCGATGTTGGTTTGCATGATATCAGCGGCTTTGAAttat                                                                                                                                                                                       >  1:169072/1‑249 (MQ=255)
                                         cgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTTCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTgccgc                                                                                                                                                                                  >  1:183358/1‑251 (MQ=255)
                                         cgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgccgc                                                                                                                                                                                  <  1:128164/251‑1 (MQ=255)
                                              aGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCCTGATTCTCGATGTTGGTTTGCCTGCTATCAGCGGCTTTGAATTATGCCGCCGGtt                                                                                                                                                                             >  2:181600/1‑251 (MQ=255)
                                                             aTGCCGTGGCGGCATAAACAGAGGCGATTTAGGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTTCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATcc                                                                                                                                                              <  2:183358/251‑1 (MQ=255)
                                                             aTGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATcc                                                                                                                                                              >  1:108331/1‑251 (MQ=255)
                                                                   tggcgcaaaTAACAGAGGCGATTTATGCAACGGGACAAGGCATGGTTAGAGGAAGATGAGCAAGGGATAGACGACACGCTGGTCTCCATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGAGAGGATTGCCGGTGCTGGATAAAGCTCGCAATCAGG‑CACCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGCATTATGCCCACTGTTCCTGGCGCTCCATCCGGCGtt                                                                                                                                                        >  2:163376/8‑251 (MQ=255)
                                                                                  agGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCt                                                                                                                                         <  1:20254/251‑1 (MQ=255)
                                                                                  agGCGATTTATGCAACGGGAAACGGTCTGGATAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCt                                                                                                                                         >  1:99509/1‑251 (MQ=255)
                                                                                               aacggggaccgttCTGGTTAGTGGAAGATGAGGAAGGGATAGCCGACACGCTGGGCTACATTTTGCAGCAGGAAGGTGTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGt                                                                                                                            <  2:181463/239‑1 (MQ=255)
                                                                                                        tggTCTGGTTAGTGGAAGATGAGAAAGGGATAGCCGACACGCTGGTCTACATGGTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGCGCTTGATAAAGCTCGCCAGCAGG‑CTCCCGACGGCATGATTGTCGATGTTGGTATGCCTGACATCAGCGGCTTTGAATACTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGAACTGTTCCTGACGGCCCGactcgaaggagtc                                                                                                                   >  2:376864/2‑239 (MQ=255)
                                                                                                           gcTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCGGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGAt                                                                                                                <  1:135023/250‑1 (MQ=255)
                                                                                                               gTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGcc                                                                                                            >  1:127271/1‑251 (MQ=255)
                                                                                                                        aGATGAGCTAGGAATAGGGGACACGCTAGGATACATGTTGCAGCAGGAAGGTTTTGACGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGATACA‑GACGTCATGAATCGCGATGTTGGTTTGCTAGATATCAGCGGCTTTACATTATGCAGCCAGTTACTGGAGATCCATACGGCGTTACCTGAACTGATCCTGACGGCCCGAAGGGGAGAGGTCGATCGCCTGCCTGGGc                                                                                                   >  2:58798/1‑251 (MQ=255)
                                                                                                                          agtaGCAAGGGATAGCCGACACGCTGGTCTACATGTGGGTGTAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCGTGATAAAGCTCGCCAGCAGG‑CTCCCGGCGTCATGATTCTCGATGATGGTTTGCCTGATACCAGCGGCTTTGAATTATGCCGCTAGTTGCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGTCTCGAAGTGAAGAGGTCCATCGCCTGCTTGGGcgg                                                                                                 <  2:33071/248‑3 (MQ=255)
                                                                                                                             ggCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTTCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGaa                                                                                              <  1:279421/250‑1 (MQ=255)
                                                                                                                                aaGGGATAGCCGCCACGATGGTCTACATGTTGCAGCAGGAAGGATTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGTCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGAAGGCCCGAAGTGAAGAGGTAGATCGCCTGCTTGGGCTGGAAAtt                                                                                           >  2:289076/1‑251 (MQ=255)
                                                                                                                                    gATAGCCGACACGCTGGGCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGATTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTg                                                                                       <  1:240390/251‑1 (MQ=255)
                                                                                                                                     aTAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGc                                                                                      <  1:67820/251‑1 (MQ=255)
                                                                                                                                       aGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTg                                                                                    <  2:150012/251‑1 (MQ=255)
                                                                                                                                                        aCATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGACTGATATCAGCGGCTTTGAATTATGCCGCAAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGCGACGAGGCCGATCGCCTGCTTGGGCTGGAAATTTGTGCGGACGACTATGTGGCTaa                                                                    >  2:105500/1‑250 (MQ=255)
                                                                                                                                                                     ggAAGGCTTCGCAGTCG‑AGGTCTTTGAGCGCGGCTTGCAGGTCCTGGATAAAGCTCGCAAGCAGG‑CTCCCGACGTCCTGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGATTTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGATACTGACGGCACGAAGTGAAGCGGTCGATCGCCTGCTTGGGATGGCACTTGGTGCTGACGACTATGTGGCAAAACCGTTTTCAcccc                                                      >  2:351702/1‑251 (MQ=255)
                                                                                                                                                                                 cGTCG‑AGGTCTTTGCGCTAGGCTTGCCGGTGCTGGATAAAGCACGCCAGCAGG‑CTCCCGACATGATGATTCTCGATGTTGTTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTTCTGGTGCTCCATCCGGTGTTACCTGCACTTTTCCTGACGGCCCGAAGTGAAGAGGTCGATCTCCTGCTTGGGCTGGAAATTGTGGCTGATGTCTATGTGGCTAAAGCGTTTTCatctcgattagtgtgc                                           <  2:332323/250‑16 (MQ=255)
                                                                                                                                                                                  gTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                                         <  1:380193/251‑1 (MQ=255)
                                                                                                                                                                                        ggTCATTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGAGGCTTTGAATTATGCCGCCAGTTACTGGCGCCCCATACGGCGTTACCTGTACTTTTCCTTACGCTCCAATGTGTAGATGGCGATCGCCTGCTTGGGCTGGAAATAGGTGCTAACCAGTATGTATCTACACCGTTTTCAcccaaagaactgctccccggtg                                    >  2:118215/1‑232 (MQ=255)
                                                                                                                                                                                               gAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCa                              <  1:169337/250‑1 (MQ=255)
                                                                                                                                                                                                      gCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTg                      >  1:279904/1‑251 (MQ=255)
                                                                                                                                                                                                                    gATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTgaagaa        <  1:143261/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTGCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTtc    <  1:337170/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTtctc  <  1:373323/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                            
GGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCG‑AGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTC  >  NC_000913/4635916‑4636383

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