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,1390AC91.0% 177.7 / 2.7 67K45Q (AAG→CAG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base C (34/27);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTT  >  NC_000913/4635895‑4636380
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
acgcaacgcGCTCGCGTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGGGGCGGCAATAACAGAGGCGATGTATGCAACGGGAAACGGTCTGGTTAGTTGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAgg                                                                                                                                                                                                                                              <  1:443010/251‑1 (MQ=255)
    aacgcGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGGGCTGGATAAAGCTCGCAAGCAGGTAcc                                                                                                                                                                                                                                          >  2:196431/1‑251 (MQ=255)
     acgcGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAggctccc                                                                                                                                                                                                                                         <  2:33296/251‑6 (MQ=255)
       gcgcTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGGGCTGGATAAAGCTCGCCAGCAggctcccg                                                                                                                                                                                                                                        >  1:302681/1‑244 (MQ=255)
          ctcgctTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGt                                                                                                                                                                                                                                    >  1:360979/1‑251 (MQ=255)
                         tACTCCGACAAAGTTGAAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGAAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGAGGTCAGGATTCTCGATgg                                                                                                                                                                                                                     >  2:103301/1‑250 (MQ=255)
                                aCAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGtttg                                                                                                                                                                                                              >  1:447001/1‑248 (MQ=255)
                                aCAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGATCCTGCCATGCCGTGGCGGCAATAACAGAGGCGATTGATGCAACGGGAAACGGTCTGGTTATTGGAAGATGAGCAAGTGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGATGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGtttg                                                                                                                                                                                                              <  2:94436/251‑4 (MQ=255)
                                   aaGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGGCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCGGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGcct                                                                                                                                                                                                           <  1:32489/251‑2 (MQ=255)
                                                        aaaaCGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAAtt                                                                                                                                                                                       >  1:291267/1‑250 (MQ=255)
                                                              ggATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTATATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgccgc                                                                                                                                                                               <  2:79231/250‑1 (MQ=255)
                                                                 ttAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgccgcccg                                                                                                                                                                            >  2:480165/1‑249 (MQ=255)
                                                                   aGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCg                                                                                                                                                                                                 <  2:428645/229‑1 (MQ=255)
                                                                   aGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCg                                                                                                                                                                                                 >  1:428645/1‑229 (MQ=255)
                                                                            ccTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTAGGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCt                                                                                                                                                                 <  2:360979/251‑1 (MQ=255)
                                                                             cTGTCCTGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTCGAATTATGCCGCCAGTTACTGGCGCTc                                                                                                                                                                >  2:346616/1‑251 (MQ=255)
                                                                             cTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGGTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATTCCGCCAGTTAATGGCGCTc                                                                                                                                                                >  2:263977/1‑251 (MQ=255)
                                                                              tGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTcc                                                                                                                                                               >  1:306133/1‑251 (MQ=255)
                                                                                tCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGACTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCAt                                                                                                                                                             >  2:403827/1‑251 (MQ=255)
                                                                                 cATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGTCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATc                                                                                                                                                            >  1:74350/1‑251 (MQ=255)
                                                                                    gCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCgg                                                                                                                                                         <  2:470257/251‑1 (MQ=255)
                                                                                         ggcggcAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTa                                                                                                                                                    <  2:423365/251‑1 (MQ=255)
                                                                                              cAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                               >  1:347478/1‑251 (MQ=255)
                                                                                               aaTAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGATCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                               >  1:255824/1‑250 (MQ=255)
                                                                                                    cAGAGGCGATTTATGCAACGGGAAACGGTCTGGTGAGTGGCAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTTCAGAAGGAAGGTTTTGACGGCGCCGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGATCCCGACGACATGATTCTCGCTGTTGGTTTGCCTGCAATCAGCGGCTTTGCATTAGTAAGCCAGTTATTGGCGCTCCATCCgtggttaacgtaacggtt                                                                                                                                         >  2:409656/1‑234 (MQ=255)
                                                                                                        ggCGATATATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTg                                                                                                                                     <  1:29865/251‑1 (MQ=255)
                                                                                                                  gCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCCTCCGGCGTTACCTGTACTGTTCCTGACGGCCCggg                                                                                                                           >  2:115016/1‑249 (MQ=255)
                                                                                                                  gCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGaa                                                                                                                           <  1:263977/251‑1 (MQ=255)
                                                                                                                   cAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGACTACATGTTGAAGAAGGAAGGTTTTGCCGTCGAGGACTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCCGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTCGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGGTACCTGTACTGTTCCTGACGGCCCGa                                                                                                                            >  2:464420/1‑249 (MQ=255)
                                                                                                                           aaCGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAgg                                                                                                                  >  1:451005/1‑251 (MQ=255)
                                                                                                                               gTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCg                                                                                                               >  1:394211/1‑250 (MQ=255)
                                                                                                                                       aGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATTGTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTg                                                                                                       <  2:95761/250‑1 (MQ=255)
                                                                                                                                        gTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCt                                                                                                     <  1:243195/251‑1 (MQ=255)
                                                                                                                                                tGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGACTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTTGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGGCGGCCCGAAGTGAAGAGGTCGACCGCCTTCTTGGGCTgg                                                                                             >  2:15205/1‑251 (MQ=255)
                                                                                                                                                      aGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCATGCTCCCGACGTCATGATTCTCGATGTTGGTTTGTCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCa                                                                                                                                                              <  1:236647/180‑1 (MQ=255)
                                                                                                                                                            aGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgcc                                                                                                                                                                                 <  1:418140/155‑1 (MQ=255)
                                                                                                                                                            aGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgcc                                                                                                                                                                                 >  2:418140/1‑155 (MQ=255)
                                                                                                                                                              ccGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                               <  1:403827/251‑1 (MQ=255)
                                                                                                                                                               cGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgacg                                                                              >  1:377718/1‑251 (MQ=255)
                                                                                                                                                                       tGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTAtgtg                                                                      >  1:264949/1‑251 (MQ=255)
                                                                                                                                                                            tACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTa                                                                  >  1:284138/1‑250 (MQ=255)
                                                                                                                                                                                         agGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCGAAACCGTTTTCAccc                                                    >  2:394652/1‑251 (MQ=255)
                                                                                                                                                                                         agGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCAccc                                                    >  2:16259/1‑251 (MQ=255)
                                                                                                                                                                                          gtaaGGTTTTGCCGTCGATGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCAcccc                                                   <  1:127509/249‑1 (MQ=255)
                                                                                                                                                                                           gAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAACGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGACGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTCAACCGTTTTCAccccc                                                  >  2:476576/1‑250 (MQ=255)
                                                                                                                                                                                                   tGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAgtgt                                          >  2:252816/1‑251 (MQ=255)
                                                                                                                                                                                                      cGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcg                                       <  1:480165/251‑1 (MQ=255)
                                                                                                                                                                                                      cGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGGCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCGCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcg                                       <  2:53615/251‑1 (MQ=255)
                                                                                                                                                                                                       gTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                                      >  1:326079/1‑251 (MQ=255)
                                                                                                                                                                                                          gAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                                     >  1:342039/1‑249 (MQ=255)
                                                                                                                                                                                                             gTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGGGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGGGCGCCAgggg                                >  2:265946/1‑250 (MQ=255)
                                                                                                                                                                                                              tCTTTGAGCGAGGCTTGCCGGTGCTGGATGAAGCTAGCCAGCAGGGTCCGGACGTCATGATTGGCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGTTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTCAAGAGGTCGCTCGCCTGCTTGGGCTGGAAATTGGTGCTGTCGACTATATGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAggggg                               >  1:203169/1‑249 (MQ=255)
                                                                                                                                                                                                                 ttGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTgcg                             <  1:250559/250‑1 (MQ=255)
                                                                                                                                                                                                                     gCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCTATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCt                        <  2:209685/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCCCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTATGAATTATGCCGCCAGTTTCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGTCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCt                        <  1:488818/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGTCGTCATGATTCTCGAGGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGGTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCtt                       <  2:451005/251‑1 (MQ=255)
                                                                                                                                                                                                                        aGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTAc                     >  1:362677/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTg                   <  1:252816/251‑1 (MQ=255)
                                                                                                                                                                                                                           cgtGCCGGTGCTGGATAAAGCCCGCCAGCAGGCTCCCGACGTCAGGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTGTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGACCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTACACCCCGCGAAGTGTGCGCCAGTGTGCGCACCTTACTGc                  <  2:255824/249‑1 (MQ=255)
                                                                                                                                                                                                                           cTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc                  >  1:479954/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGCCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTc               <  1:196431/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAGCTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGCAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCggg            >  2:150218/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                               <  2:127312/179‑1 (MQ=255)
                                                                                                                                                                                                                                      tGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                               >  1:127312/1‑179 (MQ=255)
                                                                                                                                                                                                                                         aTAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc                  >  2:109297/1‑237 (MQ=255)
                                                                                                                                                                                                                                         aTAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc                  <  1:109297/237‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaGCTCGCCAGCAGGCTCCCGTCGTCAGGATTCGCGATGTTGGTTTGCCTGATACCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCAGCCTGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGTCTGCTTGGGCTGTAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCggttcaataagtt  <  2:394031/251‑12 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTT  >  NC_000913/4635895‑4636380

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