Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,420,807 +C coding (124/423 nt) ydaT → Rac prophage; protein YdaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,420,8031.C98.2% 202.6 / ‑3.2 55coding (120/423 nt)ydaTRac prophage; protein YdaT
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (26/28);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCT  >  NC_000913/1420664‑1420937
                                                                                                                                             |                                                                                                                                      
cATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTg                                                                                                                                         <  1:334382/140‑1 (MQ=255)
  tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGcc                                                                                                                                      >  2:1212509/1‑140 (MQ=255)
   cAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGccc                                                                                                                                     >  1:179099/1‑140 (MQ=255)
    aaGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGcccc                                                                                                                                    <  2:914490/140‑1 (MQ=255)
    aaGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGcccc                                                                                                                                    >  1:1171404/1‑140 (MQ=255)
           aaaTGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAGTATCTGCGCCAGGGGGGTGGCGAATTGCCCCCTGGTTc                                                                                                                              >  2:269112/1‑140 (MQ=255)
                aTTCATGAAAATCAAGCATGATCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGc                                                                                                                         >  2:119780/1‑140 (MQ=255)
                aTTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGc                                                                                                                         <  2:1237609/140‑1 (MQ=255)
                 ttCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCa                                                                                                                        <  2:132613/140‑1 (MQ=255)
                 ttCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCa                                                                                                                        <  1:379310/140‑1 (MQ=255)
                  tCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCaa                                                                                                                       >  1:190271/1‑140 (MQ=255)
                    aTGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTcc                                                                                                                             <  2:351415/132‑1 (MQ=255)
                    aTGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTcc                                                                                                                             >  1:351415/1‑132 (MQ=255)
                        aaaTCAAGCATGATCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTg                                                                                                                 <  1:119780/140‑1 (MQ=255)
                               gCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTTCTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAAt                                                                                                          <  2:483646/140‑1 (MQ=255)
                               gCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAAt                                                                                                          <  2:544831/140‑1 (MQ=255)
                                  tGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                                       <  2:1231832/140‑1 (MQ=255)
                                    aGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATcagca                                                                                                     <  2:1268726/140‑1 (MQ=255)
                                    aGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGATTCCAGGCAAGGACTGGAACAATcagca                                                                                                     >  1:968334/1‑140 (MQ=255)
                                     gCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                                   >  1:40409/1‑140 (MQ=255)
                                      cacaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                                   <  1:765860/140‑1 (MQ=255)
                                      cacaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                                   <  2:1171404/140‑1 (MQ=255)
                                      cacaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                                   <  2:179099/140‑1 (MQ=255)
                                        caTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAat                                                                                                 <  1:95851/140‑1 (MQ=255)
                                         aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAata                                                                                                <  1:227280/140‑1 (MQ=255)
                                          tCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAatat                                                                                               >  2:1261341/1‑140 (MQ=255)
                                               tCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAAtatctatc                                                                                          <  1:532637/140‑1 (MQ=255)
                                                 aGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAc                                                                                        >  1:899389/1‑140 (MQ=255)
                                                 aGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAc                                                                                        <  1:1068498/140‑1 (MQ=255)
                                                 aGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAc                                                                                        <  1:564677/140‑1 (MQ=255)
                                                       gttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATcagcagg                                                                                                   <  1:77160/123‑2 (MQ=255)
                                                       gttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATcagcagg                                                                                                   >  2:77160/1‑122 (MQ=255)
                                                        tttGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgttggt                                                                                 >  2:877519/1‑140 (MQ=255)
                                                                      cGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACg                                                                   >  1:378294/1‑140 (MQ=255)
                                                                        aaaaaGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAGTATCTGCGCCAGGGGGGTGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGaa                                                                 <  1:269112/140‑1 (MQ=255)
                                                                               gcaggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCaa                                                          >  2:436423/1‑140 (MQ=255)
                                                                                   gcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAagag                                                      >  1:777052/1‑140 (MQ=255)
                                                                                   gcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAagag                                                      >  2:425107/1‑140 (MQ=255)
                                                                                           aGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAAtt                                              <  2:1243683/140‑1 (MQ=255)
                                                                                            gCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTc                                             <  1:1119780/140‑1 (MQ=255)
                                                                                             ccAATGCAATTACTGNAGAATANCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCa                                            >  1:659167/1‑140 (MQ=255)
                                                                                                 tGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAg                                        >  1:460692/1‑140 (MQ=255)
                                                                                                 tGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAg                                        <  2:40409/140‑1 (MQ=255)
                                                                                                              gaaTATCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCaa                           <  1:1261341/140‑1 (MQ=255)
                                                                                                                atatCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAAtt                         >  1:858987/1‑140 (MQ=255)
                                                                                                                  atCTGCGCCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCt                       <  1:877519/140‑1 (MQ=255)
                                                                                                                      gcgcCAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCa                   >  1:827138/1‑140 (MQ=255)
                                                                                                                         ccAGGGGGGCGGCGAATTACCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATc                >  1:397080/1‑140 (MQ=255)
                                                                                                                         ccAGGGGGGCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATc                >  2:796453/1‑140 (MQ=255)
                                                                                                                         ccAGGGGGGCGGCGAATTGCCCCCTGGCTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCaa                           <  1:617518/129‑1 (MQ=255)
                                                                                                                         ccAGGGGGGCGGCGAATTGCCCCCTGGCTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCaa                           >  2:617518/1‑129 (MQ=255)
                                                                                                                            ggggggCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCtt             >  2:65269/1‑140 (MQ=255)
                                                                                                                             gggggCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTc            >  2:1279062/1‑140 (MQ=255)
                                                                                                                             gggggCGGCGAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTc            >  2:1040342/1‑140 (MQ=255)
                                                                                                                                  cggcgAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCcg          >  2:520911/1‑137 (MQ=255)
                                                                                                                                  cggcgAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCcg          <  1:520911/137‑1 (MQ=255)
                                                                                                                                    gcgAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGa     >  1:199338/1‑140 (MQ=255)
                                                                                                                                    gcgAATTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGa     >  2:219329/1‑140 (MQ=255)
                                                                                                                                       aaTTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCt  <  2:827138/140‑1 (MQ=255)
                                                                                                                                       aaTTGCCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCt  <  1:436423/140‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                      
CATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCT  >  NC_000913/1420664‑1420937

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