Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,359,120 A→G G304G (GGT→GGC cadB ← lysine:cadaverine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,359,1200AG98.3% 190.6 / ‑4.7 58G304G (GGT→GGCcadBlysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (29/28);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GTCAGATGCTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCAC  >  NC_000913/4358985‑4359250
                                                                                                                                       |                                                                                                                                  
gTCAGATGCTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCcgt                                                                                                                                >  1:557780/1‑140 (MQ=255)
 tCAGATGCTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCcgtc                                                                                                                               >  2:46934/1‑140 (MQ=255)
   aGATGCTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCcgtcgt                                                                                                                             <  2:72495/140‑1 (MQ=255)
     aTGCTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTa                                                                                                                           >  2:604626/1‑140 (MQ=255)
     aTGCTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTa                                                                                                                           >  2:505540/1‑140 (MQ=255)
        cTTTACCACCGGCAGGGTTCATCAGAGTGGTAAGGATTATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCg                                                                                                                        <  1:298618/140‑1 (MQ=255)
        cTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCg                                                                                                                        <  1:505540/140‑1 (MQ=255)
          ttACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTGGCTGTCGATTTCCCCAAAAACTTTCAGGAAATTGCCGTCGTTAGCGGc                                                                                                                      >  2:224029/1‑140 (MQ=255)
          ttACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTCGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGc                                                                                                                      >  1:609622/1‑140 (MQ=255)
                  ggCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAAATTTCGGGAAGTTGCCGTGGTTAGCGGCACGTacac                                                                                                              >  2:280052/1‑140 (MQ=255)
                  ggCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTAATTTTCACCGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTacac                                                                                                              <  2:205371/140‑1 (MQ=255)
                   gCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACAcc                                                                                                             >  1:587825/1‑140 (MQ=255)
                       agTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACAcctgcc                                                                                                         >  1:86392/1‑140 (MQ=255)
                          tcatcaGAGTGATAAGGTTCATCAGGGCAGTGATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGc                                                                                                     <  2:238659/140‑1 (MQ=255)
                               agagTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTAc                                                                                                 <  2:23621/140‑1 (MQ=255)
                                   tGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTGGGGAAGTTGCCGGCGTAAGCGGCACCTACAGCTGCCTGGCCTACCAAc                                                                                             >  2:183767/1‑140 (MQ=255)
                                   tGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAAc                                                                                             >  2:11186/1‑140 (MQ=255)
                                   tGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAAc                                                                                             <  2:86392/140‑1 (MQ=255)
                                       aaGGATAATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAAcatca                                                                                         <  1:403839/140‑1 (MQ=255)
                                           atcatcAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACAGCTGCCTGGCTTACCAACATCATCCa                                                                                     >  1:168734/1‑140 (MQ=255)
                                            tcatcaGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAg                                                                                    <  2:496546/140‑1 (MQ=255)
                                            ccatcaGGGCAGTCATTTTCACTGCAGTCAGCAGCAGAACTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACATTCGGAAAGAGGACGTCGTTAGCGGCACGTCCACCTGCCTGGCCTACCAAcatcaacaac                                                                                    <  2:609622/139‑6 (MQ=255)
                                             catcaGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTAGACCAGCTTGGCCTGCCAACATCATCCAgg                                                                                   >  1:192489/1‑140 (MQ=255)
                                             catcaGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCAGCACGTAGACCAGCTTGGCCTACCAACATCATCCAgg                                                                                   >  1:221637/1‑140 (MQ=255)
                                               tcaGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCCACCAACATCATCCAGGAg                                                                                 >  2:497452/1‑140 (MQ=255)
                                                caGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGc                                                                                <  1:497452/140‑1 (MQ=255)
                                                   ggCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCa                                                                             <  1:107696/140‑1 (MQ=255)
                                                   ggCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGGACACCTGGCTGCCCTACCAACATCATCCAGGAGCCCa                                                                             >  1:293942/1‑140 (MQ=255)
                                                    gCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCag                                                                            >  1:47676/1‑140 (MQ=255)
                                                     cAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCaga                                                                           <  1:5720/140‑1 (MQ=255)
                                                        tCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGaa                                                                        >  2:609663/1‑140 (MQ=255)
                                                        tCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGaa                                                                        >  2:129067/1‑140 (MQ=255)
                                                        tCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGaa                                                                        >  1:141113/1‑140 (MQ=255)
                                                        tCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGaa                                                                        >  1:575361/1‑140 (MQ=255)
                                                             ttCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAg                                                                   >  1:364017/1‑140 (MQ=255)
                                                             ttCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACCTACACCTGACTGGCCTCCCAACATCTTCCACGAGCCCAGAGAAGTCAg                                                                   >  2:221517/1‑140 (MQ=255)
                                                              tCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAgg                                                                  <  1:112662/140‑1 (MQ=255)
                                                                    cagccagcAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTCCACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCa                                                            <  2:135218/140‑1 (MQ=255)
                                                                    cagccagcAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCa                                                            <  2:364017/140‑1 (MQ=255)
                                                                    cagccagcAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCa                                                            <  2:587825/140‑1 (MQ=255)
                                                                     agccagcAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCaa                                                           <  1:609663/140‑1 (MQ=255)
                                                                         agcagcagACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTAGCCTACCAACATCATCGAGGAGCCCGGAGAGGTCAGGCACGCAAAGGc                                                       >  2:198858/1‑140 (MQ=255)
                                                                            agcagACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGt                                                    <  2:47676/140‑1 (MQ=255)
                                                                                   cTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCa                                             >  1:356467/1‑140 (MQ=255)
                                                                                           ggAATACCGTTGCTGTCAACTTCACCATAAACTTTCGGGCAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAg                                     <  2:171136/140‑1 (MQ=255)
                                                                                                aCCGTTGCTGTCGACATCACCATAAGCTTCCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAgcggcg                                <  1:285520/140‑1 (MQ=255)
                                                                                                  cGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCa                              >  2:547639/1‑140 (MQ=255)
                                                                                                     tGCTGTCGACTTCACCATAAACATTCGGGAAGTTGACGTTGTTAGCGGTACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGcc                           <  2:221637/140‑1 (MQ=255)
                                                                                                     tGCTGTCGACTTCACAATAAACTTTCGGGAAGTTGCCTTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGcc                           <  1:224029/140‑1 (MQ=255)
                                                                                                           cGACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTa                     <  2:612534/140‑1 (MQ=255)
                                                                                                            gACTTCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTAc                    >  2:151845/1‑140 (MQ=255)
                                                                                                               ttCACCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTGCACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCgg                 >  2:130490/1‑139 (MQ=255)
                                                                                                                  aCCATAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGa              <  2:557780/140‑1 (MQ=255)
                                                                                                                   ccATAAACTTTCGGGTAGTTGCCGTCGTTTGCGGCACTTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGAt             <  1:198858/140‑1 (MQ=255)
                                                                                                                   ccATAAACTTTCGGGCAGTTGCCGTCGTTTGCGGCACATACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGAt             <  1:183767/140‑1 (MQ=255)
                                                                                                                   ccATAAACTTTCGGGCAGTTGCAGTCGTTAGCGGCACTTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGAt             <  1:280052/140‑1 (MQ=255)
                                                                                                                     aTAAACTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAg           <  1:604626/140‑1 (MQ=255)
                                                                                                                       aaaCTTTCGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAACGCGGTGAAGGCAGAAATCAGCGGCGTAGCCCAGTTACCGAGGATAGtt         >  2:261725/1‑140 (MQ=255)
                                                                                                                              cGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGTATGCAGAACCCAGCCGCGCAGCCCAATTACGGAGGATAGGTGAAGCAc  >  2:306097/1‑140 (MQ=255)
                                                                                                                              cGGGAAGTTGCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCAc  <  1:46934/140‑1 (MQ=255)
                                                                                                                              cGGGAAGTTGCCGTCCTTAGCGCCAAGTATACCTGCCGGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCAc  <  1:221517/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
GTCAGATGCTTTACCACCGGCAGAGTTCATCAGAGTGATAAGGATCATCAGGGCAGTCATTTTCACTGCAGCCAGCAGCAGACCTTTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCAC  >  NC_000913/4358985‑4359250

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