Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,077,912 G→A 6.6% P324L (CCG→CTG)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,077,9120GA 6.6% 190.9 / 2.6 60P324L (CCG→CTG) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/23);  new base A (2/2);  total (35/25)
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.18e-01

TCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCC  >  NC_000913/1077768‑1078058
                                                                                                                                                |                                                                                                                                                  
tCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGcc                                                                                                                                                 >  6:24585/1‑148 (MQ=255)
tCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGcc                                                                                                                                                 >  8:214967/1‑148 (MQ=255)
  ggCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGccc                                                                                                                                                <  2:242372/147‑1 (MQ=255)
        aaTGTGGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGGGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAg                                                                                                                                         <  6:119169/148‑1 (MQ=255)
          tGTTGATACCAATATCTTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTGTTCCATTACGCTGTCATACTGTGCTCGGCTATAACTCGGATGGAGCTCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAg                                                                                                                                         <  1:226615/146‑1 (MQ=255)
           gTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATg                                                                                                                                      >  5:115447/1‑148 (MQ=255)
            ttGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGc                                                                                                                                     <  6:278052/148‑1 (MQ=255)
            ttGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCAGCCCTTCATAGATGc                                                                                                                                     >  3:28768/1‑148 (MQ=255)
            ttGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCAGCCCTTCATAGATGc                                                                                                                                     >  7:122868/1‑148 (MQ=255)
             tGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGc                                                                                                                                     >  8:170066/1‑147 (MQ=255)
                     atatCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGt                                                                                                                            >  2:6306/1‑148 (MQ=255)
                       atCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGt                                                                                                                            >  8:332194/1‑146 (MQ=255)
                            aCTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGc                                                                                                                     >  8:93419/1‑148 (MQ=255)
                               gACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCtt                                                                                                                   >  1:242247/1‑147 (MQ=255)
                                aCGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTa                                                                                                                 <  3:141188/148‑1 (MQ=255)
                                         cagGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCg                                                                                                        <  1:233434/148‑1 (MQ=255)
                                            ggTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGa                                                                                                     <  7:332194/148‑1 (MQ=255)
                                                agtgaTTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCg                                                                                                 >  6:26499/1‑148 (MQ=255)
                                                agtgaTTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCg                                                                                                 >  2:77805/1‑148 (MQ=255)
                                                  tgaTTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGcc                                                                                               >  5:252032/1‑148 (MQ=255)
                                                    aTTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCctg                                                                                             <  5:26499/148‑1 (MQ=255)
                                                    aTTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCctg                                                                                             <  1:77805/148‑1 (MQ=255)
                                                    aTTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCct                                                                                              <  5:24585/147‑1 (MQ=255)
                                                          gACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATa                                                                                       <  7:170066/148‑1 (MQ=255)
                                                           aCGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGTCCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATat                                                                                      >  7:369320/1‑147 (MQ=255)
                                                                 gTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAg                                                                                                                                         >  8:282637/1‑91 (MQ=255)
                                                                 gTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAg                                                                                                                                         <  7:101515/91‑1 (MQ=255)
                                                                 gTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAg                                                                                                                                         >  8:101515/1‑91 (MQ=255)
                                                                 gTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAg                                                                                                                                         <  7:282637/91‑1 (MQ=255)
                                                                    aaGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCat                                                                              <  2:242247/147‑1 (MQ=255)
                                                                        tctTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGCTTGAAATGCGCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGGGAATCGCCTGCTGATAGGAAACCATATACg                                                                         <  6:11269/148‑1 (MQ=255)
                                                                          ttCCATTACCCGGGCATACTGGGCGCGGCTATAACGCGGCTGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCGTTCCCGATGGCGTGAATCGCCTGCTGCTAGGAAACCATATACGcc                                                                       >  1:337748/1‑148 (MQ=255)
                                                                              aTTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCtgtg                                                                   >  6:283956/1‑148 (MQ=255)
                                                                                   ccGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGCTGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATACGAAACCATATACGCCTGTGCATCt                                                              >  5:178921/1‑148 (MQ=255)
                                                                                    cGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGAATCTg                                                             >  6:279891/1‑148 (MQ=255)
                                                                                    cGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGGGCATCTg                                                             >  3:183775/1‑148 (MQ=255)
                                                                                     ggTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTgc                                                            >  6:302042/1‑148 (MQ=255)
                                                                                      gTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATgg                                                                                                          >  2:176409/1‑101 (MQ=255)
                                                                                      gTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATgg                                                                                                          <  1:176409/101‑1 (MQ=255)
                                                                                           aCTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGt                                                      >  5:72652/1‑148 (MQ=255)
                                                                                            cTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGt                                                      >  2:297986/1‑147 (MQ=255)
                                                                                              gggCGCGGCTATAACGCGGATGCAGCGCCGAAAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTgcggcgtgccg                                                   >  6:66484/1‑143 (MQ=255)
                                                                                               ggCGCGGCTATAACGCGGATGCAGCGCCTACAGTTTGATTGAAATGCCCGGACCTTCATTGATGCCAAGACCGTTAGACGCTTTACCGATGGGGTGAATAGCCTGCTGATAGGAAACAATATACGCCTGTGCATCTGCGGCGGTAAgc                                                  >  6:157361/1‑148 (MQ=255)
                                                                                                gcgcGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAgcg                                                 <  1:301840/148‑1 (MQ=255)
                                                                                                       tataACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCAGCGGCGGTCAGCGCGGCTTc                                          >  5:159598/1‑148 (MQ=255)
                                                                                                           aCGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCAGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGcc                                       <  8:122868/147‑1 (MQ=255)
                                                                                                           aCGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCAGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGcc                                       <  4:28768/147‑1 (MQ=255)
                                                                                                                gATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCat                                 <  1:6306/148‑1 (MQ=255)
                                                                                                                gATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCat                                 <  1:297986/148‑1 (MQ=255)
                                                                                                                  tGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCatat                               >  2:270819/1‑148 (MQ=255)
                                                                                                                    cAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCg                             >  5:64321/1‑148 (MQ=255)
                                                                                                                        gcCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAag                         <  6:159598/148‑1 (MQ=255)
                                                                                                                           gACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAagag                       >  1:181121/1‑147 (MQ=255)
                                                                                                                               gTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGATTGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACg                  >  5:247171/1‑148 (MQ=255)
                                                                                                                                tttGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACg                  >  3:136876/1‑147 (MQ=255)
                                                                                                                                 ttcattgaAATTCCCGTAAATTCATAGATGACACGAACTTTAGACGCTTTACCGATGCAGTGAATCGCCTGCTGATAGGAAACCATAAACGCATTTGCATCTTCTGCGGTCAGAGATGCTTCGCCAAGCATATCGTAAGAGTAACgt                 >  8:51822/4‑146 (MQ=255)
                                                                                                                                    attgaAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGaaa              <  5:261589/147‑1 (MQ=255)
                                                                                                                                        aaaTGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCtt          <  7:214967/147‑1 (MQ=255)
                                                                                                                                          aTGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTtct       >  2:49487/1‑148 (MQ=255)
                                                                                                                                               cGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTcc  <  7:93419/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCC  >  NC_000913/1077768‑1078058

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