Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,787,510 T→A 6.5% F558I (TTC→ATC)  bamA_1 → Outer membrane protein assembly factor BamA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,787,5100TA 6.5% 195.8 / 4.0 62F558I (TTC→ATC) bamA_1Outer membrane protein assembly factor BamA
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base A (2/2);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CGTCAGCCTGGGCTACCCGATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGACCTGTCGTT  >  ACB122/1787377‑1787634
                                                                                                                                     |                                                                                                                            
cGTCAGCCTGGGCTACCCGATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACtt                                                                                                                             <  6:22414/135‑1 (MQ=255)
cGTCAGCCTGGGCTACCCGATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACtt                                                                                                                             <  4:297665/135‑1 (MQ=255)
       cTGGGCTACCCGATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAAc                                                                                                                      >  7:311163/1‑135 (MQ=255)
       cTGGGCTACCCGATCAGCGAAACCTCGCGCCTGACCTAAGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCGCGAGGACGAAGGTGACAACTTCCTGAAc                                                                                                                      >  2:378931/1‑135 (MQ=255)
        tGGGCTACCCGATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCATATACACTGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACt                                                                                                                     <  1:352701/135‑1 (MQ=255)
                ccGATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGTACAACATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAgg                                                                                                               <  1:186421/133‑1 (MQ=255)
                  gATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCt                                                                                                            >  7:84822/1‑134 (MQ=255)
                                cgcgCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACGGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCg                                                                                             <  5:316446/135‑1 (MQ=255)
                                  cgcCTGACCTACGGCCGGCGCGCGCAGCAGGACAAGCTCAAGACCTGCAAAAACACCGTTGACGAGAGCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGAGCGGCTGGTCCg                                                                                             <  6:365010/133‑1 (MQ=255)
                                  cgcCTGACCTACGCCCTCAGCGTGCAGAAGGACAAGAGCACGACCGGCAAATAAACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAATGCCGCGATCGGCTGGTCCg                                                                                             <  3:347478/133‑1 (MQ=255)
                                     cTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCg                                                                                        <  6:69323/135‑1 (MQ=255)
                                     cTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCg                                                                                        <  5:243233/135‑1 (MQ=255)
                                     cTGACCTACGGCCTGAGCGTGCAGCAGCACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCg                                                                                        <  3:185674/135‑1 (MQ=255)
                                          cTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCt                                                                                   <  8:72032/135‑1 (MQ=255)
                                          cTACGGCCTGAGCGTGCAGCAGGACAAGAGCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCt                                                                                   <  7:81591/135‑1 (MQ=255)
                                           tACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCTAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCt                                                                                   >  6:102857/1‑134 (MQ=255)
                                               gCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACa                                                                              <  3:263201/135‑1 (MQ=255)
                                               gCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACa                                                                              <  4:140738/135‑1 (MQ=255)
                                                     gCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGc                                                                         <  8:263422/134‑1 (MQ=255)
                                                       gTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTa                                                                      >  5:368034/1‑135 (MQ=255)
                                                       gTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTa                                                                      >  7:47330/1‑135 (MQ=255)
                                                       gTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTa                                                                      >  7:265813/1‑135 (MQ=255)
                                                       gTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCATGAACAAAGGCGTa                                                                      >  6:305217/1‑135 (MQ=255)
                                                         gcagcaGGACAAGATCAAGACCGGCTAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACt                                                                    <  5:102857/135‑1 (MQ=255)
                                                             cagGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCGGAACAAAGGCGTACTGGcc                                                                >  2:286501/1‑135 (MQ=255)
                                                                   aaGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGt                                                          >  5:175687/1‑135 (MQ=255)
                                                                   aaGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGt                                                          <  5:184993/135‑1 (MQ=255)
                                                                     gATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTgg                                                        <  2:87284/135‑1 (MQ=255)
                                                                     gATCAAGACCGGCAAATACACCGTTGACGAGATCTTAGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTgg                                                        <  2:113355/135‑1 (MQ=255)
                                                                        cAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCa                                                     >  2:8188/1‑135 (MQ=255)
                                                                         aaGACCTGCAAATCCTCCGTTGACGAGATCTTGGACTTCCTCGAGGAGGAAGGTGACATCTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAAGAAAGGCGTACTGGCCACCCGTGGTCAc                                                    <  1:372947/135‑1 (MQ=38)
                                                                          aGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACt                                                   >  7:290516/1‑135 (MQ=255)
                                                                             ccGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACACAGGCGTACTGGCCACCAGTGGTCACTCg                                                 >  6:276588/1‑134 (MQ=255)
                                                                             ccGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCg                                                 >  8:223210/1‑134 (MQ=255)
                                                                              cGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGc                                                >  6:17401/1‑134 (MQ=255)
                                                                                gCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCaaa                                             <  7:104266/135‑1 (MQ=255)
                                                                                gCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCaaa                                             <  8:57328/135‑1 (MQ=255)
                                                                                gCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCaaa                                             <  2:350198/135‑1 (MQ=255)
                                                                                gCAAATACACCGTTGACGAGAGCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGGCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCaaa                                             <  1:133523/135‑1 (MQ=255)
                                                                                 cAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAg                                            <  8:84822/135‑1 (MQ=255)
                                                                                   aaTACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGcc                                          >  8:219966/1‑135 (MQ=255)
                                                                                   aaTACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGcc                                          >  7:286002/1‑135 (MQ=255)
                                                                                   aaTACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCAGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGcc                                          >  8:176252/1‑135 (MQ=255)
                                                                                   aaTACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACATCCTGAAc                                                                                                                      >  2:246486/1‑59 (MQ=255)
                                                                                   aaTACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACATCCTGAAc                                                                                                                      <  1:241000/59‑1 (MQ=255)
                                                                                   aaTACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACATCCTGAAc                                                                                                                      <  1:246486/59‑1 (MQ=255)
                                                                                   aaTACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACATCCTGAAc                                                                                                                      >  2:241000/1‑59 (MQ=255)
                                                                                     taaccggTTGAGGAGAGCTTCGACTTCCTCGAGTAGGAAGGTGACAACTTACTGAACTTCAAGTCCTCGATGGTGTGGTCCGAGTAGACCCTGAACATAGGCGTACTGGCCACCCGTGGTCACTCGCATAGCCTg                                        <  7:360158/129‑1 (MQ=255)
                                                                                              gACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTTGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAg                               >  8:141838/1‑135 (MQ=255)
                                                                                                ctagaTCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGt                              <  2:176275/132‑1 (MQ=255)
                                                                                                    aTCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCATCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAACCCTGACCCTGGAGTCCAcc                         >  1:170417/1‑135 (MQ=255)
                                                                                                       ttCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCggt                      >  5:216960/1‑133 (MQ=255)
                                                                                                       ttCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGtt                      >  6:286250/1‑135 (MQ=255)
                                                                                                          gACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTcc                    <  5:286250/134‑1 (MQ=255)
                                                                                                           aCTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCgg                  <  4:360947/135‑1 (MQ=255)
                                                                                                           aCTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGAGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCgg                  <  6:44580/135‑1 (MQ=255)
                                                                                                           aCTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTAAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGAGGTCACTCTCAAAGCCTGACCCTGGAGTCCACCGTTCCgg                  >  2:268376/1‑135 (MQ=255)
                                                                                                                 tCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCGTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGc             <  1:121219/134‑1 (MQ=255)
                                                                                                                  cGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGa           >  6:382103/1‑135 (MQ=255)
                                                                                                                      gaggaAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGACCTg       <  6:175687/135‑1 (MQ=255)
                                                                                                                      gaggaAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGACCTg       <  8:286002/135‑1 (MQ=255)
                                                                                                                      gaggaAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGACCTg       <  8:128258/135‑1 (MQ=255)
                                                                                                                      gaggaAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGACCTg       <  2:99377/135‑1 (MQ=255)
                                                                                                                           aGGTGACAACTTCCTGAACTTCAAGGCCTTGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGACCTGTCGtt  <  7:141838/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                            
CGTCAGCCTGGGCTACCCGATCAGCGAAACCTCGCGCCTGACCTACGGCCTGAGCGTGCAGCAGGACAAGATCAAGACCGGCAAATACACCGTTGACGAGATCTTCGACTTCCTCGAGGAGGAAGGTGACAACTTCCTGAACTTCAAGGCCTCGATCGGCTGGTCCGAGTCGACCCTGAACAAAGGCGTACTGGCCACCCGTGGTCACTCGCAAAGCCTGACCCTGGAGTCCACCGTTCCGGGCAGCGACCTGTCGTT  >  ACB122/1787377‑1787634

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