Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,714,436 2 bp→GC coding (261‑262/660 nt) yiaD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,4360CG100.0% 356.7 / NA 133N87K (AAC→AAGyiaDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (60/73);  total (60/73)
*NC_0009133,714,4370GC99.2% 330.5 / NA 132V88L (GTG→CTG) yiaDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (58/73);  minor base A (1/0);  total (59/73)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

TGATTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAACGTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGC  >  NC_000913/3714366‑3714506
                                                                      ||                                                                     
  aTTGGCGCGGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTAGATGGATGTGCAGGAAGCGAAGCTGcgcg                                                                 <  2:4674253/76‑1 (MQ=255)
  aTTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGcgcg                                                                 <  2:4444068/76‑1 (MQ=255)
   ttGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGa                                                                <  2:719917/76‑1 (MQ=255)
    tGGCGCAGCAGCAGGCGCAGCTCTGGGGGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGAc                                                               <  2:715594/76‑1 (MQ=255)
    tGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTCCATGGATGTGCCGGAAGCGAAGCTGcgcgcc                                                               <  2:2261011/76‑3 (MQ=255)
    tGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGAc                                                               >  1:4850976/1‑76 (MQ=255)
     ggCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGGGACa                                                              >  2:2578619/1‑76 (MQ=255)
       cgcAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGACGCTGCGCGACaaa                                                            >  2:2042788/1‑76 (MQ=255)
        gcagcagcaGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACaaaa                                                           <  2:1726977/76‑1 (MQ=255)
        gcagcagcaGGCGCAGCTCTGGGCGGCGGCGTTGGTTACTACATGGATGTGCAGGAAGCGAAGCTGCGCGACaaaa                                                           >  1:2859523/1‑76 (MQ=255)
          agcagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACAAGGATGTGGAGGAAGCGACGCTGCGCGGCaaaaac                                                         >  2:458713/1‑74 (MQ=255)
           gcagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgc                                                        >  1:1307599/1‑76 (MQ=255)
           gcagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgc                                                        >  1:2740729/1‑76 (MQ=255)
            cagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcg                                                       <  1:1888121/76‑1 (MQ=255)
             ngcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGAGAAAAAATtcgg                                                      >  1:2347168/2‑75 (MQ=255)
             agcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcgc                                                      >  1:4726199/1‑76 (MQ=255)
             agcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcgc                                                      <  1:3204734/76‑1 (MQ=255)
             agcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcgc                                                      <  1:5527334/76‑1 (MQ=255)
             agcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcgc                                                      <  1:1100157/76‑1 (MQ=255)
               cagGCGCAGCTCTGGGCGGCTNCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCgg                                                    >  1:4516472/1‑76 (MQ=255)
                agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGc                                                   >  1:1686904/1‑76 (MQ=255)
                agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGc                                                   >  2:2489906/1‑76 (MQ=255)
                agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGc                                                   <  2:3572198/76‑1 (MQ=255)
                agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGc                                                   >  2:4180984/1‑76 (MQ=255)
                 ggCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCa                                                  <  2:3857185/76‑1 (MQ=255)
                  gcgcAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCAc                                                 >  1:5822404/1‑76 (MQ=255)
                   cgcAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACt                                                <  1:3785811/76‑1 (MQ=255)
                   cgcAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACt                                                <  2:2928342/76‑1 (MQ=255)
                     cAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTgg                                              <  1:763248/76‑1 (MQ=255)
                     cAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTgg                                              <  1:4343306/76‑1 (MQ=255)
                      aGCTCTGGGCGGCGGCGTTGGTTATTCCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGgt                                             <  1:3446726/76‑1 (MQ=255)
                      aGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGgt                                             >  1:2937011/1‑76 (MQ=255)
                        ctctGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGgtgt                                           <  2:3919976/76‑1 (MQ=255)
                        ctctGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGgtgt                                           >  1:2727400/1‑76 (MQ=255)
                         tctGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGtt                                          >  1:2759500/1‑76 (MQ=255)
                           tGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAg                                        <  2:4544855/76‑1 (MQ=255)
                           ggggCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAg                                        <  2:919446/75‑1 (MQ=255)
                              gcggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCCCCGGTGTTAGCGt                                     >  2:5270229/1‑76 (MQ=255)
                              gcggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGt                                     <  1:4326932/76‑1 (MQ=255)
                              gcggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGt                                     <  2:3869454/76‑1 (MQ=255)
                              gcggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGt                                     <  2:5466089/76‑1 (MQ=255)
                                ggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTaa                                   >  1:5747634/1‑76 (MQ=255)
                                 gcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAc                                  >  2:3238568/1‑76 (MQ=255)
                                 gcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAc                                  >  1:4529230/1‑76 (MQ=255)
                                  gggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAcc                                 <  2:2545051/75‑1 (MQ=255)
                                  cggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAcc                                 <  2:1336058/76‑1 (MQ=255)
                                  cggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAcc                                 <  1:3665598/76‑1 (MQ=255)
                                  cggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAcc                                 <  2:3807000/76‑1 (MQ=255)
                                  cggcgTTGGGTATTACATCGATGTGCAGGAAGAGACGCCGCGCGACAAAACGCGCGGCACTGGTCACAGCGTAAcc                                 >  2:2411976/1‑76 (MQ=255)
                                   ggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCTCGGCACTGGTGTTAGCGTAAccc                                >  1:5523188/1‑76 (MQ=255)
                                   ggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAccc                                >  1:367104/1‑76 (MQ=255)
                                    gcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCg                               <  1:411427/76‑1 (MQ=255)
                                    gcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCg                               <  2:4409287/76‑1 (MQ=255)
                                    gcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCg                               <  1:2785507/76‑1 (MQ=255)
                                     ggTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGc                              <  2:4842883/75‑1 (MQ=255)
                                     cgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGc                              <  2:1027725/76‑1 (MQ=255)
                                     cgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGc                              <  1:2028765/76‑1 (MQ=255)
                                     cgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGc                              <  1:2809486/76‑1 (MQ=255)
                                      gTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCa                             <  2:3053606/76‑1 (MQ=255)
                                       ttGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAg                            <  2:611777/76‑1 (MQ=255)
                                       ttGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAg                            >  2:3145571/1‑76 (MQ=255)
                                        tGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCg                               <  1:821469/72‑1 (MQ=255)
                                        tGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGc                           <  2:2232800/76‑1 (MQ=255)
                                        tGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGc                           <  1:1562759/76‑1 (MQ=255)
                                         gggtattaCAGGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCg                          <  2:923024/73‑1 (MQ=255)
                                          gTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCgg                         <  1:4345038/76‑1 (MQ=255)
                                          gTTATGACAGGGATGGGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCgg                         <  2:5757478/76‑1 (MQ=255)
                                           ttattaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCggg                        <  2:2349928/76‑1 (MQ=255)
                                           ttattaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCggg                        <  2:4850975/76‑1 (MQ=255)
                                           ttattaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCggg                        >  1:3572894/1‑76 (MQ=255)
                                           ttattaCATGGATATGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCggg                        <  1:2190019/76‑1 (MQ=255)
                                            tattaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCgggg                       >  2:4777950/1‑76 (MQ=255)
                                              ttaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGAt                     <  2:2225122/76‑1 (MQ=255)
                                              ttaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGAt                     <  1:3847316/76‑1 (MQ=255)
                                              ttaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGCGTTAGCGGCACCCGCAGCggggct                     >  2:410591/1‑74 (MQ=255)
                                               taCGTGGATGTGCAGGAAGCGAAGCCGCGCGCCAAAATGCGCGGCGCTGGTGTTAGCGTATCCCGCAGCGGGGATa                    <  2:1762506/76‑1 (MQ=255)
                                               taCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGGGTTAGCGGCACCCGCAGCggggcga                    >  2:450876/1‑73 (MQ=255)
                                                aCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATaa                   >  2:1967431/1‑76 (MQ=255)
                                                aCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATaa                   <  1:3521683/76‑1 (MQ=255)
                                                 cATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAAc                  >  1:2013242/1‑76 (MQ=255)
                                                   tGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACAt                >  2:2801262/1‑76 (MQ=255)
                                                   tGGAGGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACAt                <  2:543019/76‑1 (MQ=255)
                                                   gggATGTGCGGGCAGCGAAGCTTCGCGACAAAATGCGCGTCACTGGTGTTAGCGTCACCCGCAGCGTGGATAACAt                <  2:1391089/75‑1 (MQ=255)
                                                     gATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTa              >  1:3773489/1‑76 (MQ=255)
                                                     gATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTa              >  1:1401890/1‑76 (MQ=255)
                                                     gATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTa              <  1:4158002/76‑1 (MQ=255)
                                                       tgtgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATc            >  1:5485695/1‑76 (MQ=255)
                                                       tgtgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATc            >  2:3377378/1‑76 (MQ=255)
                                                        gtgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATcc           <  2:5747633/76‑1 (MQ=255)
                                                          gCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTc         >  2:3065917/1‑76 (MQ=255)
                                                           cAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGGTAGCGTAACCCGCAGCGGGGATACAATTATCCCCa        >  2:5547530/1‑76 (MQ=255)
                                                            aGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCaa       >  1:565387/1‑76 (MQ=255)
                                                            aGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCaa       <  2:4146790/76‑1 (MQ=255)
                                                             ggAAGCGAAGCTGCGCGACACAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAat      <  1:1744881/76‑1 (MQ=255)
                                                             ggAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAat      <  2:5198356/76‑1 (MQ=255)
                                                             ggAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAat      <  1:3795625/76‑1 (MQ=255)
                                                              gaagcgaagcTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAata     >  1:4472492/1‑76 (MQ=255)
                                                              gaagcgaagcTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGAGGGGATAACATTATCCTCAata     >  1:2497395/1‑76 (MQ=255)
                                                               aaGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAatat    >  2:5177655/1‑76 (MQ=255)
                                                                 gCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGc  <  2:4299403/76‑1 (MQ=255)
                                                                      ||                                                                     
TGATTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAACGTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGC  >  NC_000913/3714366‑3714506

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