Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,978,089 A→C 18.4% M288R (ATG→AGG)  yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,978,0890AC18.4% 127.8 / 3.9 58M288R (ATG→AGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (12/34);  new base C (9/2);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAC  >  CP009273/2977941‑2978236
                                                                                                                                                     |                                                                                                                                                   
gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCCCCAATCCCAt                                                                                                                                                    >  1:475451/1‑150 (MQ=255)
 acCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCAt                                                                                                                                                    <  2:792424/149‑1 (MQ=255)
   ccGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCaa                                                                                                                                                 <  1:570818/150‑1 (MQ=255)
   ccGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCaa                                                                                                                                                 <  1:414219/150‑1 (MQ=255)
        gCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCg                                                                                                                                            <  1:279965/150‑1 (MQ=255)
            aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGAt                                                                                                                                         <  1:420942/149‑1 (MQ=255)
            aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCGCAAATCCCATCAAATCCGGATc                                                                                                                                        >  1:182490/1‑150 (MQ=255)
               ccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTac                                                                                                                                     <  1:742396/150‑1 (MQ=255)
                cACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATATAACGGCGCAATTTCCCTCAATGCCGGATCTaca                                                                                                                                    >  1:548814/1‑150 (MQ=255)
                     aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTTCGGAGTAATTCCGAGCCATGCGAGATCTACACCa                                                                                                                                 <  2:43493/148‑1 (MQ=255)
                        cATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCt                                                                                                                            <  1:933839/150‑1 (MQ=255)
                        cATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAATTTCCCTCCATGCCCGGTCTACACCCACCCt                                                                                                                            >  2:322159/1‑150 (MQ=255)
                                     gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATGCCGGAGCTACACCAACCCCGGCAAATGcgcg                                                                                                                >  1:1033983/1‑149 (MQ=255)
                                            ccTGTGCAGGAAGCGCTTCATTAGCGTCGATAAGGTCGACTTCAGCCAACTGCCGGCCTACACGCTCCAGGCAACGGCGGGTCGCGTTCTCCTTAATGATTCGGCGCACACCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCggg                                                                                                        <  2:123962/150‑1 (MQ=255)
                                               gtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCCCGCAATTTCCATCCATGCCGGCTCTACACCAACCCCCGCAAATGCGCGAAAGCGGGc                                                                                                       >  2:125864/1‑148 (MQ=255)
                                                tgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCa                                                                                                     <  1:612791/149‑1 (MQ=255)
                                                    gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACg                                                                                                <  2:789637/150‑1 (MQ=255)
                                                      aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGcgcg                                                                                                                <  1:558885/132‑1 (MQ=255)
                                                       aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGt                                                                                             <  1:792153/150‑1 (MQ=255)
                                                       aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGTGCAAACACCCGCGCAAATATCAATCATATCCGCTCTACACCAACGCTGGCCAATTCCCCAATGCGGGGCAGCACcgg                                                                                              >  2:428389/1‑149 (MQ=255)
                                                        aCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTATCCCCGCCCCATTCCCCTCAACTCCCGACCTACACCAAACCCGGCCAAAGCGCGCAAACGGGCCAGCCCgggg                                                                                             >  2:721788/1‑148 (MQ=255)
                                                                    ccTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTTCCCTCAAAGCCCGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCCGCCCCGGGAAATTCAACGccc                                                                                >  1:238573/1‑148 (MQ=255)
                                                                         aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACAACCCCGCCAATTCCCTCAACTCCGGATCTACACCCACCCCGGCCAAAGCGCGAAAGCGGGCCCGCACCGGGAAAATCAACACTcccgct                                                                           >  2:522675/1‑145 (MQ=255)
                                                                          tAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt                                                                          <  1:287137/150‑1 (MQ=255)
                                                                           aaGATCGACTTCAGCCAGCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATCCCCTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCCCGTGCttt                                                                         >  1:533404/1‑150 (MQ=255)
                                                                                 gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCtt                                                                    <  1:264082/149‑1 (MQ=255)
                                                                                 gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTc                                                                   <  1:199323/150‑1 (MQ=255)
                                                                                    ttCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCCTCAATGCCGGGTCTACACCCACCCCGGCAAAAGCGCGAAAGCGGGCCAGCCCCGGGAAAATCAACGCCCCTGCCTTTGTTTCCCt                                                                >  2:422458/1‑150 (MQ=255)
                                                                                      cAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCATCTGCGCAATTCCCAGCAATGCCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCt                                                                <  2:571992/148‑1 (MQ=255)
                                                                                       aGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcat                                                             <  1:740657/150‑1 (MQ=255)
                                                                                       aGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCATCAATGCCGGATCTACACCCACCCTGGCAAAAGCGCGAATGCGGGCCAGCACCGGGAAAATCAAAGCTCGTGCTTTGGCTTCGCtcat                                                             >  1:607016/1‑150 (MQ=255)
                                                                                         ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatca                                                           <  1:17989/150‑1 (MQ=255)
                                                                                          cAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCAATCCCCCCAAAGCCCGGTCTACACCAAAGCCGGCCAAAACGCGAATGCGGGCCCGCCCCGGGAAATTCAACGCTCCTGCTTTTGCTTCGCtcatcat                                                          >  2:981634/1‑150 (MQ=255)
                                                                                           aaCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGGATACACCCGCGCCACTCCCCTCCATGCCGGGTCTACACCAACGCCGGCAAAAGCGCGGAAGCGGGCCCGCACCGGGGAAATCAACGCCCGTGCTTTTGCTTCTCtcatcct                                                          >  2:530809/1‑147 (MQ=255)
                                                                                           aaCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTTCGGCACAATTCCCCTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcatc                                                         <  2:475451/150‑1 (MQ=255)
                                                                                             cTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCCCCATTTCCCTCAATATCCGGTCTACACCAACCCTGGCCAATATGCGAATGCGGGGCCGCACGGGGAAAATCAAAGCTCGTGCTTTTGCTTTGCTcctcatcat                                                       >  2:307685/1‑150 (MQ=255)
                                                                                                ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCAtt                                                      <  1:89989/148‑1 (MQ=255)
                                                                                                 cAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTAcc                                                   <  1:37693/150‑1 (MQ=255)
                                                                                                  aGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCg                                                  <  1:1028377/150‑1 (MQ=255)
                                                                                                        cacGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTg                                            <  1:322159/150‑1 (MQ=255)
                                                                                                          cGCTCCAGGCAACGGCGGGTCACATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgcg                                          <  1:766919/150‑1 (MQ=255)
                                                                                                               cAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATc                                     <  2:254178/150‑1 (MQ=255)
                                                                                                                        gcggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga                            >  2:770223/1‑150 (MQ=255)
                                                                                                                         cggGTCGCATACTTCGGCGCAATCCGGATCAATGCGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat                           <  2:853155/150‑1 (MQ=255)
                                                                                                                          gggTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat                           >  1:354754/1‑149 (MQ=255)
                                                                                                                             tCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGatgatg                       <  1:307685/150‑1 (MQ=255)
                                                                                                                                cATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGcgt                    <  1:478986/150‑3 (MQ=255)
                                                                                                                                      ccGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCg              <  1:833867/150‑1 (MQ=255)
                                                                                                                                       cGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGc             >  1:393206/1‑150 (MQ=255)
                                                                                                                                           gcAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGt          >  2:603141/1‑149 (MQ=255)
                                                                                                                                           gcAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGt          >  2:115234/1‑149 (MQ=255)
                                                                                                                                              aaTTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTcac       <  1:522675/150‑1 (MQ=255)
                                                                                                                                              aaTTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTca        >  1:732816/1‑149 (MQ=255)
                                                                                                                                               aTTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTcaca      <  1:13882/150‑1 (MQ=255)
                                                                                                                                                ttCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAg     <  2:556389/150‑1 (MQ=255)
                                                                                                                                                atCCCATCAAAGAAGGAACTAAACAAACGCAGGCAAAAGCGCGAATGCGGGCCAGAAAGGGTAAAAACAACGAAAGAGATATGGATTCGCTAAAAAAAAATAAAGAAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAg     <  1:886382/149‑1 (MQ=255)
                                                                                                                                                   ccATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGaa   <  1:530809/149‑1 (MQ=255)
                                                                                                                                                   ccATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAc  <  1:125864/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAC  >  CP009273/2977941‑2978236

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