Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,354,687 C→A V262F (GTT→TTT)  sapB ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,354,6870CA98.2% 217.1 / ‑3.3 56V262F (GTT→TTT) sapBantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (27/28);  total (28/28)
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.29e-01

GCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAG  >  NC_000913/1354545‑1354832
                                                                                                                                              |                                                                                                                                                 
gCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCgg                                                                                                                                           >  1:589757/1‑151 (MQ=255)
gCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCgg                                                                                                                                           <  1:998184/151‑1 (MQ=255)
 ccATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGt                                                                                                                                          <  1:759124/151‑1 (MQ=255)
  cATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTg                                                                                                                                         <  2:282290/151‑1 (MQ=255)
    tGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGat                                                                                                                                       >  1:587082/1‑151 (MQ=255)
     ggCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGatc                                                                                                                                      <  2:278051/151‑1 (MQ=255)
        aCCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGATAAAAAACATTTCGGTGatcatc                                                                                                                                   <  2:694578/151‑1 (MQ=255)
          ccAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGc                                                                                                                                 >  2:564065/1‑151 (MQ=255)
                 atCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTc                                                                                                                          <  1:539433/151‑1 (MQ=255)
                  tCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAAAATTCCGGTGTTCATCGccgcggt                                                                                                                           >  2:1053737/1‑144 (MQ=255)
                   cAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTc                                                                                                                          >  1:720658/1‑149 (MQ=255)
                   cAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCa                                                                                                                         <  1:1047949/150‑1 (MQ=255)
                           aCGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCAtggtgg                                                                                                                >  1:635189/1‑151 (MQ=255)
                            cGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGa                                                                                                               <  1:80979/151‑1 (MQ=255)
                             gTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGa                                                                                                               >  2:670156/1‑150 (MQ=255)
                                  cAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTg                                                                                                         >  1:965903/1‑151 (MQ=255)
                                  cAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTAAGCATGGTGGAAAACt                                                                                                          >  2:350089/1‑150 (MQ=255)
                                   aataatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGt                                                                                                        <  2:455643/151‑1 (MQ=255)
                                      aatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTaa                                                                                                      <  2:742049/150‑1 (MQ=255)
                                      aatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTaa                                                                                                      <  2:759740/150‑1 (MQ=255)
                                        tGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGcc                                                                                                   >  2:942897/1‑151 (MQ=255)
                                        tGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGcc                                                                                                   >  1:1009402/1‑151 (MQ=255)
                                          aCCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGccc                                                                                                  >  1:267017/1‑150 (MQ=255)
                                           ccAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGccc                                                                                                  >  2:55063/1‑149 (MQ=255)
                                           ccAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGccc                                                                                                  >  1:1077427/1‑149 (MQ=255)
                                           ccAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAg                                                                                                <  2:968817/151‑1 (MQ=255)
                                             aGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGc                                                                                              <  1:353641/151‑1 (MQ=255)
                                              gTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCg                                                                                             <  1:745353/151‑1 (MQ=255)
                                                gAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAg                                                                                           <  2:683078/151‑1 (MQ=255)
                                                    cacacaCCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGaa                                                                                        >  1:287311/1‑150 (MQ=255)
                                                    cacacaCCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAAt                                                                                       <  1:858491/151‑1 (MQ=255)
                                                     acacacCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAAt                                                                                       >  1:1009158/1‑150 (MQ=255)
                                                                              gCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGaa                                                             <  1:1041024/151‑1 (MQ=255)
                                                                              gCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGa                                                              >  1:894474/1‑150 (MQ=255)
                                                                                tGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACg                                                           >  1:675256/1‑151 (MQ=255)
                                                                                 gCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAAc                                                            <  2:402888/149‑1 (MQ=255)
                                                                                   gTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTgg                                                        >  2:737704/1‑151 (MQ=255)
                                                                                     aaTCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCg                                                      <  2:506268/151‑1 (MQ=255)
                                                                                      aTCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCg                                                      >  1:1014066/1‑150 (MQ=255)
                                                                                      aTCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGa                                                     <  2:216924/151‑1 (MQ=255)
                                                                                        ccTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGAc                                                    <  2:814715/150‑1 (MQ=255)
                                                                                           gctgGCGGATGGCGGTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAACAAAAATTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGCGGAATAACCGGCGGCAACGCGTTATGCAAAACGTTGCGACGCa                                                 >  2:853563/1‑150 (MQ=255)
                                                                                                 ggATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAg                                           >  1:587752/1‑150 (MQ=255)
                                                                                                     ggCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTaa                                        >  2:776545/1‑149 (MQ=255)
                                                                                                          taattaaCCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGcgcg                                  >  2:544020/1‑150 (MQ=255)
                                                                                                              taaCCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACaa                              <  2:622548/150‑1 (MQ=255)
                                                                                                                   aGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACcgcg                        <  1:742789/151‑1 (MQ=255)
                                                                                                                       cccTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTc                    <  1:668767/151‑1 (MQ=255)
                                                                                                                       cccTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTc                    <  2:504606/151‑1 (MQ=255)
                                                                                                                          tAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgcc                 >  1:581163/1‑151 (MQ=255)
                                                                                                                            aaCCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgccgc               >  2:793909/1‑151 (MQ=255)
                                                                                                                             aCCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgccgcc              >  1:991307/1‑151 (MQ=255)
                                                                                                                             aCCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgccgcc              <  1:1011643/151‑1 (MQ=255)
                                                                                                                                  gCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCttt         <  1:602357/151‑1 (MQ=255)
                                                                                                                                  gCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCttt         <  1:546362/151‑1 (MQ=255)
                                                                                                                                        tAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACAt    >  1:285248/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAg  <  1:523882/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAG  >  NC_000913/1354545‑1354832

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