Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,432,273 G→A A1075A (GCG→GCA stfR → Rac prophage, putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,432,2730GA98.6% 217.1 / NA 71A1075A (GCG→GCAstfRRac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/35);  minor base C (1/0);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGC  >  NC_000913/1432140‑1432403
                                                                                                                                     |                                                                                                                                  
cTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATAc                                                                                                                               >  2:481683/1‑139 (MQ=255)
     gCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCg                                                                                                                          <  1:136643/139‑1 (MQ=255)
         aCTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTAt                                                                                                                      <  2:144696/139‑1 (MQ=255)
            ccGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTgg                                                                                                                   <  1:340064/139‑1 (MQ=255)
                  cTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTc                                                                                                                           <  1:144887/125‑1 (MQ=255)
                  cTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTc                                                                                                                           >  2:144887/1‑125 (MQ=255)
                     gTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATAc                                                                                                          >  1:176722/1‑139 (MQ=255)
                         tGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCAc                                                                                                      <  1:364739/139‑1 (MQ=255)
                          gCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACt                                                                                                     >  2:363556/1‑139 (MQ=255)
                          gCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACt                                                                                                     >  2:324441/1‑139 (MQ=255)
                                         tCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                                <  1:772578/97‑1 (MQ=255)
                                         tCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                                >  2:772578/1‑97 (MQ=255)
                                              tgtgCATAATCAAAACTATGCAACATCGTCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACAt                                                                                 >  2:250839/1‑139 (MQ=255)
                                              tgtgCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACAt                                                                                 >  2:791093/1‑139 (MQ=255)
                                                  cATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                        >  2:307052/1‑128 (MQ=255)
                                                  cATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                        <  1:307052/128‑1 (MQ=255)
                                                  cATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGac                                                                             >  2:649255/1‑139 (MQ=255)
                                                   aTAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGaca                                                                            >  1:100107/1‑139 (MQ=255)
                                                     aaTCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGacaca                                                                          >  2:259587/1‑139 (MQ=255)
                                                       tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACAcc                                                                        >  1:51750/1‑139 (MQ=255)
                                                       tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACAcc                                                                        >  2:171935/1‑139 (MQ=255)
                                                               aTGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCGAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCCATTGGTTCACATGGACACACCATCACCGt                                                                >  2:543607/1‑139 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTgg                                                                                                                   >  1:681249/1‑87 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTgg                                                                                                                   <  2:681249/87‑1 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTg                                                                                                                 >  1:512985/1‑89 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTg                                                                                                                 <  2:512985/89‑1 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGtt                                                               >  2:375326/1‑139 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGtt                                                               >  1:462656/1‑139 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGtt                                                               >  2:288531/1‑139 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGtt                                                               >  2:599738/1‑139 (MQ=255)
                                                                tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACGGTCGGTATTGGTGCTCATACGCACTCCATTGCGATTGGTTCACATGGACACACCATCACCGtt                                                               >  2:484384/1‑139 (MQ=255)
                                                                  cAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTaa                                                             <  1:117843/139‑1 (MQ=255)
                                                                    aCATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACg                                                           <  1:649255/139‑1 (MQ=255)
                                                                     catcatCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgc                                                          <  1:453822/139‑1 (MQ=255)
                                                                       tcatcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACAAACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                              <  2:102835/117‑1 (MQ=255)
                                                                       tcatcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACAAACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                              >  1:102835/1‑117 (MQ=255)
                                                                               tgGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                                <  2:428615/139‑1 (MQ=255)
                                                                                ggCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAcg                                               >  1:543154/1‑139 (MQ=255)
                                                                                    cacaTACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGa                                           >  1:100227/1‑139 (MQ=255)
                                                                                        tACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAac                                       <  1:599738/139‑1 (MQ=37)
                                                                                        tACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAac                                       <  1:324441/139‑1 (MQ=37)
                                                                                         aCCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                      <  2:528379/139‑1 (MQ=37)
                                                                                                 actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCa                                                                                                       >  1:67905/1‑66 (MQ=255)
                                                                                                 actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCa                                                                                                       <  2:67905/66‑1 (MQ=255)
                                                                                                 actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                              >  2:348935/1‑91 (MQ=255)
                                                                                                 actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                              <  1:348935/91‑1 (MQ=255)
                                                                                                   tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCCCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                  >  2:584559/1‑117 (MQ=255)
                                                                                                   tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                  <  1:584559/117‑1 (MQ=25)
                                                                                                     tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTcaca                                                                                  >  2:10540/1‑83 (MQ=255)
                                                                                                     tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTcaca                                                                                  <  1:10540/83‑1 (MQ=255)
                                                                                                     tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                          >  2:267483/1‑139 (MQ=21)
                                                                                                     tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                          >  1:432064/1‑139 (MQ=25)
                                                                                                     tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                          <  2:119137/139‑1 (MQ=21)
                                                                                                        ggCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCg                       <  2:1606/139‑1 (MQ=21)
                                                                                                         gCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGc                      >  1:426161/1‑139 (MQ=25)
                                                                                                          cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                     <  1:288531/139‑1 (MQ=37)
                                                                                                          cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                     <  2:543154/139‑1 (MQ=25)
                                                                                                          cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                     <  1:267483/139‑1 (MQ=21)
                                                                                                          cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                     <  2:176722/139‑1 (MQ=21)
                                                                                                            ctgctgCAAGCGCAGGTGCACACGCACATACTGTCGGTATTCGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt                   <  1:481683/139‑1 (MQ=16)
                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                               <  2:557294/44‑1 (MQ=255)
                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                               >  1:557294/1‑44 (MQ=255)
                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCg                       >  1:202047/1‑132 (MQ=18)
                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCg                       <  2:202047/132‑1 (MQ=25)
                                                                                                                   aaGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                  <  1:621230/101‑1 (MQ=25)
                                                                                                                   aaGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                  >  2:621230/1‑101 (MQ=25)
                                                                                                                        cAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgt       <  1:259587/139‑1 (MQ=18)
                                                                                                                        cAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgt       <  2:765351/139‑1 (MQ=37)
                                                                                                                          ggTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACg                                                           <  2:694378/85‑1 (MQ=255)
                                                                                                                          ggTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACg                                                           >  1:694378/1‑85 (MQ=255)
                                                                                                                             gCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGc  <  2:100107/139‑1 (MQ=14)
                                                                                                                                     |                                                                                                                                  
CTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGC  >  NC_000913/1432140‑1432403

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