Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,201 G→A L65L (TTG→TTA 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,2010GA98.5% 209.4 / ‑5.7 68L65L (TTG→TTAcreBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (33/34);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC  >  NC_000913/4635954‑4636447
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
aCGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGCGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgc                                                                                                                                                                                                                                                     >  2:190960/1‑251 (MQ=255)
 cgctatcatcgcgGTTACTGTCATGCCGTGGCGGCAAGAACAGAGGCGATTTATGCAACGGGAAACGTTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTGGCAGCAGGAAGGTTTTGCAGGCGATGTCGTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATGCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgcc                                                                                                                                                                                                                                                    <  2:151184/242‑1 (MQ=255)
 cgcgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTGGAATTATgcc                                                                                                                                                                                                                                                    <  1:190960/251‑1 (MQ=255)
     aTTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTTCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCa                                                                                                                                                                                                                                                <  1:354792/251‑1 (MQ=255)
         gcgcGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCCGCAGGCTCCCGACGTCATGATTCTCGAAGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGtt                                                                                                                                                                                                                                             >  2:197830/1‑250 (MQ=255)
          cgcgGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTAc                                                                                                                                                                                                                                           <  1:592637/251‑1 (MQ=255)
                  cTGTCATGACGTGGCGGCAATAACAGAGACGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCCAGGGAAAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGAAGTCATGATTCTCGATGTTGGCTTGCCTGATATCAGAGGCTTTGAATTATGCCGCCCGTTACTGGCGCTc                                                                                                                                                                                                                                   >  2:116370/1‑251 (MQ=255)
                    gTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTccc                                                                                                                                                                                                                                 >  2:392681/1‑250 (MQ=255)
                      cATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTTATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATc                                                                                                                                                                                                                               <  1:254670/251‑1 (MQ=255)
                                   cAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTTGCGCTCCATCCGGCGTTACCTGt                                                                                                                                                                                                                  <  1:119765/251‑1 (MQ=255)
                                   cAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                                                                                                  >  1:390622/1‑251 (MQ=255)
                                    aaTAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGATTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                                                                                                  >  1:267930/1‑250 (MQ=255)
                                       aaCAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCAGCCGGCGTTACCTGTACTg                                                                                                                                                                                                              <  1:508113/251‑1 (MQ=255)
                                            agGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCt                                                                                                                                                                                                         <  1:198970/251‑1 (MQ=255)
                                                 aTTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACgg                                                                                                                                                                                                    <  1:294104/251‑1 (MQ=255)
                                                                aaCGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAgg                                                                                                                                                                                     >  2:589089/1‑251 (MQ=255)
                                                                        ggTTAGTGGAAGATGAGGAAGGGATAGCCGACACGCTGGTCTTCATGTTGCAGCAGGAAGGTTGTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGTTTTGCCTGATATCAGCGTCTTTGAATTATGCCGCCAGGTACTGGCGCTCCATCCTGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGc                                                                                                                                                                             <  1:627172/251‑1 (MQ=255)
                                                                        ggTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGc                                                                                                                                                                             >  1:122534/1‑251 (MQ=255)
                                                                                gAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTgg                                                                                                                                                                     <  1:310030/251‑1 (MQ=255)
                                                                                 aaGATGAGCAAGGGATAGCTGACACGCAGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCTCGACGGCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTGCCGGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTggg                                                                                                                                                                    <  1:599189/251‑1 (MQ=255)
                                                                                 aaGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTggg                                                                                                                                                                    <  2:244335/251‑1 (MQ=255)
                                                                                   gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCt                                                                                                                                                                  >  1:381017/1‑251 (MQ=255)
                                                                                        gCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACATGTACTGTTCCTGACGGCCCGAAGGGAAGAGGTCGATCGCCTGCTTGGGCTggga                                                                                                                                                              >  2:417476/1‑248 (MQ=255)
                                                                                         cAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAAt                                                                                                                                                            >  2:323248/1‑251 (MQ=255)
                                                                                          aaGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAAtt                                                                                                                                                           >  2:363218/1‑251 (MQ=255)
                                                                                                 aGCCGACACGCTTGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCt                                                                                                                                                     <  2:337972/250‑1 (MQ=255)
                                                                                                    cGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGGGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGACTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgacg                                                                                                                                                 >  2:397570/1‑251 (MQ=255)
                                                                                                       cacGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACt                                                                                                                                              <  1:132115/251‑1 (MQ=255)
                                                                                                             gttCTACATGTTGCAGCAGGAAGGTTTGGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCGCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGTGCTGGAAATTTGTGCTGACGACTATGTgg                                                                                                                                        <  2:566822/249‑1 (MQ=255)
                                                                                                                    aTGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGCCGACTATGTGGCTAAAcc                                                                                                                                 >  2:135494/1‑251 (MQ=255)
                                                                                                                      gTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGt                                                                                                                               <  1:293437/251‑1 (MQ=255)
                                                                                                                                          ccGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGAAGGCTCCCGCCGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGAACTGTTCCTGACGGCCCGAAGGGAAGAGGGCGATCGCCTGCTTGGGCTGGAAATTTGTGCTGACGACTATGTGGCTAAACGGCTTTCACCCCGCGCAGTGTgc                                                                                                           >  2:338264/1‑251 (MQ=255)
                                                                                                                                            gTCGATGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGGGAAGTGTgcgc                                                                                                         <  2:528296/251‑1 (MQ=255)
                                                                                                                                               gAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                                                                                                        >  1:339769/1‑249 (MQ=255)
                                                                                                                                                      ttGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTgcg                                                                                                <  2:441699/250‑1 (MQ=255)
                                                                                                                                                        gAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCAc                                                                                             <  2:122534/251‑1 (MQ=255)
                                                                                                                                                           cGAGGCTTGCCGGTGCAGGATAAAGCTCGCCAGCAGGCGCCCGACGTCATGATTATCGATGTTTGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCAACtt                                                                                          <  2:547477/251‑1 (MQ=255)
                                                                                                                                                            gAGGCTTGCCGGTGCTGGCTAACGCTCGCCAGCAGGCTCCCGCCGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCCCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGCAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGCGCGCACCtt                                                                                          >  2:117809/1‑250 (MQ=255)
                                                                                                                                                                cTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGTTGCTGACGACTACGTGGCTAAACCGGTTTCAACCAGCGAAGTGTGGGCAAGGGAGAGCACCTTACTGc                                                                                     >  1:385512/1‑251 (MQ=255)
                                                                                                                                                                 ttGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCCGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCATGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTACACCGTTTTCACCCCGCGAAGTGTGCGCCAg                                                                                                      >  2:227076/1‑233 (MQ=255)
                                                                                                                                                                 ttGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAg                                                                                                      <  1:227076/233‑1 (MQ=255)
                                                                                                                                                                     cGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTACACCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGCCTATGTGGCTAAACCGTTTTCACCCAGCGAAGTCGGCGCCAGGGGGCGCACCTTACTGCGTCgg                                                                                >  1:537035/1‑251 (MQ=255)
                                                                                                                                                                          cTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTga                                                                            >  1:37283/1‑250 (MQ=255)
                                                                                                                                                                          cTGGATAAAGCTCGCCAGCAGGCTCCCGACGTAATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTga                                                                            >  1:508178/1‑250 (MQ=255)
                                                                                                                                                                                          agcaggctcccgccgtactGATTCTAGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCg                                                           >  2:615555/19‑251 (MQ=255)
                                                                                                                                                                                          agcaggctcccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCg                                                           <  1:469524/243‑1 (MQ=255)
                                                                                                                                                                                            caggctcccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGtc                                                         <  1:392681/245‑1 (MQ=255)
                                                                                                                                                                                                    cGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc                                                     >  2:334560/1‑247 (MQ=255)
                                                                                                                                                                                                    cGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc                                                     <  1:334560/247‑1 (MQ=255)
                                                                                                                                                                                                         tCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCgg                                            >  2:150659/1‑250 (MQ=255)
                                                                                                                                                                                                          cATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGCAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCCGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGt                                            >  2:577571/1‑250 (MQ=255)
                                                                                                                                                                                                            tGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTc                                                                                  <  2:64471/210‑1 (MQ=255)
                                                                                                                                                                                                            tGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTc                                                                                  >  1:64471/1‑210 (MQ=255)
                                                                                                                                                                                                             gATTCTCGATGTTGGTTTGCCTGATAGCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCTCCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCAGCCGTAtt                                         <  2:37283/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggtttGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTTGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGGATCTTACTGCGTCGGGTGAAGAAGTGCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaat                            <  1:254314/246‑1 (MQ=255)
                                                                                                                                                                                                                          ggtttGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaatt                           <  1:135494/247‑1 (MQ=255)
                                                                                                                                                                                                                          ggtttGCCTGATATCAGCGGCTTTGAATTATGACGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGACGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaatt                           >  1:529579/5‑251 (MQ=255)
                                                                                                                                                                                                                              tGCCTGATATCAGCGGCTTTGAGTTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATAGCCTGCTTGGGCTGGTAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCACGGTGCGCACATTACTGCGGCGGGAGACGAAGGTCTCGTCGCCGTCCGCCGTCATCCGTATTGGTCATTttgaattgaa                        >  2:380348/1‑250 (MQ=255)
                                                                                                                                                                                                                              tGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaattgaa                        <  1:7405/250‑1 (MQ=255)
                                                                                                                                                                                                                              tGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTAACGGCCCGAAGTGCAGAGGTCGATCGCCTGCTTGGGCTGGCAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCCCCTTACTACGTCGGGTGAAGAAGTTCTAGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTtgcattgaat                       >  1:610973/1‑251 (MQ=255)
                                                                                                                                                                                                                               gcctgATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATtgaatg                      >  2:478890/5‑251 (MQ=255)
                                                                                                                                                                                                                               gcctgATATCAGCGGCTGTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAAGTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAATAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTtgaattgaat                       <  2:373643/246‑1 (MQ=255)
                                                                                                                                                                                                                                cctgATATCAGCGGCTTTGAATTATGCCGCCAGTTCCTGGCGCTCCATCCGGCGTTCCCTGTACTGTTCCTGACGGCCAGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATtgaatga                     <  2:331832/248‑1 (MQ=255)
                                                                                                                                                                                                                                    atatCAGCGGTTTTGAATTATGCCGCCAGTTACTGGCGCTCCAGCCGGCGTTACCTGTACTGTGCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAATCCGGTTTCACCCCTCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCGCCCGTCAGCCGTATGGGTCATTTTCAATTGAATGAAccc                 <  2:536295/251‑1 (MQ=255)
                                                                                                                                                                                                                                        caGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCgcgg             >  1:220967/2‑251 (MQ=255)
                                                                                                                                                                                                                                         aGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAACTTGGTGCTGACGCCTATGTGGCTAAACCGTTTTCACCCCGCGCAGTGTGCGCCAGGGTGAGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGGTGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATAAAcccacgac            >  2:337207/1‑246 (MQ=255)
                                                                                                                                                                                                                                            ggCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGTCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGgcgc          >  1:578310/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   aaTTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGc  <  2:376213/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
ACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC  >  NC_000913/4635954‑4636447

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