Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,275,815 T→C E5E (GAA→GAG narL ← response regulator in two‑component regulatory system with NarX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,275,8150TC100.0% 234.8 / NA 67E5E (GAA→GAGnarLresponse regulator in two‑component regulatory system with NarX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/35);  total (32/35)

ACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275682‑1275950
                                                                                                                                     |                                                                                                                                       
aCTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTg                                                                                                                                    >  1:2048300/1‑139 (MQ=255)
  tCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGAt                                                                                                                                  >  1:636960/1‑139 (MQ=255)
    cGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTa                                                                                                                                >  1:275734/1‑139 (MQ=255)
       cAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTc                                                                                                                             <  2:636960/139‑1 (MQ=255)
              aTACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTa                                                                                                                      <  1:44796/139‑1 (MQ=255)
                    tGTTCGCCATCACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCt                                                                                                                >  1:1066906/1‑139 (MQ=255)
                     gTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCtt                                                                                                               >  2:207214/1‑139 (MQ=255)
                      ttCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATgg                                                                                                                         >  2:875693/1‑128 (MQ=255)
                      ttCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATgg                                                                                                                         <  1:875693/128‑1 (MQ=255)
                         gCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGAc                                                                                                           <  1:259140/139‑1 (MQ=255)
                                cTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       >  1:475537/1‑120 (MQ=255)
                                cTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  2:475537/120‑1 (MQ=255)
                                cTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCtgtg                                                                                                    >  1:1734236/1‑139 (MQ=255)
                                           aaCCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         <  2:1734236/139‑1 (MQ=255)
                                                  gTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATa                                                                                  <  1:207214/139‑1 (MQ=255)
                                                      tatCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAgg                                                                              <  2:1066906/139‑1 (MQ=255)
                                                           ggTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGacc                                                                         >  2:1999553/1‑139 (MQ=255)
                                                             tGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGaccac                                                                       <  2:275734/139‑1 (MQ=255)
                                                             tGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGaccac                                                                       >  2:1989980/1‑139 (MQ=255)
                                                             tGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGa                                                                           <  1:964903/135‑1 (MQ=255)
                                                             tGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGa                                                                           >  2:964903/1‑135 (MQ=255)
                                                               ccATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGAccacca                                                                     >  2:2017350/1‑139 (MQ=255)
                                                                  tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       >  2:1303338/1‑86 (MQ=255)
                                                                  tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  1:1303338/86‑1 (MQ=255)
                                                                  tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACtt                                                                  <  1:912038/139‑1 (MQ=255)
                                                                      gATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGaa                                                                                    <  2:81164/117‑1 (MQ=255)
                                                                      gATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGaa                                                                                    >  1:81164/1‑117 (MQ=255)
                                                                       aTAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGCCGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTg                                                             >  1:367720/1‑139 (MQ=255)
                                                                       aTAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATaa                                                                                 <  2:1471703/119‑1 (MQ=255)
                                                                       aTAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATaa                                                                                 >  1:1471703/1‑119 (MQ=255)
                                                                       aTAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTg                                                             >  2:1787139/1‑139 (MQ=255)
                                                                       aTAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTg                                                             >  2:1933935/1‑139 (MQ=255)
                                                                        tAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGc                                                            <  2:1792643/139‑1 (MQ=255)
                                                                          aGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCa                                                          <  1:2017350/139‑1 (MQ=255)
                                                                            cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATa                                                                                  <  2:664322/113‑1 (MQ=255)
                                                                            cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATa                                                                                  >  1:664322/1‑113 (MQ=255)
                                                                            cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCAcc                                                        >  1:957527/1‑139 (MQ=255)
                                                                            cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCAcc                                                        <  2:2048300/139‑1 (MQ=255)
                                                                                 ttACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            >  2:540261/1‑66 (MQ=255)
                                                                                 ttACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            <  1:540261/66‑1 (MQ=255)
                                                                                 ttACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGAtt                                                   >  1:90372/1‑139 (MQ=255)
                                                                                         gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTc                                                                                                                                       <  1:339611/47‑1 (MQ=255)
                                                                                         gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTg                                                             <  2:405639/121‑1 (MQ=255)
                                                                                         gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTg                                                             >  1:405639/1‑121 (MQ=255)
                                                                                         gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGAcgg                                           >  2:1921209/1‑139 (MQ=255)
                                                                                         gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGAcgg                                           >  1:634109/1‑139 (MQ=255)
                                                                                         gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGCGACCACCACTTCGGTGCCACCTGATTCACGAcgg                                           <  1:612237/139‑1 (MQ=255)
                                                                                                    cGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCa                                >  2:525685/1‑139 (MQ=255)
                                                                                                     gggTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCaa                               <  1:1787139/139‑1 (MQ=255)
                                                                                                      ggTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACg                                                                                                          <  1:1118191/63‑1 (MQ=255)
                                                                                                      ggTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACg                                                                                                          >  2:1118191/1‑63 (MQ=255)
                                                                                                      ggTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAAt                              >  2:1177181/1‑139 (MQ=255)
                                                                                                       gTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTGTGTGAAAGTTTTTTCGTGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATc                             <  2:367720/139‑1 (MQ=255)
                                                                                                       gTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATc                             <  1:1177181/139‑1 (MQ=255)
                                                                                                       gTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATc                             <  2:90372/139‑1 (MQ=255)
                                                                                                                   cagcagGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGGCGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACt                 <  1:1496002/139‑1 (MQ=255)
                                                                                                                   cagcagGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACt                 <  1:525685/139‑1 (MQ=255)
                                                                                                                        ggATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATaa                                                                                 <  1:1086540/70‑1 (MQ=255)
                                                                                                                        ggATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATaa                                                                                 >  2:1086540/1‑70 (MQ=255)
                                                                                                                          aTAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCa          <  1:1933935/139‑1 (MQ=255)
                                                                                                                          aTAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCa          <  2:634109/139‑1 (MQ=255)
                                                                                                                              tagCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         <  1:855383/56‑1 (MQ=255)
                                                                                                                              tagCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         >  2:855383/1‑56 (MQ=255)
                                                                                                                              tagCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGAt      >  2:1303865/1‑139 (MQ=255)
                                                                                                                                gCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcg    >  2:1255606/1‑139 (MQ=255)
                                                                                                                                  cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  1:1999553/139‑1 (MQ=255)
                                                                                                                                  cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  1:761884/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
ACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275682‑1275950

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