Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,513 A→G S169S (TCA→TCG)  creB → DNA‑binding transcriptional regulator CreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,5130AG98.3% 175.8 / NA 60S169S (TCA→TCG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (38/21);  minor base T (0/1);  total (38/22)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCA‑CCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTC  >  NC_000913/4636270‑4636758
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
cgagcGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGcc                                                                                                                                                                                                                                                <  1:350398/247‑1 (MQ=255)
 gATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCcg                                                                                                                                                                                                                                               >  2:421100/1‑251 (MQ=255)
     cccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTc                                                                                                                                                                                                                                           <  1:511383/250‑1 (MQ=255)
      ccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTc                                                                                                                                                                                                                                           >  2:511383/1‑250 (MQ=255)
       cTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTg                                                                                                                                                                                                                                         >  1:204090/1‑251 (MQ=255)
           ttGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTc                                                                                                                                                                                                                                      <  2:117364/250‑1 (MQ=255)
            tGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCCGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTccc                                                                                                                                                                                                                                    >  2:303407/1‑251 (MQ=255)
             gggCTGGAAATTGGTGCTTACGACTATGTGGCTAACCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATGGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTccc                                                                                                                                                                                                                                    <  1:615555/250‑1 (MQ=255)
                    aaaTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCa                                                                                                                                                                                                                             >  1:312264/1‑250 (MQ=255)
                                    ctatgTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGtgt                                                                                                                                                                                                            >  2:149188/5‑249 (MQ=255)
                                               aaCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATgcg                                                                                                                                                                                                 >  1:329279/1‑251 (MQ=255)
                                                      ttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGac                                                                                                                                                                                           >  1:328963/1‑250 (MQ=255)
                                                           cccGCGAAGTGTGCGCCAGTGTGCCCCCTTTCCTGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCTTCATCCGTATTGGTCATTTTGAATGGAATGACCCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTGCTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTAc                                                                                                                                                                                     <  1:270189/251‑1 (MQ=255)
                                                                 aagtggcccccaggtggcgcCCCTTACTGGGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGCCACCCGATTTGGTATTTTTAAATTGAATGACCCCGGGGCCCAGATCAGTTGG‑TTTGACACGCCATTAACGCTGTCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGc                                                                                                                                                                               <  1:460151/235‑1 (MQ=255)
                                                                  aGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCCTTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACg                                                                                                                                                                                    >  2:90094/1‑245 (MQ=255)
                                                                       gcgcCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                                                                                                                         >  2:283526/1‑250 (MQ=255)
                                                                                 gcgcACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCcac                                                                                                                                                                >  2:24130/1‑250 (MQ=255)
                                                                                  cgcACCTTACTGCGTCGGGTGAAGACGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGACTGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGcc                                                                                                                                                                                                                                 >  2:604425/1‑184 (MQ=255)
                                                                                           cTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgtt                                                                                                                                                     >  2:166902/1‑249 (MQ=255)
                                                                                            tGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgttg                                                                                                                                                    >  1:95485/1‑248 (MQ=255)
                                                                                              cgtgggttgaagaatgtcgcgtcGCCGTCCCCCGTCATCCGTATTGGTTATTTTGATTTAATTGAACCCGCGGCGCAAATCAGTTGT‑TTTGCCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGTTTCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGc                                                                                                                                                   <  1:417476/228‑1 (MQ=255)
                                                                                               gTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATAAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGt                                                                                                                                                 >  2:391575/1‑251 (MQ=255)
                                                                                                     tGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATgg                                                                                                                                                                                                                     >  1:587863/1‑177 (MQ=255)
                                                                                                         gaatttCTCGTCGCCGTCTCCCTTCATCCTTATTGGTCATTTTAATTTGAATGAACCCGCGGCGCAGATCAGCGGG‑TTTAACACGCCATTACCGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTgc                                                                                                                                       <  1:397570/247‑1 (MQ=255)
                                                                                                         gaagtgccctcccccgcctcCCTTCATCCTTATCGGTTATTTTGAATGGAATGACCCCGCTGCGCAGAACA‑CTGGGTTTCACCCGCCATTAACGCTGACGCGGTATGAGTTCTTATTGCTGAAGACGTTACCCAAGCCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAACCACATTAAAACGTTGCGTGCCAAGCTgc                                                                                                                                       <  1:338264/234‑1 (MQ=255)
                                                                                                                     gCCGTCGCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCCCAGATCAGCTGG‑TTTGACACGCCATTACCGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCAACATTAAAACGTTGCGTGCCAAGCTGCGCGCAATCAAcc                                                                                                                           >  1:257142/1‑251 (MQ=255)
                                                                                                                         tctcCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCg                                                                                                                        >  1:89430/1‑250 (MQ=255)
                                                                                                                            cccGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGctct                                                                                                                    >  1:249341/1‑251 (MQ=255)
                                                                                                                            cccGTCATCCGTATGGGTCATTTTGAATTGAATGACCCCGCGGCGAAGATCAGTTGG‑TTTGACACGCCATTAACGCTGACGCGGTAGGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGGTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCt                                                                                                                    <  1:480694/251‑1 (MQ=255)
                                                                                                                                 cATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCAc                                                                                                               >  1:92914/1‑251 (MQ=255)
                                                                                                                                        aTTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGTCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGCTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTa                                                                                                         >  2:143138/1‑250 (MQ=255)
                                                                                                                                          tggtcATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCGGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAAt                                                                                                       >  1:431702/5‑250 (MQ=255)
                                                                                                                                                     aatttaatgaccctgcgtgccAGATCACCGGG‑TTTGACACGCCATTAACGCTGACCCGGTTGGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCGGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCTCAGGACACCTACGATCGCTCCTTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCTCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGc                                                                                           <  1:30458/231‑1 (MQ=255)
                                                                                                                                                     aattgaatGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAg                                                                                                                                                                                                      >  2:453737/1‑144 (MQ=255)
                                                                                                                                                         gaatgaaCCCGCGGCGAAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCAGATCTTTCACCGATTAATACTCATCGCGGCATgg                                                                                       >  1:79234/1‑251 (MQ=255)
                                                                                                                                                                 ccGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGc                                                                               >  1:130873/1‑251 (MQ=255)
                                                                                                                                                                    cggcgCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCCTGTGATATAGCCTg                                                                            >  2:346002/1‑251 (MQ=255)
                                                                                                                                                                         tagacacggtGG‑TTTGACCAGCCATTAGCGCTGCGGCGGTATGAGTTTGTATTGCGGAAGAGGTTATTCAA‑TCGGCCAGGCCGCTATTGGGACCGCCATCAACTAAGGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGAGCGCACCGTTGATACCCGCATTAAAACGGTGTGTGCCAAGCTGCGTGCCATCAGCCCCGAGCTTTCGCCGATTAATACTCATCGCGGCATGGGATATAGCCTGAggg                                                                        <  1:270683/241‑1 (MQ=255)
                                                                                                                                                                              cAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGt                                                                  >  2:154153/1‑251 (MQ=255)
                                                                                                                                                                                gCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAAt                                                                                                       >  2:43706/1‑212 (MQ=255)
                                                                                                                                                                                gCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTaa                                                                <  1:119306/251‑1 (MQ=255)
                                                                                                                                                                                  tGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGAGGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGAGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATAGCGGCATGGGATATAGCCTGAGGGGCCTGTAATg                                                              >  1:304237/1‑251 (MQ=255)
                                                                                                                                                                                    g‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTgcgcg                                                                                                                                    >  1:119331/1‑179 (MQ=255)
                                                                                                                                                                                            acGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCCGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTACTGCGTATCGGc                                                     >  1:560012/1‑251 (MQ=255)
                                                                                                                                                                                              gcctttaacgcggacgcGGTATGAGTTTTTAGTGCTAAAGACTTTACTCAAGTCG‑CCAGGCCGCGTCTGGGCCCGCCAGCAGCTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAt                                                   <  2:339769/235‑1 (MQ=255)
                                                                                                                                                                                              gCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAt                                                   >  2:502871/1‑251 (MQ=255)
                                                                                                                                                                                                    aacgcTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGgtt                                             >  1:130886/3‑251 (MQ=255)
                                                                                                                                                                                                      cgcTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGgttgc                                           >  1:107142/1‑250 (MQ=255)
                                                                                                                                                                                                          gacgcGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACTGTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCGTGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTgctgct                                       <  1:391575/247‑1 (MQ=255)
                                                                                                                                                                                                                 tATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTATGGTCCCGCCAGCAACTAATGGATAGCGAGTGGGAAGATGCGCAGGACCCCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCACGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAt                                >  1:388960/1‑251 (MQ=255)
                                                                                                                                                                                                                   tGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAttt                              >  1:361029/1‑251 (MQ=255)
                                                                                                                                                                                                                    gggTTTTTATGGCTGAAGACGTCACTCAGGTCG‑CCAGGCCGGGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTTTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAttt                              <  2:95485/248‑1 (MQ=255)
                                                                                                                                                                                                                    gAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttt                             <  1:154153/251‑1 (MQ=255)
                                                                                                                                                                                                                    gAGTTTTTATTGCTGAAGACGTTACTCAAGTAG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACAt                                                                                                                                                              >  1:184260/1‑122 (MQ=255)
                                                                                                                                                                                                                               gctgAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCGGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATGTGTTACtggtgg                   <  1:117809/248‑1 (MQ=255)
                                                                                                                                                                                                                                ctgAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcg                 >  1:510953/2‑251 (MQ=255)
                                                                                                                                                                                                                                  gAAGACGTTACTCAAGTCG‑CCGGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcg                 <  1:421100/249‑1 (MQ=255)
                                                                                                                                                                                                                                        gtaaCTCAAGGCT‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcggtggcgg          <  2:508178/247‑3 (MQ=255)
                                                                                                                                                                                                                                           aCTCAAGTCG‑CCAGGCCGCGTCTGGTCCCTCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTg      <  1:273640/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aaGTCG‑CCAGGCTGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCTGGTCACCTACGATCGCACCGTTGATTCCCACATTAAAACGTTGCGCGCCAAGCTGCGCGCCATCAACCCCGATCTTGCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTc  <  2:578310/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCA‑CCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTC  >  NC_000913/4636270‑4636758

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: ‑