Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,610 C→T L202L (CTG→TTG)  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,6100CT93.2% 121.1 / ‑1.5 44L202L (CTG→TTG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (24/17);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAGCGCTGACTCG‑GTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGA‑TGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCG  >  NC_000913/4636364‑4636854
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
cGTCGGGTGAAGAAGTTCTCGACGCAGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGAT‑GGATTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTAGCCAGTCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGGTACCGGGGGGGAAGATGCGAAGGACCCCTAAGATAGCACAGTTGATACCCACATTAAAACGTTGCg                                                                                                                                                                                                                                                  >  2:387618/1‑251 (MQ=255)
 gtggggtaaaaaagttccggaggcggTCTCCCGTCACCCGTGTGGGAAATTTTAAATGAAATGACCCCGGGGCGCAGA‑CAACTTGGTTGGACACGCCATTAGCGCTGACTCGTGT‑TGAGTTTTTATTGGGGAAGCCGTTCCTCAAGTCCCCGGGCCGCGTTGGGTCCCGCCAGCAACTGA‑TGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGt                                                                                                                                                                                                                                                 <  2:73296/226‑1 (MQ=255)
       caccccccttttccgatctgtcccccGTCATCCGTTTTGGTAATTTGAAATTAAATAAACCCGGGGCGCATATCAGCT‑GTTTTGCCACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                                                                                                                                                                                                                                           <  2:202321/228‑1 (MQ=255)
                       gCCGTCTCCCGTCCTCCGTATTGGTCATTTTGAATTGAATGAACCAGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GAATGAGTTTTTATTGCTGAAGACGTTACTAAAGTCGCCAGGCCGCGTATGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCCGGACACCTACGATCGACCCGTTGATCCCCACATTAAAACGTTGCGTGCCAAGCGGCGCGCCATCAAcc                                                                                                                                                                                                                           >  2:448636/1‑251 (MQ=255)
                          gTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                                                                                                                                                                                                                                           >  1:202321/1‑232 (MQ=255)
                                    aTCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAGCGCTGACTCG‑GTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGA‑TGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCAACATTAAAACGCTGCGTGCCAAGCTGCGCGCCCTCAACCCCGCTCTTTCAcc                                                                                                                                                                                                              >  2:450372/1‑251 (MQ=255)
                                             ggtcATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGCGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATAc                                                                                                                                                                                                     >  1:288939/4‑251 (MQ=255)
                                                      gaagtgaatgaaCCCGCGGCACAGATCGGCG‑GGTTTGACACGCCATTAACGCTGACGCG‑GTTTGATTTTTTATTGCGGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATcgcg                                                                                                                                                                                             <  1:342061/246‑1 (MQ=255)
                                                            aatgaaCCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTAACCGATTAATACTCATCGCGGCATggg                                                                                                                                                                                      >  1:278632/1‑251 (MQ=255)
                                                                accccGCGGCGCAGATCAGTT‑GGTTTGACACGCCATTAACGCTGCCGCG‑TTATGATTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGCCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGatat                                                                                                                                                                                  <  1:257047/249‑1 (MQ=255)
                                                                              aTCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGTCATGGGGTCTAGCCTGGGGGGCCt                                                                                                                                                                    >  2:178870/1‑251 (MQ=255)
                                                                                 aagttggttggccacGCCATTAACGCGGACGGG‑GTATGAGTTTTTATTGCTGAAGACGTTCCTCAAGCCGCGAGGCCGCGTCTGGCCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGt                                                                                                                                                                  <  1:167476/239‑1 (MQ=255)
                                                                                                            gacgcG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCACGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTCTCGGCATGCGGTTgctgct                                                                                                                                       >  1:72510/5‑251 (MQ=255)
                                                                                                                             ttATTGCTAAAGACTTTCCTCAAGTCGCCGGGCCGGGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACt                                                                                                                        <  1:329550/250‑1 (MQ=255)
                                                                                                                               aTTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACtggt                                                                                                                     >  2:465967/1‑251 (MQ=255)
                                                                                                                                  gctgAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTAACGGCATGCGGTTCCTGCTGGGCTATTTTTTACtggtggg                                                                                                                  >  2:103341/3‑250 (MQ=255)
                                                                                                                                      aaGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACAGTTGCTACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGAAATGCGTATCGGCATGCGGTTGCTGCTGGGCAATTTTTTACTggtggaggtg                                                                                                              >  1:321551/1‑251 (MQ=255)
                                                                                                                                        gACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGAGTTTTTTTACTggtggcggtggg                                                                                                            >  2:89134/1‑250 (MQ=255)
                                                                                                                                           gTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGatat                                                                                                                                                                                  >  2:164229/1‑178 (MQ=255)
                                                                                                                                           gTTACTAAAGTCGCCAGGCGGCGTCGGGCCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTGTCACCGATTAATACTCATCGCGGCATGGGatat                                                                                                                                                                                  <  1:164229/178‑1 (MQ=255)
                                                                                                                                                 cAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTgg                                                                                                     <  1:15205/249‑1 (MQ=255)
                                                                                                                                                             ccGCGTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGctat                                                                                       <  1:11497/251‑4 (MQ=255)
                                                                                                                                                               gcgTCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAACCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTAttt                                                                                     >  1:303880/1‑251 (MQ=255)
                                                                                                                                                                 gTCTGGTCCCGCCAGCAA‑TGAATGGATAGCGTGTGGGAAGTTGCGCAGGACGCCTACGATCGCACCGTTGAAACCCACATTAAAACGTTGCGTGGCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATGTTTTACTGGTGGCGGTTGCGGCCTGGTTCGTACTGGCTAttttt                                                                                   <  1:311463/251‑1 (MQ=255)
                                                                                                                                                                  tCTGGTCCCGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGGGGCGGCCTGGTTCGTACTGGctaattttt                                                                                  >  1:99186/1‑243 (MQ=255)
                                                                                                                                                                    tGGTCCCGCCAGCAACTAA‑TGGTAAGCGTGTGGGACGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCCTGGGATATAGCCTGAGGGGACTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTAATTGTGGCGGTGgccgccagtcttggacaggccttttttgct                                                                                >  2:124615/1‑227 (MQ=255)
                                                                                                                                                                     ggccccGCCAGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGt                                                                                 <  1:309595/246‑1 (MQ=255)
                                                                                                                                                                             aaGCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAACCGTTGCGTGCCAAGCTGCGTGCCATCAACCCTGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGGATCGGCATGCGGTTTCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGTTGCGTACTGGCTATTTTTGTCAAAGAAGtt                                                                       <  1:409656/250‑1 (MQ=255)
                                                                                                                                                                               gCAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTa                                                                      >  1:397315/1‑250 (MQ=255)
                                                                                                                                                                                cAACTAA‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGACTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTGGCTATTTTTGTCAAAGAAGTTaa                                                                     >  1:485483/1‑250 (MQ=255)
                                                                                                                                                                                  aCTGA‑TGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCTGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAAcc                                                                  <  1:450372/251‑1 (MQ=255)
                                                                                                                                                                                   ctaa‑TGGCTAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTTGTTCTGGCTATTTTTGTTAAAGAAGTTAAACCg                                                                 >  2:456642/4‑251 (MQ=255)
                                                                                                                                                                                     aa‑TGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCg                                                                 <  1:456642/248‑1 (MQ=255)
                                                                                                                                                                                                    gggAAGCTGCGCATGACACCTCCGATCGCACCGTTGATACCCAAATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTCTTTCACCGATTACTACTAATCGCGGCATGGGCTATAGCCGGAGGGGCCTGTAATGCGTATCGGCATGCGGGTGCTGCTGGGCTATTTTTTACTGGTGGAGGTGGCGGCCTGGTTCGGACTGGCACTTTTTTTCAaagaaggtcacacggggcggcgacgggccg                                                 >  2:432334/1‑227 (MQ=255)
                                                                                                                                                                                                     ggAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAAc                                                <  1:150218/251‑1 (MQ=255)
                                                                                                                                                                                                      gAAGATGCGCCGGACACCTACGATCGCTCCGTTGATACCCACATTAAAACGTGGCGGGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACg                                               <  2:172694/251‑1 (MQ=255)
                                                                                                                                                                                                         gATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCACAGAAGTTAAACCGGGCGTGCGAAGAGCAACgg                                              >  1:189281/1‑249 (MQ=255)
                                                                                                                                                                                                              gcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTGGCTATTTTTGTCAAAGACGTTAAACCGGGCGTGCGAAGAGCCACGgggggg                                         >  2:500954/1‑249 (MQ=255)
                                                                                                                                                                                                               cAGGACACCTACGAGCGCACCGTGGACACCCACATTAAAACGTTGCGTTCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCTCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTCTTACTGGTGGCGGTGGCGGCCTGGTTCGTATTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACg                                      <  1:170/251‑1 (MQ=255)
                                                                                                                                                                                                                    cacatACGATCGCACCGTTGATACCCACATTAACCCGTTGCGTACCAAGCTGCGCGCCATCAACCCCGATCTTCCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGGACTGGCCATTTTTGTTAACGTTGTTAAAACGGGCGGGCGGAGAGCAACG‑Aggggggagttta                                 >  2:326001/5‑242 (MQ=255)
                                                                                                                                                                                                                     acCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCATGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGGGGAGGTGACGGCCTGGTTAGTAATGGCTAATTTTGTCAAAGCAGGTAAATCGGGCGTGCGAATAGCAAAGGTGGGGACGttaat                                 >  1:300371/1‑247 (MQ=255)
                                                                                                                                                                                                                               gCACCGTAGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTGGCTATTTTTGTCACAGAAGGTAAACCGGGGGTGCGAAGGGCAAAGGGGGGGGCGttactagataccacc                       >  2:358622/1‑237 (MQ=255)
                                                                                                                                                                                                                                   cgttgATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGtt                  <  1:494703/247‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaaaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGgcg  <  1:326222/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
CGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAGCGCTGACTCG‑GTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGA‑TGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCG  >  NC_000913/4636364‑4636854

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