Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 969,352 G→A *329* (TAG→TAA lpxK → lipid A 4'kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913969,3520GA100.0% 212.5 / NA 66*329* (TAG→TAAlpxKlipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/34);  total (32/34)

GACTGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTGGAGGA  >  NC_000913/969217‑969486
                                                                                                                                       |                                                                                                                                      
gACTGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTa                                                                                                                                     <  2:560733/139‑1 (MQ=255)
gACTGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTa                                                                                                                                     >  2:806134/1‑139 (MQ=255)
     aaaaaGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTAcgccg                                                                                                                                >  1:731453/1‑139 (MQ=255)
     aaaaaGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTAcgccg                                                                                                                                >  1:1730246/1‑139 (MQ=255)
     aaaaaGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTAcgccg                                                                                                                                >  1:297547/1‑139 (MQ=255)
          gATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCa                                                                                                                           >  1:1292011/1‑139 (MQ=255)
           aTGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAg                                                                                                                          <  1:314213/139‑1 (MQ=255)
            tGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCTGCGGCAGc                                                                                                                         >  2:597429/1‑139 (MQ=255)
            tGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGc                                                                                                                         >  2:1098644/1‑139 (MQ=255)
             gCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGACACTGCTTACGCCACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCg                                                                                                                        <  1:1269679/139‑1 (MQ=255)
                       gCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGAt                                                                                                              <  1:806134/139‑1 (MQ=255)
                        ccGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                             <  1:1292862/139‑1 (MQ=255)
                           ggCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAg                                                                                                          >  1:1745326/1‑139 (MQ=255)
                                           aaTTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGc                                                                                          <  1:1432973/139‑1 (MQ=255)
                                           aaTTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGc                                                                                          <  2:297547/139‑1 (MQ=255)
                                            aTTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGcc                                                                                         <  2:731453/139‑1 (MQ=255)
                                            aTTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGcc                                                                                         <  2:885103/139‑1 (MQ=255)
                                             ttGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCg                                                                                        <  2:1745326/139‑1 (MQ=255)
                                               ggtggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCTGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCa                                                                                      <  1:597429/139‑1 (MQ=255)
                                                gtggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGAt                                                                                                                  <  1:928806/110‑1 (MQ=255)
                                                gtggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGAt                                                                                                                  >  2:928806/1‑110 (MQ=255)
                                                  ggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACct                                                                                   >  1:150826/1‑139 (MQ=255)
                                                  ggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACct                                                                                   <  2:808450/139‑1 (MQ=255)
                                                   gtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACctc                                                                                  <  1:1098644/139‑1 (MQ=255)
                                                     aTTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCtctc                                                                                >  1:726458/1‑139 (MQ=255)
                                                              tAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGa                                                                       >  2:371985/1‑139 (MQ=255)
                                                               aGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGAt                                                                      >  1:1545230/1‑139 (MQ=255)
                                                                    cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                                 >  1:1613370/1‑139 (MQ=255)
                                                                    cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                                 <  2:1292011/139‑1 (MQ=255)
                                                                     acaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACct                                                                                   >  1:1088992/1‑120 (MQ=255)
                                                                     acaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACct                                                                                   <  2:1088992/120‑1 (MQ=255)
                                                                        gCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAAt                                                             >  1:1459625/1‑139 (MQ=255)
                                                                        gCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAAt                                                             <  2:515031/139‑1 (MQ=255)
                                                                        gCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAAt                                                             >  2:1261956/1‑139 (MQ=255)
                                                                         cTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGc                                                                                          <  1:141675/109‑1 (MQ=255)
                                                                         cTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGc                                                                                          >  2:141675/1‑109 (MQ=255)
                                                                         cTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATc                                                            >  2:743616/1‑139 (MQ=255)
                                                                            tCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATG‑‑GTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGAt                                                                      <  1:1533835/124‑1 (MQ=255)
                                                                            tCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATG‑‑GTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGAt                                                                      >  2:1533835/1‑124 (MQ=255)
                                                                             cAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCa                                                        >  2:1685242/1‑139 (MQ=255)
                                                                             cAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCa                                                        <  2:150826/139‑1 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTc                                                                                                                     <  1:279037/72‑1 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTc                                                                                                                     >  2:279037/1‑72 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATg                                                                     >  1:1720431/1‑120 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATg                                                                     <  2:1720431/120‑1 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGc                                                  >  2:595267/1‑139 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCGTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGc                                                  >  1:997369/1‑139 (MQ=255)
                                                                                           gCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAg                                          <  2:1545230/139‑1 (MQ=255)
                                                                                             gAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGc                                        <  1:1685242/139‑1 (MQ=255)
                                                                                                  tGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGtt                                   >  2:1665375/1‑139 (MQ=255)
                                                                                                       aCGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTaaaca                              <  1:371985/139‑1 (MQ=255)
                                                                                                             cTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGtt                                                                                                                      <  1:926162/45‑1 (MQ=255)
                                                                                                             cTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGtt                                                                                                                      >  2:926162/1‑45 (MQ=255)
                                                                                                              tAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCccg                       >  1:1693061/1‑139 (MQ=255)
                                                                                                              tAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCccg                       <  1:475716/139‑1 (MQ=255)
                                                                                                                     gCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGc                                        >  2:1714896/1‑115 (MQ=255)
                                                                                                                     gCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGc                                        <  1:1714896/115‑1 (MQ=255)
                                                                                                                       tGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGca                                              <  1:92056/107‑1 (MQ=255)
                                                                                                                       tGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGca                                              >  2:92056/1‑107 (MQ=255)
                                                                                                                              tGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCg                                                                                              <  2:644052/52‑1 (MQ=255)
                                                                                                                              tGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCg                                                                                              >  1:644052/1‑52 (MQ=255)
                                                                                                                              tGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTg       <  1:1665375/139‑1 (MQ=255)
                                                                                                                              tGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTg       <  2:1693061/139‑1 (MQ=255)
                                                                                                                                 cAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTggag    <  2:1459625/139‑1 (MQ=255)
                                                                                                                                   aCTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTggagga  <  2:1613370/139‑1 (MQ=255)
                                                                                                                                   aCTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTggagga  >  1:111330/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GACTGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTGGAGGA  >  NC_000913/969217‑969486

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: