Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 588 G→A G84G (GGG→GGA thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135880GA97.2% 102.6 / ‑5.7 36G84G (GGG→GGAthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/17);  minor base C (0/1);  total (18/18)
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.44e-01

CAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGAT  >  NC_000913/365‑834
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
aaggggCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATCCGGGAAAGCAATGCCAGGCAGGGGCGGGTGGCCCCCGTCCTCCCTGCCCCCGCCAAAAGCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACttgc                                                                                                                                                                                                                             <  2:289801/247‑3 (MQ=255)
    ggCAAATTCAAAACTTTTTCTGCGGGTTCCCGATATTCTGGAAAGCATTGCCAGGCAGGGGCAGTTGGCCACCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGTATTCCCGCTGGCGCTATTGAAAACTTtcgttg                                                                                                                                                                                                                         <  2:337889/251‑3 (MQ=255)
     gCAAAGGCAGAACGTTTCCTGCGGGTTGCCGATATGCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGa                                                                                                                                                                                                                        <  1:222143/251‑1 (MQ=255)
          tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGAAGCCAAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCACCCAGCCGGGATTCCCGCTGGCGCCAtt                                                                                                                                                                                                                                       >  1:74543/1‑231 (MQ=255)
             aGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAAt                                                                                                                                                                                                                >  1:97414/1‑251 (MQ=255)
                           gggtTGCCGATATTCTGGAAAGCAATGCCGGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATaaaa                                                                                                                                                                                                  <  1:454427/249‑1 (MQ=255)
                                 ccGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCCGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTc                                                                                                                                                                                            >  1:67925/1‑251 (MQ=255)
                                      atTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGc                                                                                                                                                                                        >  1:135177/1‑250 (MQ=255)
                                        tCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCAt                                                                                                                                                                                      <  1:217485/250‑1 (MQ=255)
                                                    ccccAGGCAGGGGCTGGTGGCCTCCGCCCTACCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACCCGCCGCCGCCCAGCCGTGATTCCCGCTGGCGCAATTGAAAACTGTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTttgt                                                                                                                                                                         <  2:92312/249‑1 (MQ=255)
                                                      ccAGGCAGGGGCAGGGGGCCACCGTCCACTCTGCCCCCGCCCAAATCACCAACCACCTGGTGGAGATGATTGAAAAAACCATTCGCGGCCAGGATGGTTTAACCAATATCAGCGGAGCCGAACGTATTTTTGCCGAACTTTTGACGGGGCTCGCCGCCGCCCAGCCGGGAATCACGTTGGCGCAATTGACAACTTTCGTCGACCAGGAATTTGCCCAAATCAAACAGGACCTGAATGGCCTTAGTttgttg                                                                                                                                                                       >  2:475045/1‑251 (MQ=255)
                                                        gggcaggggaaggTGCCCACCGTCCTCCCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCACGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTggg                                                                                                                                                                     <  2:141340/250‑1 (MQ=255)
                                                             ggtgcaggTGGCCACCGTCTCCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGt                                                                                                                                                                <  2:328225/248‑1 (MQ=255)
                                                                       gCCACCATCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGCAAAAACCATTAGCGGCCAGGATGCTTTACCACATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGCATCGCCGCCGCACAGCCGGGATTCCCGCCGGCGCAATTGAAAACTTCCGTCGATCCGGAATTTGCCCCAATAAAACATGTCCCGCATGGCCTTAGTTTGGTGGGGAAGTGACCgacg                                                                                                                                                        >  2:85033/1‑246 (MQ=255)
                                                                                                 aaTCACCAACCACCTGGTGGCGATGATTGAAAAAACCATAAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCCGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGCTCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACCCTTCGTTGATTTGCCGTg                                                                                                                            >  2:181346/1‑251 (MQ=255)
                                                                                                                          atTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTAt                                                                                                   >  1:136239/1‑251 (MQ=255)
                                                                                                                                     cATTAACGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCggtt                                                                                        >  2:391195/1‑248 (MQ=255)
                                                                                                                                           cGGCCAGGATGCTTTACCCAATAGCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTAGTAGAAg                                                                                  <  2:67925/251‑1 (MQ=255)
                                                                                                                                               cAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGcc                                                                                                                               >  1:313887/1‑202 (MQ=255)
                                                                                                                                               cAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGcc                                                                                                                               <  2:313887/202‑1 (MQ=255)
                                                                                                                                                                tatCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTAccg                                                             >  1:192687/1‑249 (MQ=255)
                                                                                                                                                                        aTGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGAt                                                                                                                                                        >  2:235992/1‑152 (MQ=255)
                                                                                                                                                                        aTGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGAt                                                                                                                                                        <  1:235992/152‑1 (MQ=255)
                                                                                                                                                                          gcccaaCGTATTTTTTCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATcc                                                   <  2:136239/247‑1 (MQ=255)
                                                                                                                                                                           ccGAACGTATTTTTGCCGACCTTTGGACGGGACTGGCCGCCTCCCAGCCGTGCGTCCAGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGGCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTTGCGAGAAAAGGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCg                                                  <  2:384585/251‑1 (MQ=255)
                                                                                                                                                                             gAACGTCTTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCCATTGAAAACTTTCGTCGATCAGGAATTTACCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCCTTATGGCCGGCGTATTAGAAGCGCGCGGTCACACCGTTACCGTTATCGATCCggg                                                >  2:433532/1‑250 (MQ=255)
                                                                                                                                                                               aCGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCg                                              <  1:310703/251‑1 (MQ=255)
                                                                                                                                                                                        tGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAActg                                      <  1:360430/250‑1 (MQ=255)
                                                                                                                                                                                           cGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAACCATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTgg                                  >  1:215785/1‑251 (MQ=255)
                                                                                                                                                                                             aaCTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTCTTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGc                                 >  1:131156/1‑250 (MQ=255)
                                                                                                                                                                                                          aCTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTAcc                   <  1:348599/251‑1 (MQ=255)
                                                                                                                                                                                                               ccgccgccCAGCCGGAATTCCCGCTGGCGCAATTGCAAACTTTCGTCGATCAGGACTTTGCCCAAATCAACAATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATCGCATCAACGCTGCGCTGCTTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGACGCGAGCGGGCACAACGTTACCGTTATCGATCCGGTCGAAAAACTTCTGGCAGTGGGGCATTACCTCGa               >  2:66342/1‑250 (MQ=255)
                                                                                                                                                                                                                cgcagccCCGCCGGGATTCCCGCCGGCCCAATTGAAAACTTTCGTCGATCAGGAAGTTGCACAAATAAAACATGTCCTGCATGGCAGTAGTTTGTGGGGGCAGGGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGAGCCGGTCGAAAAACTGCTGGCAGTGGGGCATTCCCTCGAAt             <  2:327292/247‑1 (MQ=255)
                                                                                                                                                                                                                 gccgccCATCCGGGATTCCCGCTGGCGCAAGTGAAAACTTTCGGCGATCAGGAATTTGCCCAAATAAAACAGGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGATAAAATGTCGATCGCCATTATGGCCGGCGTATTATAAGCGCGCGGTCACAACGTTACCATTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATc            <  1:332090/251‑1 (MQ=255)
                                                                                                                                                                                                                         gCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGGTCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTTCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACACAGTTGCCGTTATCGATCCGGTCGAAAAAATGCTGGCAGAGAGGCATTACCTCGAATCTACCGTcc    >  2:336790/1‑250 (MQ=255)
                                                                                                                                                                                                                           cGGGATTCCCGATGGCGCAATTGAAAACTTTCGTCGCTCAGGAATTTGACCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAATCGCGCGGTCACAACGTTACCGTTATCGCTCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCggt  >  2:88003/1‑249 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
CAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGAT  >  NC_000913/365‑834

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