Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,483,113 C→A S18* (TCG→TAG)  hrpA → ATP‑dependent RNA helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,483,1130CA100.0% 216.2 / NA 62S18* (TCG→TAG) hrpAATP‑dependent RNA helicase HrpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (33/29);  total (33/29)

ACGCATAGAAAGCCGCGATGTGGTACTCTATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTCGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCC  >  NC_000913/1482979‑1483244
                                                                                                                                      |                                                                                                                                   
acGCATAGAAAGCCGCGATGTGGTACTCTATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGa                                                                                                                                <  2:850699/140‑1 (MQ=255)
  gCATAGAAAGCCGCGATGTGGTACTCTATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGAtg                                                                                                                              <  1:1118857/140‑1 (MQ=255)
    aTAGAAAGCCGCGATGTGGTACTCTATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGAtgct                                                                                                                            <  1:1101732/140‑1 (MQ=255)
    aTAGAAAGCCGCGATGTGGTACTCTATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGAtgct                                                                                                                            <  2:797409/140‑1 (MQ=255)
            ccGCGATGTGGTACTCTATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACa                                                                                                                    >  1:700019/1‑140 (MQ=255)
                         ctctATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTtc                                                                                                       >  2:987109/1‑140 (MQ=255)
                            tataTCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCg                                                                                                    >  1:1010554/1‑140 (MQ=255)
                            tataTCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCg                                                                                                    >  1:844744/1‑140 (MQ=255)
                            tataTCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCg                                                                                                    >  1:701099/1‑140 (MQ=255)
                             atatCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGc                                                                                                   <  2:425027/140‑1 (MQ=255)
                              tatctatcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGcc                                                                                                  >  1:177175/1‑140 (MQ=255)
                                tctatcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGTGTTTTTCTCGCCGt                                                                                                <  1:51897/140‑1 (MQ=255)
                                 ctatcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTc                                                                                               >  2:587019/1‑140 (MQ=255)
                                 ctatcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTc                                                                                               >  1:738119/1‑140 (MQ=255)
                                 ctatcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTc                                                                                               >  1:1116331/1‑140 (MQ=255)
                                   atcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTg                                                                                             >  2:557822/1‑140 (MQ=255)
                                    tcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGc                                                                                            >  1:1203523/1‑140 (MQ=255)
                                    tcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGc                                                                                            >  1:1020725/1‑140 (MQ=255)
                                    tcATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGc                                                                                            >  1:831605/1‑140 (MQ=255)
                                     cATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCa                                                                                           <  2:738119/140‑1 (MQ=255)
                                     cATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCa                                                                                           <  2:277121/140‑1 (MQ=255)
                                     cATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCa                                                                                           >  2:515329/1‑140 (MQ=255)
                                     cATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTACGTGACAGACTGCGTTTTTCTCGCCGTCTGCa                                                                                           <  2:1051459/140‑1 (MQ=255)
                                      aTTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCAc                                                                                          >  1:1114297/1‑140 (MQ=255)
                                        ttAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACgg                                                                                        >  1:191067/1‑140 (MQ=255)
                                        ttAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACgg                                                                                        >  2:219426/1‑140 (MQ=255)
                                               aaaaTTAATCAGGCAGAGTACTGCCCACTAACGTTGTGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTgaag                                                                                 <  2:262120/140‑1 (MQ=255)
                                               aaaaTTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTgaag                                                                                 <  2:502332/140‑1 (MQ=255)
                                                aaaTTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTgaaga                                                                                <  2:177175/140‑1 (MQ=255)
                                                 aaTTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTgaagaa                                                                               <  2:1116331/140‑1 (MQ=255)
                                                        cAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTaaa                                                                        <  1:199030/140‑1 (MQ=255)
                                                             aGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATcc                                                                   <  2:1114297/140‑1 (MQ=255)
                                                                       ccACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCacaa                                                         >  2:1018069/1‑140 (MQ=255)
                                                                        cACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCacaac                                                        >  2:286578/1‑140 (MQ=255)
                                                                         aCTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCacaaca                                                       >  1:1023781/1‑140 (MQ=255)
                                                                         aCTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCacaaca                                                       >  2:874030/1‑140 (MQ=255)
                                                                         aCTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCacaaca                                                       >  2:992891/1‑140 (MQ=255)
                                                                         aCTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCacaaca                                                       >  1:389095/1‑140 (MQ=255)
                                                                            aaCGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGc                                                    <  1:1180656/140‑1 (MQ=255)
                                                                             aCGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGcc                                                   >  1:620919/1‑140 (MQ=255)
                                                                              cGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCa                                                  <  1:1018069/140‑1 (MQ=255)
                                                                              cGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCa                                                  <  2:149848/140‑1 (MQ=255)
                                                                              cGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCa                                                  <  1:419990/140‑1 (MQ=255)
                                                                                  aTGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCAtttt                                              >  1:1040029/1‑140 (MQ=255)
                                                                                        gAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGga                                        <  2:1020725/140‑1 (MQ=255)
                                                                                        gAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGga                                        <  2:844744/140‑1 (MQ=255)
                                                                                         aacaacAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGgag                                       >  2:785655/1‑140 (MQ=255)
                                                                                            aacaaAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATg                                    <  2:1023781/140‑1 (MQ=255)
                                                                                            aacaaAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATg                                    <  2:389095/140‑1 (MQ=255)
                                                                                                aaaaTTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGa                                <  1:587019/140‑1 (MQ=255)
                                                                                                 aaaTTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGaa                               <  1:515329/140‑1 (MQ=255)
                                                                                                             aCGGCCTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCa                   <  2:191067/140‑1 (MQ=255)
                                                                                                                  cTTGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAggcgg              <  1:785655/140‑1 (MQ=255)
                                                                                                                   ttGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAggcggc             >  1:725658/1‑140 (MQ=255)
                                                                                                                   ttGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAggcggc             >  1:85156/1‑140 (MQ=255)
                                                                                                                    tGCAGCAGCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCa            >  1:1101951/1‑140 (MQ=255)
                                                                                                                      cagcagCGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAgg          >  1:325936/1‑140 (MQ=255)
                                                                                                                        gcagcGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGACATGGCGAAAGAGATTGACCAGGCGGCAGGGa        <  2:1010554/140‑1 (MQ=255)
                                                                                                                           gcGGCTGGATTAGCTGATGCTGCTTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAg     >  1:1065852/1‑140 (MQ=255)
                                                                                                                           gcGGCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAg     >  1:1018389/1‑140 (MQ=255)
                                                                                                                             ggCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTc   <  2:1203523/140‑1 (MQ=255)
                                                                                                                              gCTGGATTAGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTcc  >  2:529553/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                   
ACGCATAGAAAGCCGCGATGTGGTACTCTATATCTATCATTTAAAAGAAAATTAATCAGGCAGACTACTGCCCACTAACGTTATGACAGAACAACAAAAATTGACCTTTACGGCCTTGCAGCAGCGGCTGGATTCGCTGATGCTGCGTGACAGACTGCGTTTTTCTCGCCGTCTGCACGGCGTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCC  >  NC_000913/1482979‑1483244

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