Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5720AC29.5% 72.7 / 12.2 44V556G (GTG→GGG) ybaLputative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (30/1);  new base C (0/13);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 2.70e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAT  >  NC_000913/501433‑501708
                                                                                                                                           |                                                                                                                                        
aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca                                                                                                                                          >  1:384056/1‑140 (MQ=255)
        ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                                                                  >  1:1228935/1‑140 (MQ=255)
                 cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt                                                                                                                         >  1:1094227/1‑140 (MQ=255)
                         gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                 >  2:1071676/1‑140 (MQ=255)
                            tccccgccgcTTCTGGCTTTTGCCTCTCGATCCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTCCACCCGCCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCcag                                                                                                              <  2:1278440/136‑1 (MQ=255)
                            tccccgccgcCTCTGGCTTTGGGCTCTCGATGCCTGAGGCGACGCTGGCCCGTTTTTTTACCCTTACATTTTTCTAATTTTTTACCCCCGACACCCTCCTTTTCCCCCTTCACCCCCTCCCCCGCCGGTGGCGTTTCcag                                                                                                              <  2:709111/136‑1 (MQ=255)
                             ccccgccgcATCTGTCTTTTGGCTCTCGTTGCCTGATGCGCCGCGGGCGCGTTTTATCATCCCTAAATTTTTTTCATTTTTTACACCGGGCAACCCCCTTTTCCCCCTTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                                                                             <  2:1228935/137‑1 (MQ=255)
                              cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                            >  1:858143/1‑140 (MQ=255)
                                  ccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc                                                                                                        >  2:1352339/1‑140 (MQ=255)
                                     cATCTGTCTTTTGCCTCCCGATGCGTTACGCGACGCTGGCGCGTCTTATCTTGCCTACATTTTTTTAATATTTTCCACCCAGCAACCACCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                     <  2:672693/140‑1 (MQ=255)
                                       tCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                   <  2:946075/140‑1 (MQ=38)
                                       tCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                   >  1:299325/1‑140 (MQ=255)
                                                     ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgctgtctc                                                                                     >  1:373743/1‑133 (MQ=255)
                                                          tGCCTGTCGCGCCGCTGGCGCCTTTTATCATCCCTACTTTTTTTTCATATTTTACTCCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACg                                                                                <  2:384056/140‑1 (MQ=255)
                                                           gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                               <  1:1213230/140‑1 (MQ=38)
                                                            cttagtGCGACGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATTTTTCACACCCGCCAACCCCATTTTCCCCCAACACCCCCACCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                              <  2:589101/135‑1 (MQ=255)
                                                            ccTGATGCGCCGCTGGCGCTTCTTATCATCCCTACTTTTTTTTCATATTTTCCCCCCGCCAACCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                              <  2:689881/140‑1 (MQ=255)
                                                            ccTGATGCGCCGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                              <  2:124648/140‑1 (MQ=255)
                                                               gATGCGGCGCTGGCGCGTCTTATCATGCCTATATTTTTTTAATATTTCACATCCGGCAACCCCCTTTTCCCCCCTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                           <  2:1269747/140‑1 (MQ=255)
                                                                aTGCGCCGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACCCCCGCCACCCACCTTTTCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                          <  2:1123462/140‑1 (MQ=255)
                                                                aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                          >  1:627295/1‑140 (MQ=255)
                                                                    gACGCTGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                      <  1:1293188/140‑1 (MQ=255)
                                                                     aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                     >  1:84965/1‑140 (MQ=255)
                                                                            cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                              >  1:662011/1‑140 (MQ=255)
                                                                                     tCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                     >  1:1261602/1‑140 (MQ=255)
                                                                                     tCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                     >  2:1155740/1‑140 (MQ=255)
                                                                                      cATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                                    <  2:858143/140‑1 (MQ=255)
                                                                                         gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                 >  1:1090607/1‑140 (MQ=255)
                                                                                         cccTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                 <  2:299325/139‑1 (MQ=255)
                                                                                            tACATTTTTTTCATATTTTACACCCGAAACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                              <  2:1161299/140‑1 (MQ=255)
                                                                                            tACATATTTTTCATATTTTACATCCGGAAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                              >  1:1161299/1‑140 (MQ=255)
                                                                                             aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                             >  1:336832/1‑140 (MQ=255)
                                                                                              cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                            >  1:1268855/1‑140 (MQ=255)
                                                                                              cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                            >  2:517333/1‑140 (MQ=255)
                                                                                              cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                            >  1:799132/1‑140 (MQ=255)
                                                                                               atatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                           >  2:479520/1‑140 (MQ=255)
                                                                                                  tttttCATTTTTTCCACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                        <  2:84965/140‑1 (MQ=255)
                                                                                                   ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                       >  1:832145/1‑140 (MQ=255)
                                                                                                     ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                     >  1:160138/1‑140 (MQ=255)
                                                                                                      tCATTTTTTACATCCGGCACCCACCGTTTCCCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                    <  1:1352339/140‑1 (MQ=255)
                                                                                                        atatTTTACATCCGGCAACCCCCGTTTCCCCCTTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                  <  2:627295/140‑1 (MQ=255)
                                                                                                            tttACACCCGGCAACCCCCTTTCCCCCCGCCACCCCCTCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                              <  2:720369/140‑1 (MQ=255)
                                                                                                              tACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                            <  1:713638/140‑1 (MQ=255)
                                                                                                              tACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                            <  2:1455550/140‑1 (MQ=255)
                                                                                                              tACACCCGGCAACCCCCTTTTCCCCCCTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                            <  2:585577/140‑1 (MQ=255)
                                                                                                                  ccgggCACCCACCGTTACCCCCGTCCCCCCCTCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                        <  2:925176/137‑1 (MQ=255)
                                                                                                                    cGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa                      >  1:1356687/1‑140 (MQ=255)
                                                                                                                     ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGAGATATACGCCACTTCATCGTCAt                     >  1:925176/1‑140 (MQ=255)
                                                                                                                             ctttttCCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg             <  2:1145987/137‑1 (MQ=255)
                                                                                                                             ccGTTTCCCCCGCCCCCCCCCCAACCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg             <  2:1135772/140‑1 (MQ=255)
                                                                                                                             ccGTTTACCCCGTCACCACCTCACCCGCCGGGGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCCTCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGTgcg             >  2:165919/1‑140 (MQ=255)
                                                                                                                             ccGTTCCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg             <  2:1268855/140‑1 (MQ=255)
                                                                                                                              ctttccccccGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc            <  2:160138/134‑1 (MQ=255)
                                                                                                                              cGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc            <  1:517333/140‑1 (MQ=255)
                                                                                                                               ttttCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg           <  1:1512111/139‑1 (MQ=255)
                                                                                                                               gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGTAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg          >  2:1094833/1‑140 (MQ=255)
                                                                                                                                tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg          >  2:560345/1‑140 (MQ=255)
                                                                                                                                tttACCCCGTCACCAACTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg          >  2:133798/1‑140 (MQ=255)
                                                                                                                                  tcccccGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc        <  2:600691/138‑1 (MQ=255)
                                                                                                                                     cccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat     >  1:1145987/1‑140 (MQ=255)
                                                                                                                                     cccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat     >  1:1575424/1‑140 (MQ=255)
                                                                                                                                        gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat  >  1:1135772/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                        
ATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAT  >  NC_000913/501433‑501708

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