Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 714,998 T→A 7.0% intergenic (+64/+13) ispH → / ← PP_0607 4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947714,9980TA 7.0% 169.7 / 2.8 57intergenic (+64/+13)ispH/PP_06074‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base A (2/2);  total (28/29)
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.76e-01

ACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCC  >  NC_002947/714851‑715143
                                                                                                                                                   |                                                                                                                                                    
aCTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTg                                                                                                                                                     <  3:5087/149‑1 (MQ=255)
    gcgcTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGAGACACAAGACCGCTCCTGCAAt                                                                                                                                                 >  8:303279/1‑149 (MQ=255)
                  gACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTagcag                                                                                                                                 >  1:46048/1‑149 (MQ=255)
                  gACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTagcg                                                                                                                                  >  6:230921/1‑148 (MQ=255)
                    cGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGCGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCg                                                                                                                              >  1:255626/1‑149 (MQ=255)
                       gcgcGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGt                                                                                                                            >  3:348908/1‑148 (MQ=255)
                        cgcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCa                                                                                                                          <  5:181553/149‑1 (MQ=255)
                          cgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGCGCGGTCAtt                                                                                                                        <  5:269161/149‑1 (MQ=255)
                             agagAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCt                                                                                                                     <  8:230880/149‑1 (MQ=255)
                               agaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCa                                                                                                                   <  3:325700/149‑1 (MQ=255)
                                gaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTAGCGCGGTCATTGcagga                                                                                                                  <  5:230921/149‑5 (MQ=255)
                                gaACATCACCTTATCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAg                                                                                                                  <  7:63015/149‑1 (MQ=255)
                                     tCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcct                                                                                                             >  7:306998/1‑149 (MQ=255)
                                     tCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcct                                                                                                             <  2:46048/149‑1 (MQ=255)
                                      cACCTTTTCAATGCCCAAAGAGCTGCTGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTAGCGCGGTCAtt                                                                                                                        <  5:194068/137‑1 (MQ=255)
                                      cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTAGCGCGGTCAtt                                                                                                                        >  6:194068/1‑137 (MQ=255)
                                            ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                       <  7:329760/148‑1 (MQ=255)
                                            ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                       <  6:162495/148‑1 (MQ=255)
                                                   ccgatagagGTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCTGTCATTGCTCAGGTCCTCCGTAAGc                                                                                                     <  8:610/138‑1 (MQ=255)
                                                   cccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGTGGCTGCTGTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGATCAGGTCCTCCGTAAGc                                                                                                     >  7:610/1‑143 (MQ=255)
                                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGTCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                              <  8:183061/149‑1 (MQ=255)
                                                     cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                             >  6:360816/1‑149 (MQ=255)
                                                      aaaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                             >  5:180201/1‑148 (MQ=255)
                                                       aaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGActc                                                                                           <  7:303279/149‑1 (MQ=255)
                                                              gCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCt                                                                                    <  4:3380/149‑1 (MQ=255)
                                                               cgtgTTCTCTGGCTGATCTGACGGGCAGTACCGCGGACAATGCAGTGGCTGCTTTGGAGCCGTTGCGAGACACAAGACCTCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAATCATCCTGACTCTGCCGCtt                                                                                   <  4:348908/146‑1 (MQ=255)
                                                               cGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGc                                                                                     <  2:255626/147‑1 (MQ=255)
                                                               cGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCtt                                                                                   <  3:199167/149‑1 (MQ=255)
                                                                      gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACAAAAGACCGCTCCTGCAATGACCCCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGc                                                                              >  8:151279/1‑147 (MQ=255)
                                                                       ctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGCGGCTGCTTTGCAGCCCTTGCGCGACACAACACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAACCAACCTCACTCTGCCGCTTGATCCCac                                                                           >  8:291972/1‑148 (MQ=255)
                                                                          gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGca                                                                         <  7:291972/148‑1 (MQ=255)
                                                                                    cGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGatat                                                              >  6:140626/1‑149 (MQ=255)
                                                                                           tACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                       <  6:180201/149‑1 (MQ=255)
                                                                                           tACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCCGc                                                       >  4:176040/1‑149 (MQ=255)
                                                                                             ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcaccac                                                                     <  8:220725/133‑1 (MQ=255)
                                                                                             ccgccgCCAATGCAGGGGCTGCTTGGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcaccac                                                                     >  7:220725/1‑133 (MQ=255)
                                                                                              cgccgcCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                              >  2:302283/1‑107 (MQ=255)
                                                                                              cgccgcCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                              <  1:302283/107‑1 (MQ=255)
                                                                                                   ccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAgcg                                               <  5:360816/149‑1 (MQ=255)
                                                                                                     aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTGCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCCGCACCTGATATCGGCTGCGAGCCgcggc                                             >  2:45645/1‑149 (MQ=255)
                                                                                                     aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCt                                                         >  8:69807/1‑137 (MQ=255)
                                                                                                     aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCt                                                         <  7:69807/137‑1 (MQ=255)
                                                                                                          aagggCTGCAAAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCTGGACATTGCTCAGGTCCTCCGTAAGCAACATGACTCTGCCGCTTGATGCCAGCACAACCTGATATCGGATGCGAGCAGCGGCGCGt                                         >  1:286289/3‑149 (MQ=255)
                                                                                                            gggCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCg                                      >  7:71173/1‑149 (MQ=255)
                                                                                                                     tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATAGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCACCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTc                              >  1:239654/1‑149 (MQ=255)
                                                                                                                      gCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                           >  6:225660/1‑118 (MQ=255)
                                                                                                                      gCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                           <  5:225660/118‑1 (MQ=255)
                                                                                                                      gCAGCCCTAGCGCGACACAAGCCCGCTAATGCAAAGACCGCGCGAGCAGCGGTCATTGCTCAGGGCCTCCGTAACCAACCTGCCTCTCCCCCTAGATCCcaccac                                                                        >  8:153745/1‑105 (MQ=255)
                                                                                                                            cTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGc                                     <  2:282972/134‑1 (MQ=255)
                                                                                                                            cTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGc                                     >  1:282972/1‑134 (MQ=255)
                                                                                                                            cTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTg                      >  5:175718/1‑149 (MQ=255)
                                                                                                                             ttGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCGTGATGCCAGCCCCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgc                     >  7:168532/1‑149 (MQ=255)
                                                                                                                               gcgcGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgcg                    <  8:306998/148‑1 (MQ=255)
                                                                                                                               gcgcGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgc                     <  2:239654/147‑1 (MQ=255)
                                                                                                                                      acaAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAAc              >  2:274110/1‑147 (MQ=255)
                                                                                                                                             cGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTg     >  2:164546/1‑149 (MQ=255)
                                                                                                                                                tCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc  <  8:71173/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCC  >  NC_002947/714851‑715143

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: