Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,712,106 A→G intergenic (‑99/+130) yhjX ← / ← yhjY putative pyruvate transporter/conserved protein YhjY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,712,1060AG98.4% 219.0 / ‑4.7 63intergenic (‑99/+130)yhjX/yhjYputative pyruvate transporter/conserved protein YhjY
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (31/31);  total (32/31)
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

TGTCAGCCAGCGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGAACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAATTGTGACTGTTACGTTAA  >  NC_000913/3711969‑3712239
                                                                                                                                         |                                                                                                                                     
tgtCAGCCAGCGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGAc                                                                                                                                     >  1:259427/1‑140 (MQ=255)
        agcGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGaa                                                                                                                             <  1:820191/140‑1 (MQ=255)
        agcGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGaa                                                                                                                             >  1:524906/1‑140 (MQ=255)
         gcGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGaaa                                                                                                                            <  1:92521/140‑1 (MQ=255)
          cGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGGCGGAATGAAAt                                                                                                                           >  1:294377/1‑140 (MQ=255)
          cGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAAt                                                                                                                           >  2:235796/1‑140 (MQ=255)
          cGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAAt                                                                                                                           >  1:210713/1‑140 (MQ=255)
          cGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAAt                                                                                                                           >  1:2467/1‑140 (MQ=255)
              tACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCt                                                                                                                       <  2:114938/140‑1 (MQ=255)
                    gATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAAACTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGGAAAAAGGCATGAAAAGCCATCAGGACGGAATGAAATGCCTTAGTTc                                                                                                                 >  2:92521/1‑140 (MQ=255)
                     aTAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCa                                                                                                                <  2:536767/140‑1 (MQ=255)
                     aTAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCa                                                                                                                <  1:235796/140‑1 (MQ=255)
                           tGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGgaatg                                                                                                          >  2:109302/1‑140 (MQ=255)
                               ggTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGaatgaatg                                                                                                      >  2:859555/1‑140 (MQ=255)
                                   tCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGAcgc                                                                                                  >  2:785747/1‑140 (MQ=255)
                                    cATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGAcgcg                                                                                                 >  1:692751/1‑140 (MQ=255)
                                           gTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGa                                                                                                >  1:863780/1‑134 (MQ=255)
                                           gTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGa                                                                                                <  2:863780/134‑1 (MQ=255)
                                           gTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCCGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAatta                                                                                          >  1:696390/1‑140 (MQ=255)
                                            tATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCCGGAATGAATGACGCGACAattat                                                                                         >  2:591433/1‑140 (MQ=255)
                                                 cTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCAcc                                                                                    <  2:784201/140‑1 (MQ=255)
                                                  tGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCGATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCa                                                                                   >  2:758012/1‑140 (MQ=255)
                                                         ataaAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATttat                                                                            >  1:587227/1‑140 (MQ=255)
                                                               aGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTg                                                                      >  2:400953/1‑140 (MQ=255)
                                                                 gcaaaaTCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGt                                                                    <  2:587227/138‑1 (MQ=255)
                                                                 gAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCa                                                                                   >  1:439761/1‑125 (MQ=255)
                                                                 gAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCa                                                                                   <  2:439761/125‑1 (MQ=255)
                                                                 gAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGt                                                                    >  2:124590/1‑140 (MQ=255)
                                                                  aaaaaTCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTa                                                                   <  1:400953/140‑1 (MQ=255)
                                                                  aaaaaTCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTa                                                                   <  1:327986/140‑1 (MQ=255)
                                                                    aaaTCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAg                                                                 <  1:17295/140‑1 (MQ=255)
                                                                          ctGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCa                                                           <  1:109302/140‑1 (MQ=255)
                                                                           tGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCAc                                                          <  1:121353/140‑1 (MQ=255)
                                                                           tGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCAc                                                          >  2:740456/1‑140 (MQ=255)
                                                                            gCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCAcc                                                         <  2:131441/140‑1 (MQ=255)
                                                                                aCAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTa                                                                   >  1:824714/1‑126 (MQ=255)
                                                                                aCAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTa                                                                   <  2:824714/126‑1 (MQ=255)
                                                                                aCAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACt                                                     <  1:740456/140‑1 (MQ=255)
                                                                                      aTAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGcc                                               <  2:839220/140‑1 (MQ=255)
                                                                                       tAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCca                                              >  1:128197/1‑140 (MQ=255)
                                                                                         aCGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCcaca                                            <  1:124590/140‑1 (MQ=255)
                                                                                          cGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACaa                                           <  1:785747/140‑1 (MQ=255)
                                                                                           gAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAAc                                          <  2:351056/140‑1 (MQ=255)
                                                                                             aGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTa                                        <  2:524906/140‑1 (MQ=255)
                                                                                               tGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCAccc                                                        <  1:59413/122‑1 (MQ=255)
                                                                                               tGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCAccc                                                        >  2:59413/1‑122 (MQ=255)
                                                                                               tGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGc                                                >  2:510648/1‑130 (MQ=255)
                                                                                               tGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGc                                                <  1:510648/130‑1 (MQ=255)
                                                                                               tGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTAtt                                      >  1:116602/1‑140 (MQ=255)
                                                                                               tGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTAtt                                      >  1:638464/1‑140 (MQ=255)
                                                                                                       aTGCCTGGGGAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTa                              <  1:648994/140‑1 (MQ=255)
                                                                                                             ggggAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGtt                        >  2:157625/1‑140 (MQ=255)
                                                                                                             ggggAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGtt                        >  2:744978/1‑140 (MQ=255)
                                                                                                              gggAAAAACGGCATGTAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGttt                       >  2:738676/1‑140 (MQ=255)
                                                                                                              gggAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGttt                       >  1:572829/1‑140 (MQ=255)
                                                                                                              gggAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGttt                       >  1:848148/1‑140 (MQ=255)
                                                                                                              gggAAAAACGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGttt                       >  2:242294/1‑140 (MQ=255)
                                                                                                                  aaaaCGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCCGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAAt                   <  1:591433/140‑1 (MQ=255)
                                                                                                                   aaaCGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAAtt                  <  1:140387/140‑1 (MQ=255)
                                                                                                                     aCGGCATGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGCCGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAATtgt                <  2:572829/140‑1 (MQ=255)
                                                                                                                           tGAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAATTGTGACTGt          <  2:116602/140‑1 (MQ=255)
                                                                                                                            gAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAATTGTGACTGtt         <  1:744978/140‑1 (MQ=255)
                                                                                                                            gAAATGCCATCAGGACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAATTGTGACTGtt         <  2:848148/140‑1 (MQ=255)
                                                                                                                                   cATCAGAACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAATTGTGACTGTTACGTTaa  >  2:267803/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                     
TGTCAGCCAGCGGGTACGCTGATAATTTGAAGGTGTCATGGCAGTATTCCTGCAGTAATAAAAAGGAAAAATCTCTGCCGACAACCATAACGAAGTGCACTTTATGCCTGGGGAAAAACGGCATGAAATGCCATCAGAACGGAATGAAATGCCTTAGTTCAGGAATGAATGACGCGACAATTATTCACCAAATTTATTATTTGGTAAGTTAGCACCCACTATCGCCACAACTATTGTTATTTAATGGTTTCAATTGTGACTGTTACGTTAA  >  NC_000913/3711969‑3712239

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