Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,557 C→A 18.4% G407* (GGA→TGA)  pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,372,5570CA18.4% 122.9 / 16.1 62G407* (GGA→TGA) pbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base A (4/8);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGG  >  USA300TCH1516_ALE/2372419‑2372694
                                                                                                                                          |                                                                                                                                         
cAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCca                                                                                                                                          >  1:598912/1‑140 (MQ=255)
cAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAAGATGATGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCca                                                                                                                                          >  1:224888/1‑140 (MQ=255)
  gttgTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCcaca                                                                                                                                        <  2:30525/137‑1 (MQ=255)
      ttttttGTGTATAAAATTAACAATCATCCCATCCGAATTACAGAATGAATAGTACAAAAAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCa                                                                                                                                    <  2:354138/140‑1 (MQ=255)
              gtCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATCCCAGt                                                                                                                            >  1:321949/1‑140 (MQ=255)
              gtCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGt                                                                                                                            <  2:528410/140‑1 (MQ=255)
               tCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGtt                                                                                                                           <  1:44762/140‑1 (MQ=255)
               tCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTATTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGAAATAAAAGtt                                                                                                                           >  1:65/1‑140 (MQ=255)
                cTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTg                                                                                                                          <  2:159212/140‑1 (MQ=255)
                  ggATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCa                                                                                                                        <  1:460121/140‑1 (MQ=255)
                        aaCCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGaa                                                                                                                  >  1:741161/1‑140 (MQ=255)
                            aTGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATACGCATTGTAATCGGATAGAAAATAAATCCACACGCAATCCCAGTTGCAATTGAATATg                                                                                                              >  2:46461/1‑140 (MQ=255)
                                 cAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATCAATCAACACGACATATCAGTTGCAATTGAATATGAGaaa                                                                                                         >  2:475446/1‑140 (MQ=255)
                                  aaaTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAg                                                                                                        <  2:62490/140‑1 (MQ=255)
                                  aaaTATGAAGTACAGAATGAATAGTACCAATAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAg                                                                                                        <  2:464484/140‑1 (MQ=255)
                                   aaTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAgg                                                                                                       <  2:414995/140‑1 (MQ=255)
                                    atatGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGc                                                                                                      >  1:730846/1‑140 (MQ=255)
                                           gTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATaa                                                                                               <  1:190580/140‑1 (MQ=255)
                                               aGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCACTTGAATATGATAAAGGCATCATAATGAt                                                                                           >  1:683994/1‑140 (MQ=255)
                                                 aatgaatAGTACCGTTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTg                                                                                         <  2:646422/140‑1 (MQ=255)
                                                 aaagaatTGTACCATCAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGATATAAGCATTGTAATTGGATAGAAAATAAAGCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTg                                                                                         <  2:764733/137‑1 (MQ=255)
                                                         gTACCATTAATTCATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAg                                                                                 <  2:511514/140‑1 (MQ=255)
                                                          tACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGACATCATAATGATGGTAATGAAgg                                                                                >  2:198385/1‑140 (MQ=255)
                                                             cATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTg                                                                             >  2:124169/1‑140 (MQ=255)
                                                                   tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAATGCTGGAACtc                                                                       >  2:210898/1‑139 (MQ=255)
                                                                   tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTg                                                                       >  1:519518/1‑140 (MQ=255)
                                                                   tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTg                                                                       >  1:683455/1‑140 (MQ=255)
                                                                   tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAAGTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTg                                                                       >  2:204499/1‑140 (MQ=255)
                                                                        agatTATCGGTTGTACTTCTTTATGTTTTTTAGATATAAGCATTGTAATTGGATAGAAAATAAATCTACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACt                                                                  <  2:310644/138‑1 (MQ=255)
                                                                        acatTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACt                                                                  <  1:356972/140‑1 (MQ=255)
                                                                          atTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTc                                                                <  1:547577/140‑1 (MQ=255)
                                                                                ggATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACt                                                          >  1:671216/1‑140 (MQ=255)
                                                                                ggATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTCTAATGAATGCCGGCACTTCAACTTCAAAACt                                                          >  2:198193/1‑140 (MQ=255)
                                                                                  aTGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACgtt                                                        >  1:347983/1‑140 (MQ=255)
                                                                                     tACTTCTTTATATTTTTGAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTcc                                                     <  2:311669/140‑1 (MQ=255)
                                                                                       cttcttTATGTTTTTTAGAGATAAGCATTGTAATGGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAg                                                   <  2:598912/140‑1 (MQ=255)
                                                                                         tcttTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGAACTGAAACTTCAAAACTTTTCCAGtt                                                 >  1:735722/1‑140 (MQ=255)
                                                                                          cttTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTATTGAAGGCTGGCACTGCAACTTCTAAACTTTTCCAGttt                                                >  2:277274/1‑140 (MQ=255)
                                                                                             tATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATACACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATAATTGTAATGAAGGCTGGCAATGGAACTTCAAAAATTGTCCAgatgagt                                             >  2:587403/1‑134 (MQ=255)
                                                                                                 ttttttAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAg                                         >  1:302534/1‑140 (MQ=255)
                                                                                                  tttttAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGc                                        >  1:441960/1‑140 (MQ=255)
                                                                                                   ttttAGAGATAAGCATTGTAATGGGATAGAAAATATATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCa                                       <  2:689861/140‑1 (MQ=255)
                                                                                                      tatagaTAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg                                    <  1:478135/137‑1 (MQ=255)
                                                                                                      tAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAACGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg                                    <  2:741161/140‑1 (MQ=255)
                                                                                                      tAGAGATAAGCATAGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg                                    <  1:216646/140‑1 (MQ=255)
                                                                                                        gagaTAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGtt                                  <  2:671216/140‑1 (MQ=255)
                                                                                                              aGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGCTTGTAATGAAGGCTGGCACTGCCACTTCAAAACTCTACCAGCTTATCTCTGCAAAGTCTGCAGc                            >  1:711925/1‑140 (MQ=255)
                                                                                                              aGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGc                            <  1:501710/140‑1 (MQ=255)
                                                                                                               gCATTGTAATTGTATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGcc                           <  1:157758/140‑1 (MQ=255)
                                                                                                               gCATTGTAATTGGATAGAAATTGTTTGCGCACACAACACCAGATGCAGAAAAAGACTATAAAGGCATTATTGTGATGGGAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGcc                           <  1:181649/140‑1 (MQ=255)
                                                                                                               gCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGcc                           <  1:531799/140‑1 (MQ=255)
                                                                                                                cATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCa                          <  2:135866/140‑1 (MQ=255)
                                                                                                                  ttGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCACAATGAGTGGAATGACGGCTGGCACTGCAACTTCAAAACTCTCCCAGCTTATCTCAGCAAAGTTCGCAGCCAtt                        >  2:359791/1‑140 (MQ=255)
                                                                                                                  ttGTAATTGGATAGAAAATAAAGCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCAtt                        >  2:492470/1‑140 (MQ=255)
                                                                                                                    gTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTaa                      >  2:222642/1‑140 (MQ=255)
                                                                                                                       ttatgATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCACCTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGAAGCCATtactac                   >  2:536160/5‑140 (MQ=255)
                                                                                                                       aTTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATac                   >  1:497077/1‑140 (MQ=255)
                                                                                                                          ggATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACAcc                <  1:188938/140‑1 (MQ=255)
                                                                                                                             tAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACCCCTac             >  2:617124/1‑140 (MQ=255)
                                                                                                                             tAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTac             <  1:176143/140‑1 (MQ=255)
                                                                                                                             tAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTCCCGCCCTTAATCGACCTag             >  2:153420/1‑139 (MQ=255)
                                                                                                                              aGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGCTTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCTTCAATGGTTGCAGCCATTACTACACCTccc            >  1:146408/1‑140 (MQ=255)
                                                                                                                              aGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACCCCTacc            >  1:392369/1‑140 (MQ=255)
                                                                                                                                aaaaTAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCGGCCATTAATACACCTaccac          >  1:42155/1‑140 (MQ=255)
                                                                                                                                        tCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAgg  <  1:124169/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
CAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGG  >  USA300TCH1516_ALE/2372419‑2372694

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: