Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,564,444 T→C intergenic (+72/‑79) yffR → / → yffS CPZ‑55 prophage; uncharacterized protein YffR/CPZ‑55 prophage; uncharacterized protein YffS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,564,4440TC100.0% 223.3 / NA 61intergenic (+72/‑79)yffR/yffSCPZ‑55 prophage; uncharacterized protein YffR/CPZ‑55 prophage; uncharacterized protein YffS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/34);  total (27/34)

AACGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGA  >  NC_000913/2564313‑2564576
                                                                                                                                    |                                                                                                                                    
aaCGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCtt                                                                                                                                <  1:879260/138‑1 (MQ=255)
aaCGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCtt                                                                                                                                >  2:879260/1‑138 (MQ=255)
aaCGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTcc                                                                                                                              <  1:936434/140‑1 (MQ=255)
  cGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGt                                                                                                                            <  1:592177/140‑1 (MQ=255)
  cGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGt                                                                                                                            <  2:202062/140‑1 (MQ=255)
        aaaaaaCCACCGCTTTTTTTTTGGTAGTTGAGTTATTGAAAAACCTTAATGAATTCCAATAGCAACCACCAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCt                                                                                                                      <  2:202312/138‑1 (MQ=255)
        aGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCt                                                                                                                      <  1:817091/140‑1 (MQ=255)
        aGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCt                                                                                                                      <  1:647861/140‑1 (MQ=255)
            aacaacGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTAcc                                                                                                                  <  1:406449/140‑1 (MQ=255)
             acaacGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTAccc                                                                                                                 >  1:367791/1‑140 (MQ=255)
               aacGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTc                                                                                                               <  2:189646/140‑1 (MQ=255)
                 cGGCTTTGTTTTGGTAGGTAGAGAATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATGTGCGTCAATATTCCCTTCCGTTACGCTTACCCTcca                                                                                                             >  2:461669/1‑140 (MQ=255)
                         gTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATAc                                                                                                      <  2:71047/140‑1 (MQ=255)
                             tGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCgt                                                                                                  <  1:945247/140‑1 (MQ=255)
                                   gtagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAAc                                                                                            >  2:46641/1‑140 (MQ=255)
                                    tagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCACTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTATGCTTGCCCTCCACCAATACGCGTGTTAAcc                                                                                           >  2:819927/1‑140 (MQ=255)
                                      gagTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTc                                                                                         <  2:980178/140‑1 (MQ=255)
                                      gagTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTc                                                                                         <  2:35438/140‑1 (MQ=255)
                                      gagTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTc                                                                                         <  2:800155/140‑1 (MQ=255)
                                       agTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTcc                                                                                        >  2:542671/1‑140 (MQ=255)
                                       agTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTcc                                                                                        >  2:625460/1‑140 (MQ=255)
                                       agTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTcc                                                                                        >  2:617367/1‑140 (MQ=255)
                                       agTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTcc                                                                                        >  1:140870/1‑140 (MQ=255)
                                            ttGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATg                                                                                   <  2:553331/140‑1 (MQ=255)
                                            ttGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATg                                                                                   >  2:1090883/1‑140 (MQ=255)
                                              gAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACCCCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATTCTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGa                                                                                 <  1:461669/140‑1 (MQ=255)
                                              gAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGa                                                                                 <  2:367791/140‑1 (MQ=255)
                                                 aaaCGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTc                                                                              >  1:929237/1‑140 (MQ=255)
                                                 aaaCGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTc                                                                              >  1:512214/1‑140 (MQ=255)
                                                   aCGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTa                                                                            >  2:312982/1‑140 (MQ=255)
                                                       aaTTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATgaga                                                                       <  2:881341/140‑1 (MQ=255)
                                                          tgattgCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATg                                                                     <  1:153805/140‑1 (MQ=255)
                                                          tgattgCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATg                                                                     <  2:20926/140‑1 (MQ=255)
                                                           gattgCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATgg                                                                    <  2:546918/140‑1 (MQ=255)
                                                            attgCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATggg                                                                   >  1:1016694/1‑140 (MQ=255)
                                                               gCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCGATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGgaga                                                                >  2:1096447/1‑140 (MQ=255)
                                                               gCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGgaga                                                                <  1:46641/140‑1 (MQ=255)
                                                                       caacaaCTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATc                                                        <  1:690929/140‑1 (MQ=255)
                                                                               aaaCAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTa                                                <  1:166209/140‑1 (MQ=255)
                                                                                aacaacCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTAc                                               >  1:540738/1‑140 (MQ=255)
                                                                                        aTACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCaa                                       >  2:740961/1‑140 (MQ=255)
                                                                                             cccGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAAttt                                  >  2:83179/1‑140 (MQ=255)
                                                                                               cGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTACCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGaa                                            >  2:1170245/1‑128 (MQ=255)
                                                                                               cGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTACCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGaa                                            <  1:1170245/128‑1 (MQ=255)
                                                                                                  ccccTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGAc                             >  2:649850/1‑140 (MQ=255)
                                                                                                   cccTCGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATTGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACa                            >  2:991849/1‑140 (MQ=255)
                                                                                                       cGCCGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCt                        >  2:314113/1‑140 (MQ=255)
                                                                                                         ccGGGTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCa                      <  1:83179/140‑1 (MQ=255)
                                                                                                           gggTTTTTTGTGATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCaaa                    <  1:314113/140‑1 (MQ=255)
                                                                                                                      gATCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt         >  1:714847/1‑140 (MQ=255)
                                                                                                                       aTCTGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg        >  2:348168/1‑140 (MQ=255)
                                                                                                                          tGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGaa     >  2:209053/1‑140 (MQ=255)
                                                                                                                          tGCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGaa     >  2:18738/1‑140 (MQ=255)
                                                                                                                           gCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg    <  2:714847/140‑1 (MQ=255)
                                                                                                                           gCGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg    <  1:312982/140‑1 (MQ=255)
                                                                                                                            cGTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg   >  2:815930/1‑140 (MQ=255)
                                                                                                                             gTCAATTCTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa  <  1:1075077/140‑1 (MQ=255)
                                                                                                                             gTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa  <  2:263826/140‑1 (MQ=255)
                                                                                                                             gTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa  <  1:625460/140‑1 (MQ=255)
                                                                                                                             gTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa  <  1:542671/140‑1 (MQ=255)
                                                                                                                             gTCAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa  <  2:499539/140‑1 (MQ=255)
                                                                                                                               cAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGa      <  2:592337/134‑1 (MQ=255)
                                                                                                                               cAATACTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGa      >  1:592337/1‑134 (MQ=255)
                                                                                                                                    |                                                                                                                                    
AACGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGA  >  NC_000913/2564313‑2564576

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