Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,014,746 T→C E104E (GAA→GAG yqeC ← putative selenium‑dependent hydroxylase accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,014,7460TC98.2% 197.3 / ‑4.2 57E104E (GAA→GAGyqeCputative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (33/23);  total (34/23)
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.00e-01

CATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCTTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATAT  >  NC_000913/3014605‑3014891
                                                                                                                                             |                                                                                                                                                 
cATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGt                                                                                                                                           >  2:300395/1‑150 (MQ=255)
cATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGt                                                                                                                                           >  1:591329/1‑150 (MQ=255)
   cacaGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTaaa                                                                                                                                        <  2:350941/150‑1 (MQ=255)
     caGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATcc                                                                                                                                     >  1:157910/1‑151 (MQ=255)
      aGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATccc                                                                                                                                    >  1:338963/1‑151 (MQ=255)
          aTCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGc                                                                                                                                 >  1:108641/1‑150 (MQ=255)
                  gcaACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTccc                                                                                                                        <  2:525416/151‑1 (MQ=255)
                   caACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCt                                                                                                                       <  2:34874/151‑1 (MQ=255)
                    aaCTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCTTCTGGCGTAAATCCCTGCACTTTTccc                                                                                                                        >  2:300545/1‑149 (MQ=255)
                      cTGCTTTTGGGTATGCAAGGTTCGTGATCATCAGGCGTTTTTAACGGCATTCCACGCGAGACCTCTGCTTCAATGAGAATACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGtt                                                                                                                   >  1:576321/1‑151 (MQ=255)
                        gCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTc                                                                                                                  <  2:274346/151‑1 (MQ=255)
                        gCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTc                                                                                                                  >  1:698502/1‑151 (MQ=255)
                            ttGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCtt                                                                                                              >  1:85435/1‑151 (MQ=255)
                             tGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCtt                                                                                                              <  2:23566/150‑1 (MQ=255)
                                gTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCa                                                                                                          <  2:22641/151‑1 (MQ=255)
                                 tATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTcc                                                                                                           >  1:787649/1‑149 (MQ=255)
                                    gCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCt                                                                                                       <  2:209023/150‑1 (MQ=255)
                                           tCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAAc                                                                                               >  1:330646/1‑151 (MQ=255)
                                            cGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAAc                                                                                               <  1:24187/150‑1 (MQ=255)
                                            cGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAAc                                                                                               <  2:483871/150‑1 (MQ=255)
                                                   atcaGGCGCTTTCAACGGCATTCCACCCCATCCATCTGCTTCAATGAGAATTACGTCACTCTCTGGGCGTTGCACCAGTGCATCAATCGCCTCTGTCGCAAATCCCTGCCCTTTTCCCTGGTTCGCTTTCTAGCTGTGAAAACAaactg                                                                                         <  1:33686/149‑1 (MQ=255)
                                                      aGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAaaatgaaatg                                                                                    >  2:299969/1‑151 (MQ=255)
                                                      aGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAaaatgaaatg                                                                                    <  2:690503/151‑1 (MQ=255)
                                                        gcgcTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATggg                                                                                  >  1:200425/1‑151 (MQ=255)
                                                         cgcTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATgggg                                                                                 >  1:584689/1‑151 (MQ=255)
                                                            ttttAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGAt                                                                              >  2:236746/1‑151 (MQ=255)
                                                              ttAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGt                                                                            >  1:808914/1‑151 (MQ=255)
                                                                aaCGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTaa                                                                          <  2:689409/151‑1 (MQ=255)
                                                                 aCGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAag                                                                         <  2:76670/151‑1 (MQ=255)
                                                                    gCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAag                                                                         >  1:704696/1‑148 (MQ=255)
                                                                     cATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGAc                                                                     >  2:774271/1‑151 (MQ=255)
                                                                          cACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATg                                                                <  2:603361/151‑1 (MQ=255)
                                                                            cgcgAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAg                                                              <  2:113881/151‑1 (MQ=255)
                                                                                   cATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCAt                                                       <  2:503049/151‑1 (MQ=255)
                                                                                    aTCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATg                                                      >  1:374580/1‑151 (MQ=255)
                                                                                     tCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATg                                                      >  2:522778/1‑150 (MQ=255)
                                                                                      cTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGc                                                    >  2:111564/1‑151 (MQ=255)
                                                                                          ttCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCggg                                                 >  1:320262/1‑150 (MQ=255)
                                                                                               tGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACg                                                                    >  1:186039/1‑126 (MQ=255)
                                                                                               tGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACg                                                                    <  2:186062/126‑1 (MQ=255)
                                                                                                 agaATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGAc                                         >  2:370033/1‑151 (MQ=255)
                                                                                                   aaTTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAg                                       <  1:72390/151‑1 (MQ=255)
                                                                                                        cGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAAc                                  >  1:558114/1‑151 (MQ=255)
                                                                                                        cGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAAc                                  >  2:470205/1‑151 (MQ=255)
                                                                                                         gTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAAcc                                 <  1:810683/151‑1 (MQ=255)
                                                                                                           cacaCGCTGGTCGTGGCACCAGTGCATCAATCGCCTCTGGCATAAACCCCTGCCCGTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAAATAAATGGGGGAAGTAGGAGGCGCATGAAGAAGCATGGCGGGATCCCGGCAGa                                      <  1:430904/144‑1 (MQ=255)
                                                                                                              acTCTGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACggg                            >  1:320658/1‑151 (MQ=255)
                                                                                                                  tGGTCGTTGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGcc                          >  2:549533/1‑149 (MQ=255)
                                                                                                                         tgtaCCAGTGCATCAATCGTCTATGGCGTAAATTCCTGCACTTTTCGCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGTAATGGGGGATGTAAGAGACGCATGAGGAAGTATAGGGGGATCACTGAAGAAAACCACGGGGCAgaga                     >  1:360449/4‑147 (MQ=255)
                                                                                                                         tGCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACg                 >  2:496919/1‑151 (MQ=255)
                                                                                                                          gCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACg                 <  2:585972/150‑1 (MQ=255)
                                                                                                                          gCACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGt                >  2:321949/1‑151 (MQ=255)
                                                                                                                           cACCAGTGCATCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGtt               >  2:632422/1‑151 (MQ=255)
                                                                                                                                    aTCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATaaa       >  1:511861/1‑150 (MQ=255)
                                                                                                                                    aTCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATaaa       >  2:160261/1‑150 (MQ=255)
                                                                                                                                    aTCAATCGCCTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAAc      <  2:422130/151‑1 (MQ=255)
                                                                                                                                        aTCGCCTCTGGCTTAAATCCCTGCACTTTTCCCTTGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGTGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACatgt  <  1:856231/151‑3 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAATTACGTCACACTCTGGTCGTTGCACCAGTGCATCAATCGCTTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATAT  >  NC_000913/3014605‑3014891

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