Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,126,6451.C96.5% 97.9 / ‑3.6 29intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (20/8);  total (20/9)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCAGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTGTGCGGTGCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGTCCTCCACCGCGCGGATTTC  >  AE015451/1126384‑1126940
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        
ccAGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACCCAGAGCCCGCCACCCAATCCAGGTggg                                                                                                                                                                                                                                                                    >  1:186108/1‑300 (MQ=255)
  aGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg                                                                                                                                                                                                                                                                  >  1:561377/1‑301 (MQ=255)
                       ttGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTg                                                                                                                                                                                                                                             <  1:276255/301‑1 (MQ=255)
                                  ggTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGAGAATGCGGAGGTATAGCGACCATCGCATTCCCGGGTAAACCCGCTCCTACACAGCGCCCGCCACCCAAACCAGGTGGCGGGGAGGCCTATCAGACGCCCTGTTAGCTCtcga                                                                                                                                                                                                                                  >  2:270675/1‑301 (MQ=255)
                                                         ggCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                                                                                                                                                                                           >  1:459107/1‑301 (MQ=255)
                                                                     cGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAGGACGCCCAGTGACCTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGATTACCCGCGAACGCGGAGCTATAGCCACCATCGCATTCGCGGTTAAACCCGCCCCTACACAGACCCCGCCACCCAAGCCACCTGGCGGGTCCGGCTATACGGCGTCCTGACAACACTCGATCGACGGGCAGGCCCACACCAGGTTGCCGTCGCCg                                                                                                                                                                                               >  2:709929/1‑301 (MQ=255)
                                                                                    cggcgCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGCCCGACACCAGGTTGCGGTCGCCAGACACGGTGTCGAcc                                                                                                                                                                                >  2:137986/1‑301 (MQ=255)
                                                                                      gcgcggctCGGGGCGACCGTGTGGGGGCTGGAGGGGCCGTAGCGGGTCGTGAACAGGTCGAAGACGCTCAGGGGCATGGGGACGTCCTGGCAGGTTGTTGGGTGGACATTGGGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAGCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                                                                                                                                                                              <  2:208671/294‑1 (MQ=255)
                                                                                      gcgcgcATCGAGCCGTGCGTGTGGGAGCTGGAGGGGGCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCGCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                                                                                                                                                                              <  2:458917/301‑1 (MQ=255)
                                                                                                    gACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACAAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGcc                                                                                                                                                               >  2:608509/1‑301 (MQ=255)
                                                                                                         ggtgGGGGCTGGATGGGCCACTGTCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGATGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc                                                                                                                                                          <  1:16727/300‑1 (MQ=255)
                                                                                                                        ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTCATTGCCTTccacccag                                                                                                                                           >  1:120834/1‑298 (MQ=255)
                                                                                                                            gatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAGCCAAGCCAGCTGGCGGGGATGGCTCTCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGGACGCGTTGTCGACCCGCCCGACCGGCGTGCAACACTTCCGTTCCACCAGGCTcc                                                                                                                                       >  2:351929/1‑300 (MQ=255)
                                                                                                                              tgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTTCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGCGTGCCCGTACTTGCCTTCCACCAGTCTCGcc                                                                                                                                     >  2:313016/1‑301 (MQ=255)
                                                                                                                                         aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATAc                                                                                                                          >  1:113720/1‑301 (MQ=255)
                                                                                                                                          aCAGGTCGAAGACGCTCAGTGCCATGTGGCCCTCCGGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATAcc                                                                                                                         <  1:138045/301‑1 (MQ=255)
                                                                                                                                                                         ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCCGCCAGTACTTGCCTTCCCCCCGGCCCGCCTGGGTGTATACCGCCCGCGCCCGGCTGTAGCCATGGGCCCAtt                                                                                          >  2:405943/1‑301 (MQ=255)
                                                                                                                                                                         ctactGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCAGTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCAtt                                                                                          <  1:319087/298‑1 (MQ=255)
                                                                                                                                                                                           gTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGccc                                                                        >  1:317579/1‑300 (MQ=255)
                                                                                                                                                                                            tGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCGCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGATGTCGACCCGGCCGACCGGAGGCCAGTACTTGCCTGCCACCAGGCTCGCCAGGGGGGATACCGCCTGCTCGCGGCTGTAGCCATGGGGCCATTCGCCAAACAGCtgagaagc                                                                       >  1:80225/1‑294 (MQ=255)
                                                                                                                                                                                                 aTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAgtgt                                                                  >  1:616626/1‑301 (MQ=255)
                                                                                                                                                                                                    gAGCACTGTAGGCGCGGGTTTAGCCGCGAATGGGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTgtgcg                                                               <  1:710229/301‑1 (MQ=255)
                                                                                                                                                                                                            tAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTGTGCGGTGCGttt                                                       <  1:608770/301‑1 (MQ=255)
                                                                                                                                                                                                               gcgcGGGTTTACCCGCGACTGCGGATGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGGGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTGTGCGGTGCGTTTTTc                                                    <  2:179827/301‑1 (MQ=255)
                                                                                                                                                                                                                        tACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTGTGCGGTGCGTTTTTCAGCGGGTTg                                           >  1:572527/1‑301 (MQ=255)
                                                                                                                                                                                                                               gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCCCAGCCGCAGGGTGCGGTGCGTTTTTCAGCGCGTTGTCGGCCt                                    >  2:22095/1‑301 (MQ=255)
                                                                                                                                                                                                                                                 accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTGTGCGGTGCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGtcct                  >  1:65427/1‑301 (MQ=255)
                                                                                                                                                                                                                                                 accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGCGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCTGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTCTACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCCCCAACCAGCCCAGCCGCAGCGTGCGGTGCGCTTCTCAGCGGTTTTTCGTCCCGGTCCAGCCTACCGtccc                  >  2:684491/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTGTGCGGTGCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGTCCTCCACCGCGCGGAtgtt  >  1:28252/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        
CCAGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCTCAGCCGCAGTGTGCGGTGCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGTCCTCCACCGCGCGGATTTC  >  AE015451/1126384‑1126940

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