Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,373,506 A→T 30.0% N90K (AAT→AAA pbuG ← NCS2 family nucleobase:cation symporter‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,373,5060AT30.0% 81.5 / 48.0 63N90K (AAT→AAApbuGNCS2 family nucleobase:cation symporter‑2
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base T (11/8);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

GCATTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAAC  >  CP000730/2373375‑2373646
                                                                                                                                   |                                                                                                                                             
gCATTAATAATAACTTCCCTAAAAAGTCTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATccc                                                                                                                                      >  2:760554/1‑141 (MQ=255)
       aaataaCTGCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTg                                                                                                                               <  1:387546/140‑1 (MQ=255)
         ataaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATATATCCACTTGTTAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGt                                                                                                                             <  1:573626/141‑1 (MQ=255)
         ataaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGt                                                                                                                             <  2:604687/141‑1 (MQ=255)
         ataaCCTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGt                                                                                                                             <  1:313598/141‑1 (MQ=255)
          taaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTg                                                                                                                            >  2:216651/1‑141 (MQ=255)
           aaCTTCCCTAAACCATGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGc                                                                                                                           >  2:536484/1‑141 (MQ=255)
                  cTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGACTGCTTTTAATCCCATACCTGGTGCTAACGCa                                                                                                                   >  1:264708/1‑141 (MQ=255)
                    aaaCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAAt                                                                                                                 >  2:749874/1‑141 (MQ=255)
                     aaCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCCACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAAtt                                                                                                                >  1:139207/1‑141 (MQ=255)
                      aCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTg                                                                                                               <  1:312737/141‑1 (MQ=255)
                      aCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTg                                                                                                               >  2:15338/1‑141 (MQ=255)
                       cccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTCCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTgc                                                                                                              <  2:694480/141‑2 (MQ=255)
                        ccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGa                                                                                                             <  1:9846/141‑1 (MQ=255)
                        ccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATACTTTTAATCCCATACCTGGTGCTAACGCAATTGGa                                                                                                             <  2:217618/141‑1 (MQ=255)
                             gACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGagcttt                                                                                                        >  1:487912/1‑136 (MQ=255)
                              aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTa                                                                                                       >  2:670072/1‑141 (MQ=255)
                              aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAAATTGCCAAGGAATACCCATGGTTAACAAAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTggcgacgga                                                                                                       >  2:118795/1‑134 (MQ=255)
                               cGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAg                                                                                                      >  2:349437/1‑141 (MQ=255)
                                 gTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCt                                                                                                    <  1:544438/141‑1 (MQ=255)
                                    aaTATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCGAACGAAATTGGATATTTAGCTAtt                                                                                                 >  2:750108/1‑141 (MQ=255)
                                      tatCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAg                                                                                               <  1:36666/141‑1 (MQ=255)
                                      tatCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAg                                                                                               <  2:785624/141‑1 (MQ=255)
                                             aaGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATg                                                                                        <  1:642526/141‑1 (MQ=255)
                                               gAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGaa                                                                                      >  1:277289/1‑141 (MQ=255)
                                                    taCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTGGCCCCATGAATAGCg                                                                                 >  1:123969/1‑141 (MQ=255)
                                                    taCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCg                                                                                 >  1:123915/1‑141 (MQ=255)
                                                        tGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGcc                                                                             <  1:503784/141‑1 (MQ=255)
                                                        tGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGcc                                                                             <  1:443516/141‑1 (MQ=255)
                                                        tGACAATACAACGGCTTGTAAACCAACTTGCCAAGGAATACCTCTTGTTACCACAACAGTAAATGCAAAGAATGCTTTTAATCCCATTCCTGTTGCTAACTCAATTGGATATTTAGCTATTATCCCCATGAATAGCGAGcc                                                                             <  2:141350/141‑1 (MQ=255)
                                                         gagaATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCt                                                                            <  1:270242/141‑1 (MQ=255)
                                                          agaATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTa                                                                           <  2:282114/141‑1 (MQ=255)
                                                           gaATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTAc                                                                          >  2:188293/1‑141 (MQ=255)
                                                             aTAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACaa                                                                        <  2:85842/141‑1 (MQ=255)
                                                                      cTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTa                                                               <  1:135593/141‑1 (MQ=255)
                                                                       tGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTaa                                                              <  1:475659/141‑1 (MQ=255)
                                                                        gTTATACCAACTTGCAAAGGAATACCCATGGTTAACACCACAGTAACTGCAAAGAATGCATTCAATCCCATACCTGATGCCAACGCAATTGGATGTTTAGCTATTAGCCCCATGAATAGCGAGCCTCCAAATGCTGCCAAt                                                             >  2:220159/1‑141 (MQ=255)
                                                                                  cTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGACTACAAGTTCTGCTAATGCAGTCGCTa                                                   >  1:212707/1‑141 (MQ=255)
                                                                                    tGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCAATGAATAGCGAGCCTACAAATGGTGCTAATGCAGTCGCGACa                                                 >  2:551581/1‑141 (MQ=255)
                                                                                     gCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaa                                                >  2:711031/1‑141 (MQ=255)
                                                                                     gCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaa                                                >  1:285384/1‑141 (MQ=255)
                                                                                             attcccATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCa                                        <  1:743743/138‑1 (MQ=255)
                                                                                             aTACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCa                                        <  1:412127/141‑1 (MQ=255)
                                                                                                   aTGGTTAACACAACAGTAAATGTAAAGAAGGCATTTAATCCCTACTTTGTACCAAGCTCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTgg                                  <  1:771188/141‑1 (MQ=255)
                                                                                                   aTGGTTAACACAACAGTAAACGCAAACAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTgg                                  <  1:457185/141‑1 (MQ=255)
                                                                                                    tGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGt                                 <  1:760554/141‑1 (MQ=255)
                                                                                                     ggTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATACCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTc                                <  1:223000/141‑1 (MQ=255)
                                                                                                        tAACACAACAGTAAATGCAATGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAt                             <  1:15338/141‑1 (MQ=255)
                                                                                                        tAACACAACAGTAAATGCAAAGAATGCTTTTTATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAt                             >  1:219931/1‑141 (MQ=255)
                                                                                                        tAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAt                             >  2:129882/1‑141 (MQ=255)
                                                                                                          acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                           >  2:11382/1‑141 (MQ=255)
                                                                                                          acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCCACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                           >  2:671824/1‑141 (MQ=255)
                                                                                                              aacaGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCAGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCat                       >  2:404544/1‑141 (MQ=255)
                                                                                                                   tAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtctt                  <  1:418058/141‑1 (MQ=255)
                                                                                                                   aaaaTGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtctt                  <  2:568526/140‑1 (MQ=255)
                                                                                                                    aaaTGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCACCATGAATAGCGAGCCTACAAACGCTGCGAATGCAGTCGCTACAAAAATGGCACCTTGGTACATTCTCATAtcttc                 >  2:278243/1‑141 (MQ=255)
                                                                                                                    aaaTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttc                 >  1:677173/1‑141 (MQ=255)
                                                                                                                     caTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATGTTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTACTATTTCAGTCGCTACAAACATGGCACCTTGGTCCATTTTCATAtgttct                <  2:479737/140‑1 (MQ=255)
                                                                                                                     aaTGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttct                >  1:629449/1‑141 (MQ=255)
                                                                                                                       tGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTCGCTATTAGCCACATGAATAGCGAGCCTACCAATGCTGCTAATGCAGTCGCTGCAGAAATGGCACCTTGGTCCATTTTCATATATTCTGa              >  1:194294/1‑141 (MQ=255)
                                                                                                                            aGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGc         >  1:517771/1‑141 (MQ=255)
                                                                                                                            aGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGc         <  1:87968/141‑1 (MQ=255)
                                                                                                                                 gCATCTAATCCAAGACCTCGTGCAAACGCAATTAGATATTTAGCTATTAGCCCCAGGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCATCTTGGTCCATTTTCATATCTTCTGATACGCCTTTa    <  1:430065/141‑1 (MQ=255)
                                                                                                                                   aTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTCAAAAAATGGCACCTTGGCCCATTTTAATACCTTCTCATaaccatataac  >  2:312481/1‑131 (MQ=255)
                                                                                                                                   |                                                                                                                                             
GCATTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAAC  >  CP000730/2373375‑2373646

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