Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,726,619 T→C 7.1% I255V (ATC→GTC)  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,6190TC 7.1% 155.9 / 2.7 56I255V (ATC→GTC) brnQbranched‑chain amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base C (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

ACTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGATA  >  NC_002947/4726472‑4726767
                                                                                                                                                   |                                                                                                                                                    
aCTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATg                                                                                                                                                     >  3:230042/1‑149 (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)
                                       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)
                                                                                                                                                 gATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGa    >  7:134277/1‑149 (MQ=255)
                                                                                                                                                  aTGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGAt   <  4:89957/149‑1 (MQ=255)
                                                                                                                                                   tGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGATa  >  5:26446/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGATA  >  NC_002947/4726472‑4726767

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: ‑