Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,208,517 C→T P129L (CCT→CTT) 
*201* (TAG→TAA
tfaP →
tfaE ←
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,5170CT100.0% 211.3 / NA 68P129L (CCT→CTT) 
*201* (TAG→TAA
tfaP
tfaE
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/30);  total (38/30)

ATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTAC  >  NC_000913/1208382‑1208638
                                                                                                                                       |                                                                                                                         
aTAGCAGACAGTGGCAGGGTAAGGCTACATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTAc                                                                                                                        <  2:425519/139‑1 (MQ=255)
           tGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                 >  1:558316/1‑135 (MQ=255)
           tGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                 <  2:558316/135‑1 (MQ=255)
           tGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAg                                                                                                             >  1:791927/1‑139 (MQ=255)
             gCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGACTTACGCCTCCGGCAGtt                                                                                                           >  2:988216/1‑139 (MQ=255)
                 ggTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttaggc                                                                                                       <  2:1026089/139‑5 (MQ=255)
                     aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGa                                                                                                  <  1:180738/139‑1 (MQ=255)
                      ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGAt                                                                                                 <  2:737444/139‑1 (MQ=255)
                             ttGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATcc                                                                                          >  2:875960/1‑139 (MQ=255)
                                 gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                      <  1:875960/139‑1 (MQ=255)
                                 gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                      <  1:488599/139‑1 (MQ=255)
                                   aGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                    <  2:791927/139‑1 (MQ=255)
                                     aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                  >  2:384997/1‑139 (MQ=255)
                                               gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                        >  2:560084/1‑139 (MQ=255)
                                               gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGAGATTGAATTGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                        >  1:930492/1‑139 (MQ=17)
                                                 tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGt                                                                                                            <  2:284673/102‑1 (MQ=255)
                                                 tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGt                                                                                                            >  1:284673/1‑102 (MQ=255)
                                                 tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                  <  2:1080665/127‑1 (MQ=255)
                                                 tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                  >  1:1080665/1‑127 (MQ=255)
                                                 tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGt                                                                      >  1:1426149/1‑139 (MQ=255)
                                                 tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGt                                                                      >  2:1384260/1‑139 (MQ=255)
                                                       gAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGt                                                                >  2:664911/1‑139 (MQ=255)
                                                               ataATTTGTGGCTGGATTACCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTg                                                                                 >  2:1060167/1‑114 (MQ=255)
                                                               ataATTTGTGGCTGGATTACCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTg                                                                                 <  1:1060167/114‑1 (MQ=255)
                                                                  aTTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                          >  1:1448355/1‑118 (MQ=255)
                                                                  aTTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                          <  2:1448355/118‑1 (MQ=255)
                                                                  aTTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCaa                                                     <  1:651652/139‑1 (MQ=255)
                                                                   tttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt                                                    >  1:1436816/1‑139 (MQ=255)
                                                                   tttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt                                                    >  1:747557/1‑139 (MQ=255)
                                                                     tgtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATc                                                                                           >  2:1367189/1‑98 (MQ=255)
                                                                     tgtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATc                                                                                           <  1:1367189/98‑1 (MQ=255)
                                                                      gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                      <  1:997210/102‑1 (MQ=255)
                                                                      gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                      >  2:997210/1‑102 (MQ=255)
                                                                       tgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTaa                                                <  1:988216/139‑1 (MQ=255)
                                                                        ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATATATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                               >  2:738043/1‑139 (MQ=255)
                                                                             gATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                   <  1:1371366/130‑1 (MQ=255)
                                                                             gATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                   >  2:1371366/1‑130 (MQ=255)
                                                                              aTTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                         <  2:747557/139‑1 (MQ=255)
                                                                              aTTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                         <  2:1436816/139‑1 (MQ=255)
                                                                               ataagagacaGCACAGGAACTGGTTGATACATCCAGTGCGCCAGATATTGACTGGGTTACGCCTCCGGCAGGTAGGCCAGATGACATCCCGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                        <  1:300815/130‑1 (MQ=14)
                                                                                      gAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTa                                 >  1:888341/1‑139 (MQ=255)
                                                                                       aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTaa                                >  1:417867/1‑139 (MQ=255)
                                                                                         gCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGt                              >  1:756121/1‑139 (MQ=255)
                                                                                            cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                 <  1:1409320/54‑1 (MQ=255)
                                                                                            cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                 >  2:1409320/1‑54 (MQ=255)
                                                                                            cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttagg                                                                                                        <  1:1165799/63‑4 (MQ=18)
                                                                                            cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttagg                                                                                                        >  2:1165799/1‑60 (MQ=18)
                                                                                            cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTg                           >  1:250513/1‑139 (MQ=255)
                                                                                                   tGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGttt                    >  1:903697/1‑139 (MQ=255)
                                                                                                     gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcgc                                                                                     >  2:713058/1‑72 (MQ=14)
                                                                                                     gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcgc                                                                                     <  1:713058/72‑1 (MQ=14)
                                                                                                     gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGttctgtctcttat                  >  2:300815/1‑127 (MQ=14)
                                                                                                        gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc               <  2:1181586/139‑1 (MQ=255)
                                                                                                         aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTa                                                 <  1:586863/104‑1 (MQ=255)
                                                                                                         aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTa                                                 >  2:586863/1‑104 (MQ=255)
                                                                                                         aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcct              >  2:296508/1‑139 (MQ=255)
                                                                                                             aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt          <  1:384997/139‑1 (MQ=255)
                                                                                                             aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt          <  1:738043/139‑1 (MQ=255)
                                                                                                             aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt          >  1:50264/1‑139 (MQ=255)
                                                                                                               ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                          <  2:74128/105‑1 (MQ=255)
                                                                                                               ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                          >  1:74128/1‑105 (MQ=255)
                                                                                                               ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa        >  2:1057043/1‑139 (MQ=255)
                                                                                                                  gTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGtt     <  1:664911/139‑1 (MQ=255)
                                                                                                                   tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                        >  2:1176085/1‑87 (MQ=11)
                                                                                                                   tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                        <  1:1176085/87‑1 (MQ=11)
                                                                                                                   tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt    >  1:1235379/1‑139 (MQ=255)
                                                                                                                   tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCAGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt    >  2:1406552/1‑139 (MQ=255)
                                                                                                                     cgcCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTAc  >  2:644365/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                         
ATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTAC  >  NC_000913/1208382‑1208638

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: