Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,726,615 G→A 12.7% A256V (GCG→GTG)  brnQ ← branched‑chain amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,726,6150GA12.7% 133.1 / 8.1 55A256V (GCG→GTG) brnQbranched‑chain amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base A (4/3);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

ACTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAAT  >  NC_002947/4726472‑4726760
                                                                                                                                               |                                                                                                                                                 
aCTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATg                                                                                                                                              >  3:230042/1‑149 (MQ=255)
  tCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCACGacatc                                                                                                                                            >  1:144740/1‑145 (MQ=255)
    ggCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGATGGTGGCGTCAGCCGCGATGTCGt                                                                                                                                          <  7:47303/149‑1 (MQ=255)
          ggCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGCCACCAGCCACGACATCGCGGCtga                                                                                                                                    >  1:135649/1‑148 (MQ=255)
           gCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCACGACATCGCGGCtgac                                                                                                                                   <  2:144740/149‑3 (MQ=255)
             acaGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCtggtgg                                                                                                                                 >  2:173884/1‑149 (MQ=255)
             acaGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGAGCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCTCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCtggtgg                                                                                                                                 <  8:66835/149‑1 (MQ=255)
                ggTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGCCACCAGCCACGACATCGCGGCtgacgccac                                                                                                                              <  2:135649/149‑8 (MQ=38)
                 gTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCAcc                                                                                                                             <  6:200928/149‑1 (MQ=255)
                    aGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAg                                                                                                                          >  1:135326/1‑149 (MQ=255)
                     ggCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAgg                                                                                                                         <  1:173884/149‑1 (MQ=255)
                      gCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGc                                                                                                                        <  4:230042/149‑1 (MQ=255)
                      gCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGc                                                                                                                        <  3:313796/149‑1 (MQ=255)
                         aaCCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCgg                                                                                                                      >  8:210635/1‑148 (MQ=255)
                         aaCCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGa                                                                                                                     <  8:267240/149‑1 (MQ=255)
                            cgcgGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCACGacatcgc                                                                                                                                          <  4:72778/125‑7 (MQ=255)
                            cgcgGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCACGacatcgc                                                                                                                                          >  3:72778/1‑119 (MQ=255)
                                       ggcaggcCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAg                                                                                                       <  7:231996/149‑1 (MQ=255)
                                              cAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAggg                                                                                                <  7:92839/149‑1 (MQ=255)
                                                gcgcGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGc                                                                                               >  3:163295/1‑148 (MQ=255)
                                                         gCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCCCCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGCACAGGCTGATGTAGACCAGGGCCCGCCCGaca                                                                                     >  5:331456/1‑149 (MQ=255)
                                                         gCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGaca                                                                                     >  4:68333/1‑149 (MQ=255)
                                                          cACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGacac                                                                                    >  3:110606/1‑149 (MQ=255)
                                                               ccAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCCCCCAGGCGGAACAGGCTGATGTAGCCCAGGGCCAGCCCGACACCGGc                                                                                >  7:154044/1‑148 (MQ=255)
                                                                        cGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCg                                                                      >  4:167922/1‑149 (MQ=255)
                                                                             cccAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCgat                                                                 >  3:253213/1‑149 (MQ=255)
                                                                             cccAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCgat                                                                 >  4:276221/1‑149 (MQ=255)
                                                                                 gggAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGatgatg                                                             <  5:11292/149‑1 (MQ=255)
                                                                                   gAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGc                                                           <  4:163295/149‑1 (MQ=255)
                                                                                       cGAAGGTATGCTGTACGTAGGCATGCAGCACCGCTGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAg                                                       >  2:44585/1‑149 (MQ=255)
                                                                                       cGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCTCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGCCCAGGGCCAGCCCGACACCGGCGATCCGGCCGGCGATGCTGGCGTAg                                                       >  6:144486/1‑149 (MQ=255)
                                                                                        gAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGatgatg                                                             <  7:31007/142‑1 (MQ=255)
                                                                                        gAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGatgatg                                                             >  8:31007/1‑142 (MQ=255)
                                                                                         aaGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCg                                                     <  4:253213/149‑1 (MQ=255)
                                                                                          aGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCgg                                                    <  4:110606/149‑1 (MQ=255)
                                                                                             tATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGt                                                  >  6:199709/1‑148 (MQ=255)
                                                                                              aTGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACACAGGCGGAACATGCTGATGTAGACCCGGGCCAGCCAGACACCGGCGATAAGGCAGGCGATGATGGCGTAGCGGGTa                                                 >  4:249609/1‑148 (MQ=255)
                                                                                               tGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTaa                                                <  3:68333/148‑1 (MQ=255)
                                                                                                 cTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCa                                             <  5:199709/149‑1 (MQ=255)
                                                                                                  tGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAg                                            >  1:47093/1‑149 (MQ=255)
                                                                                                      cGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGCAGCGGGTAATCAGCCTc                                        >  7:225220/1‑149 (MQ=255)
                                                                                                       gTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGc                                                           <  3:295186/129‑1 (MQ=255)
                                                                                                       gTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGc                                                           >  4:295186/1‑129 (MQ=255)
                                                                                                             aTGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACt                                 >  2:169897/1‑149 (MQ=255)
                                                                                                                  gcaCCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGaca                            >  3:89957/1‑149 (MQ=255)
                                                                                                                      cgccgcACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACcgc                        >  2:119131/1‑149 (MQ=255)
                                                                                                                       gccgcACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACcgcg                       <  3:276221/149‑1 (MQ=255)
                                                                                                                          gcACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGCAGCGGGTAATCAGCCTCGGCGACTCGACACCGCgcga                    <  8:225220/149‑1 (MQ=255)
                                                                                                                            aCCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCgcgagc                  <  7:210635/149‑1 (MQ=255)
                                                                                                                                    tggCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCa          <  2:47093/149‑1 (MQ=255)
                                                                                                                                        gTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGa      <  2:135326/149‑1 (MQ=255)
                                                                                                                                        accAGCCACGACATCGCGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGa      >  5:185600/3‑149 (MQ=255)
                                                                                                                                             ccGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACaa   >  6:288658/1‑147 (MQ=255)
                                                                                                                                             ccGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACaa   <  5:288658/147‑1 (MQ=255)
                                                                                                                                             ccGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCGCGGCGACTCGGCACCGCGCGAGCGAATGGCATTGACAAt  >  5:286756/1‑148 (MQ=255)
                                                                                                                                              cgcgATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCAtt        <  7:56135/141‑1 (MQ=255)
                                                                                                                                              cgcgATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCAtt        >  8:56135/1‑141 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
ACTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAAT  >  NC_002947/4726472‑4726760

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: