Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,117,630 C→T 8.1% M1M (ATG→ATA) † USA300TCH1516_01040 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,117,6300CT 8.1% 167.4 / 4.5 61M1M (ATG→ATA) †USA300TCH1516_01040hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base T (2/3);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

CCAAAGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTCCGCAATATAGCAA  >  USA300TCH1516_ALE/1117490‑1117770
                                                                                                                                             |                                                                                                                                            
ccAAAGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTc                                                                                                                                              >  1:168775/1‑141 (MQ=255)
  aaaGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCAt                                                                                                                                            >  1:397813/1‑141 (MQ=255)
  aaaGGTGGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAAAAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCAt                                                                                                                                            <  2:938322/141‑1 (MQ=255)
   aaGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTTGTTTCAt                                                                                                                                            >  1:506953/1‑141 (MQ=255)
          cAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATACTTTAAGATTTTTTTTGTTTCATCCTATccc                                                                                                                                    >  2:50732/1‑141 (MQ=255)
              gAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCAcc                                                                                                                                >  1:786897/1‑141 (MQ=255)
               aGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTATTGTTTCATCCTATCCCCACCt                                                                                                                               <  1:694197/141‑1 (MQ=255)
               aGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCAAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCt                                                                                                                               <  1:103816/141‑1 (MQ=255)
                 tAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCttt                                                                                                                             <  2:75254/141‑1 (MQ=255)
                  aaTATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCtttt                                                                                                                            >  2:502824/1‑141 (MQ=255)
                     ataAACGTTAAAATCATCATTAGCACACTTCATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCt                                                                                                                         <  1:611056/141‑1 (MQ=255)
                         aggTTAAAATAATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATc                                                                                                                     <  2:367725/139‑1 (MQ=255)
                         aTGTTAAAATCATCATTGGTAAATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATc                                                                                                                     <  2:343501/141‑1 (MQ=255)
                            ttAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCttt                                                                                                                  <  1:542924/141‑1 (MQ=255)
                               aaaTCATCATTTGTACATTTGATGCACCAACTGCAAAACTAAAAAATACTAAGAACACAACAGCTAATATAGAAAATACCCAAACGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTtct                                                                                                               <  2:892912/141‑1 (MQ=255)
                                   catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCCAAAGATAATTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt                                                                                                           <  1:743578/141‑1 (MQ=255)
                                   catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTTATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt                                                                                                           <  1:454704/141‑1 (MQ=255)
                                   catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTTATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt                                                                                                           <  1:387729/141‑1 (MQ=255)
                                   catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTTATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt                                                                                                           <  1:710259/141‑1 (MQ=255)
                                         tGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTTTCTTTTCTCAATTATATg                                                                                                     >  2:929182/1‑141 (MQ=255)
                                         tGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTatata                                                                                                     >  1:185358/1‑140 (MQ=255)
                                         tGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATg                                                                                                     >  1:623245/1‑141 (MQ=255)
                                                  tGATGCACCAACAGCAAAACTAAATAATAATAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATAtaataa                                                                                            <  1:791581/141‑1 (MQ=255)
                                                        accaccAGCAAAACTAAACCATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAATATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTGAATTATATGATATAATAAAAAAGt                                                                                      <  1:601390/141‑1 (MQ=255)
                                                         ccAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATTAAAAAGtt                                                                                     >  2:925831/1‑141 (MQ=255)
                                                            aCAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAATTACCAAAAAGAAATTTTTAGAATTTTTTTGTTTCTTTCTATACCCACCGTTTGCTTATCTTTACACAATAATATGATATAATACAAAAGTTgaa                                                                                  >  2:53802/1‑139 (MQ=255)
                                                             cAGGAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTaa                                                                                 <  2:376775/141‑1 (MQ=255)
                                                                   aaCTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTaaaa                                                                           <  2:929085/141‑1 (MQ=255)
                                                                   aaCTAAATAAGACTAAAAAAACAACAGCTAATATAGAAAATACCCAAAAGAAATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTaaaa                                                                           <  1:194937/141‑1 (MQ=255)
                                                                       aaaTACTACTAACAAAACAACAGCGAATATAGAAAAAACCCAAAAGATATTTTGAGATTTTTTTTGTTTCATCCTAGCCCCACCTTTTGCTTATCTTTTCTCAATTATAGGATATAATAAAAAAGTTGTAATTAAAAGTgg                                                                       <  2:530431/141‑1 (MQ=255)
                                                                           aataCTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGAt                                                                    <  2:506953/141‑1 (MQ=255)
                                                                                 aaaaaaaGAACAGCGAATATAGAAAATATCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACtt                                                             <  1:874547/138‑1 (MQ=255)
                                                                                 aaGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACtt                                                             >  1:618466/1‑141 (MQ=255)
                                                                                  aGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTa                                                            <  1:240486/141‑1 (MQ=255)
                                                                                      aacaacAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTATTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAGAGTGGGATTTTACTTaagaa                                                        >  1:533906/1‑141 (MQ=255)
                                                                                       acaacaGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTaagaaa                                                       <  1:539457/141‑1 (MQ=255)
                                                                                         aacaGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATACTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTACTTAAAAGTGGGATTTGACTTtagaaaga                                                     >  2:267138/1‑141 (MQ=255)
                                                                                           aaGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAaagaag                                                   <  1:917082/140‑1 (MQ=255)
                                                                                                aaTATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTTTTTTTTTTTATCCCCACCTTTTGCTTATCTTTTCTCAATTTTATGATATAATAAAAAAGTTGTAATTAACAGTTGGATTCTACTTAAGAAAGAAAGAAAc                                              >  2:774399/1‑141 (MQ=255)
                                                                                                  tatataAAATACCCAAAAGATATTGTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTTACTTATATGGTTTGATAAAAGAGGTGTAATTTAAAGTGGGATTTTACTTAAGGATGGAGGAAACTa                                            <  2:199125/141‑1 (MQ=255)
                                                                                                  tataGAAAATACCCAAAAGATATTTTTAGATTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTAt                                           <  1:592481/141‑1 (MQ=255)
                                                                                                    taGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTAtt                                          <  1:606180/141‑1 (MQ=255)
                                                                                                       aaaaTACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGACATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTtat                                       <  2:403067/141‑1 (MQ=255)
                                                                                                        aaaTACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTtata                                      >  1:874051/1‑141 (MQ=255)
                                                                                                           tACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGa                                   <  2:168775/141‑1 (MQ=255)
                                                                                                           aaCCCAAAAGATATTTATAGATTTTTCTTGTTTCATCCTATCCCCACCTTTTGCTTAACATATCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGa                                   <  1:925831/140‑1 (MQ=255)
                                                                                                                aaaaGATATTTTTAGATTTTTTTTTTTTCTTCTTATCCCCACATTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGGCGGATTTTACTTAAGAAAGAAGAAAACTATTTATATGACtaat                              >  2:813010/1‑141 (MQ=255)
                                                                                                                aaaaGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGaataat                              >  1:553047/1‑141 (MQ=255)
                                                                                                                aaaaGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACtaat                              <  2:191342/141‑1 (MQ=255)
                                                                                                                aaaaGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACtaat                              <  2:786897/141‑1 (MQ=255)
                                                                                                                aaaaGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTAATATGACtaat                              >  1:564748/1‑141 (MQ=255)
                                                                                                                    gATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATaaaa                          >  1:335741/1‑141 (MQ=255)
                                                                                                                       atTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGCCTAATaaaaaaa                       >  2:495514/1‑138 (MQ=255)
                                                                                                                       atTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAACTTGTAATTAAAAGTGGGATTTTACTTAACAAAGAAGGAAACTATTTATATAACTAATAAAagaa                       >  2:125890/1‑140 (MQ=255)
                                                                                                                        tttttAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga                      >  1:161931/1‑141 (MQ=255)
                                                                                                                        tttttAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga                      >  1:864091/1‑141 (MQ=255)
                                                                                                                        tttttAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga                      >  1:847188/1‑141 (MQ=255)
                                                                                                                        tttttAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAAATATTTATATGACTAGTAAAAgaga                      >  1:161795/1‑141 (MQ=255)
                                                                                                                        tttttAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGAGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga                      >  1:491048/1‑141 (MQ=255)
                                                                                                                          tttAGATTTTTTTTTTTTTTTTCTATCACCAACTTTTGCTTGTCTTTTCTCAATTATATGATATAATAAAAAAGTTGCAATTCAAAGTCGGATTTTACTTAAGAAAGATTGAAACTTTTTATATGACTAATAAAAGagaaa                    >  2:493494/1‑140 (MQ=255)
                                                                                                                                         gTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATACAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATAGCTACTAAAAAAGAAGATATCCGCacaata      >  2:122034/1‑136 (MQ=255)
                                                                                                                                          tttCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTTCGCAATATAg     <  2:638276/141‑1 (MQ=255)
                                                                                                                                             cATCCTATCCCCACCTTTTGCTTATCTTTTCTAAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTCCGCAATATAGCaa  <  2:87142/141‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CCAAAGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTCCGCAATATAGCAA  >  USA300TCH1516_ALE/1117490‑1117770

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: