Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,513,930 A→T 100% L413F (TTA→TTT USA300HOU_1379 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,513,9300AT100.0% 208.1 / ‑6.5 66L413F (TTA→TTTUSA300HOU_1379hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (28/37);  minor base G (0/1);  total (28/38)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACAA  >  CP000730/1513794‑1514066
                                                                                                                                        |                                                                                                                                        
ttttatttGTATGTATCGCTTTATTATTAGCGATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTg                                                                                                                                      <  2:399154/140‑1 (MQ=255)
 tttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAATTACCAAATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGt                                                                                                                                     <  2:294479/141‑1 (MQ=255)
    atttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTAtt                                                                                                                                  <  1:576936/141‑1 (MQ=255)
    atttGTATGTATCGCTTTATGAGTAGCCATTGGCATGGGGTTTGACTCTTGGATAAAGAGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCATCATCTATGTCACCCATAAGATGTTTTATTGTAtt                                                                                                                                  >  1:459688/1‑141 (MQ=255)
      ttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAg                                                                                                                                >  1:645663/1‑141 (MQ=255)
       tgtatgtaTCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAgg                                                                                                                               <  1:285436/141‑1 (MQ=255)
       tgtatgtaTCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGCACGCAACCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAgg                                                                                                                               <  1:399719/141‑1 (MQ=255)
                   tttATTAGTAGCCATTGGCATGGGGTTTGTCCCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCTTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTCATGcc                                                                                                                   <  2:277235/141‑1 (MQ=255)
                    ttattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                                                                  <  1:193950/141‑1 (MQ=255)
                      attaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatt                                                                                                                >  2:549753/1‑141 (MQ=255)
                       ttaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTCAAAAAACCAGCAGCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatta                                                                                                               >  2:519446/1‑141 (MQ=255)
                       ttaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatta                                                                                                               >  2:421568/1‑141 (MQ=255)
                           tagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTg                                                                                                           >  1:291787/1‑141 (MQ=255)
                           tagCCACTGGCATGGGCTTTCTCTCTCGGATAAAGCGTCCAAAAAACACTGCAATAAAAGTACCACATCGCGTAAAAAAACCCTTATCTATATCCCTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTg                                                                                                           <  1:633268/141‑1 (MQ=255)
                            agCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTgg                                                                                                          <  1:107619/141‑1 (MQ=255)
                              ccATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACCAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAt                                                                                                        <  1:325812/141‑1 (MQ=255)
                                agtGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                                                                      <  2:137563/139‑1 (MQ=255)
                                aTTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                                                                      <  2:137527/141‑1 (MQ=255)
                                   ggCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                                                                   >  1:224623/1‑141 (MQ=255)
                                      aTGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACtt                                                                                                >  2:635846/1‑141 (MQ=255)
                                        ggggTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTAt                                                                                              <  1:224107/141‑1 (MQ=255)
                                             ttGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTCTCACTTATCCttt                                                                                         >  2:664799/1‑140 (MQ=255)
                                                tctcTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAt                                                                                      <  2:16524/141‑1 (MQ=255)
                                                 ctctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAtt                                                                                     <  1:145051/141‑1 (MQ=255)
                                                     tGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTtata                                                                                 >  1:357929/1‑141 (MQ=255)
                                                        aTAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                                                                              >  2:445793/1‑141 (MQ=255)
                                                        aTAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                                                                              <  2:716351/141‑1 (MQ=255)
                                                         tAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                                                                             >  1:575237/1‑141 (MQ=255)
                                                          acaGCCTACAAAAAACACTGCAATAATAGTACTACATCTCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATGAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGa                                                                            <  1:655731/139‑1 (MQ=255)
                                                             gcatACAAAAAATACTGCAATAAAAATACCACATCGCGTAAAAAAACCAGCATCGATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                                                                         <  2:16355/138‑1 (MQ=255)
                                                             gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTGAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATAACTTATCCTTGTATTTATAATTGAAtt                                                                         >  1:281261/1‑141 (MQ=255)
                                                                  cAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAatat                                                                    <  1:630842/141‑1 (MQ=255)
                                                                     aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat                                                                 >  1:498486/1‑141 (MQ=255)
                                                                      aaaTGCTGCCCTAAAAGTAGCACATCCCGTAAAAAAACCAGCATCTATATCACTCATAATGTGTTTTATTGTATTAGGATTATTAATGCCATTCTTTGGATTGTCACTTATCCTTGTATTTATAATTGCATTAAttttatt                                                                >  2:71636/1‑140 (MQ=255)
                                                                          aCTGCAATAAAAGTAGCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATAtt                                                            >  1:382137/1‑141 (MQ=255)
                                                                           cTGCAATAAAAGTACCACACCGTAGAAAAAAACCAGCACCTATATCACTCATAATATGCTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTa                                                           <  2:567700/141‑1 (MQ=255)
                                                                           cTGCAATAAAAGTACAACAACGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTa                                                           <  2:363083/141‑1 (MQ=255)
                                                                            tGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaa                                                          >  2:659294/1‑141 (MQ=255)
                                                                             gCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaaa                                                         >  1:209225/1‑141 (MQ=255)
                                                                                        aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCCTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGc                                              >  1:602003/1‑141 (MQ=255)
                                                                                           acaTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATGATTAATCCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTaa                                           <  1:552800/141‑1 (MQ=255)
                                                                                             aTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTTCTacac                                         >  2:204986/1‑141 (MQ=255)
                                                                                             aTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAAc                                         <  2:291787/141‑1 (MQ=255)
                                                                                                 ctcaaaaaaaCCATTATATATCTCACTCCTTATATGTTTTATTGTATTAGGATTATTAATGCCATTTTTTGCATTATCCCTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaata                                     <  2:129372/138‑1 (MQ=255)
                                                                                                 cgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaata                                     <  2:435755/141‑1 (MQ=255)
                                                                                                  gTAAAAAAACCAGCATCTATATCACTCAAAATATGTTTTATTGTATTAGGATTAGTAATGCCATTATTTGGATTATCACGTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaataa                                    >  1:490322/1‑141 (MQ=255)
                                                                                                     aaaaaaGCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATAACATATTAAAGATCGTCGTGCTAAACaataatac                                 >  2:701114/1‑139 (MQ=255)
                                                                                                     aaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaattattc                                 >  2:445039/1‑139 (MQ=255)
                                                                                                       aaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCAc                               <  2:30093/141‑1 (MQ=255)
                                                                                                        aaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACt                              <  1:187220/141‑1 (MQ=255)
                                                                                                        aaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACt                              <  2:472399/141‑1 (MQ=255)
                                                                                                        aaaCCAGCACCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACt                              <  1:354088/141‑1 (MQ=255)
                                                                                                           ccAGCATTCATATCACTTATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                           <  2:182321/141‑1 (MQ=255)
                                                                                                           ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTGATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                           >  1:146686/1‑141 (MQ=255)
                                                                                                           ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTACACAATAATGCACTTaa                           <  1:421568/141‑1 (MQ=255)
                                                                                                           ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATAGTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                           <  1:519446/141‑1 (MQ=255)
                                                                                                           ccAGCATCCATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATTCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                           >  2:559610/1‑141 (MQ=255)
                                                                                                           ccAGCATATATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                           <  2:182324/141‑1 (MQ=255)
                                                                                                              gCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGt                        >  2:134612/1‑141 (MQ=255)
                                                                                                                aTCTATATCCTTCTTAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGttt                      <  2:194552/141‑1 (MQ=255)
                                                                                                                    ttatcatctataataTGTTTTTATGTATTAGGATTATTAATGCCATTATTGTGATCATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATGATGCACTTAAAGTTTTGaa                  <  2:158089/132‑1 (MQ=255)
                                                                                                                        cACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGa              <  2:759358/141‑1 (MQ=255)
                                                                                                                                atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAttt      >  2:325649/1‑141 (MQ=255)
                                                                                                                                  ttGTTTGATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAc    <  2:645663/140‑1 (MQ=255)
                                                                                                                                  atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAc    <  1:445793/141‑1 (MQ=255)
                                                                                                                                    gTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACaa  >  1:669044/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                        
ATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACAA  >  CP000730/1513794‑1514066

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: