Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,455,698 2 bp→AT coding (1216‑1217/1260 nt) lyrA → Lysostaphin resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,455,6980CA77.8% 153.9 / 47.1 63coding (1216/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (22/27);  minor base . (9/5);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
*USA300TCH1516_ALE2,455,6990AT95.2% 177.8 / ‑4.3 64H406L (CAT→CTT) lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base T (31/29);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTA  >  USA300TCH1516_ALE/2455570‑2455826
                                                                                                                                ||                                                                                                                               
ttAAACAAACATCACAAAGTGGCCAATTAACAAATGATGCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTAATGCAAAAAAAAATATCg                                                                                                                            <  2:80886/135‑1 (MQ=255)
    aCAAGCACCACAAAGTTACCATTTAACAAACGATTCAAAAGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATTCAAAAAAATTATCGaaaaa                                                                                                                       <  1:356226/135‑1 (MQ=255)
    aCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                                                                       <  1:112169/135‑1 (MQ=255)
    aCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                                                                       >  2:389099/1‑135 (MQ=255)
    aCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                                                                       >  1:455959/1‑135 (MQ=255)
    aCAAGCACCACAAAGTGACGAATTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                                                                       >  1:242420/1‑135 (MQ=255)
        gCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAATTATCGAAAAAACTga                                                                                                                  >  1:204743/1‑135 (MQ=255)
          accacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaa                                                                                                                 >  1:82224/1‑135 (MQ=255)
           tcacaATGTGATGATGAAACAAACGATTCAAGTGTAGATGAAATAGAGCAAGCATTAAAAGAACCTGCGACGTATAAAGAAGACAGACGTTCATCATTTTTAATTGATGCAAAAAAATTTTCGTAAATACTgaag                                                                                                                <  2:82227/134‑1 (MQ=255)
             acaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGACAAAACTgaagaa                                                                                                              >  2:38932/1‑135 (MQ=255)
             acaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTAACGAAAAAACTgaagaa                                                                                                              >  2:279171/1‑135 (MQ=255)
                aaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATTGAAAAAACTGAAGAACaa                                                                                                           <  1:22861/135‑1 (MQ=255)
                aaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACaa                                                                                                           >  1:672745/1‑135 (MQ=255)
                 aGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAAt                                                                                                          <  1:40463/135‑1 (MQ=255)
                 aGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAAAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAAt                                                                                                          >  1:666582/1‑135 (MQ=255)
                   tgccGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTAAAGAACAATCt                                                                                                        >  1:593897/4‑135 (MQ=255)
                   tGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                                                        >  1:497132/1‑135 (MQ=255)
                   tGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                                                        <  1:168279/135‑1 (MQ=255)
                      cGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                                                                     >  1:646169/1‑135 (MQ=255)
                        aTTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                   >  2:600153/1‑135 (MQ=255)
                                    aTTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCTACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaaaga                                                                                       <  1:610605/135‑3 (MQ=255)
                                      tCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAATACAGACGGTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                     <  2:156387/135‑1 (MQ=255)
                                         aaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAAGAATAAc                                                                                  >  2:615457/1‑135 (MQ=255)
                                                 tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCTAAAAATCTGTAGAACAATTTTCAGATAAAAACAAATAATAACAAAGCGCt                                                                          <  1:9343/135‑1 (MQ=255)
                                                            aaTCATTAAAAGAACCGGCGACTTATAAATAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                                                               <  2:480103/135‑1 (MQ=255)
                                                              tctttAAAAGAACCTGCGACTTATAAAGAAAACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                             <  2:35615/132‑1 (MQ=255)
                                                              tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGATGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTGGCTAGTACcact                                                             >  2:400576/1‑135 (MQ=255)
                                                              tCATGAAAAGAACCTGCTACTTCTAAAGAAGACAGACGTTAATCAGGTCTAATTGATGCAAAAAAATTAACGAAAAAACTGAAGAATAATCTTCAGATAAAAACAAATAGTAACAAAGCGCTTGCTAGTACctct                                                             <  2:350166/135‑1 (MQ=255)
                                                               cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTGGCTAGTACCTCtt                                                            >  2:573350/1‑135 (MQ=255)
                                                               cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAGTAACAAAGCGCTCGCTAGCACCTCtt                                                            >  2:388377/1‑135 (MQ=255)
                                                               cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCCCATGCTAGTACCTCtt                                                            >  2:141048/1‑135 (MQ=255)
                                                                  tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCATGCTAGTACCTCTTaaa                                                         >  2:265828/1‑135 (MQ=255)
                                                                    aaaGAACCTGCGACTTATAGAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAACGCGCTTGCTAGTACCTCTTaaaaa                                                       >  2:157476/1‑135 (MQ=255)
                                                                    aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAGAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                       >  2:470133/1‑135 (MQ=255)
                                                                    aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTacaaa                                                       >  2:300521/1‑135 (MQ=255)
                                                                    aaaGAACCTGCGACATATAAAGCAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAAAAAAAAAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                       <  2:123871/135‑1 (MQ=255)
                                                                     aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg                                                      >  2:389654/1‑135 (MQ=255)
                                                                      aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat                                                    <  2:204743/135‑1 (MQ=255)
                                                                       gAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat                                                    <  2:175406/135‑1 (MQ=255)
                                                                        atccTGCGACTTATCAAGAAGACAGATGTTCATCAGCTGTAACTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAATGCGCTTGCTAGTCCCTATTAAAAAgatg                                                   <  1:428839/133‑1 (MQ=255)
                                                                        aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCATATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatg                                                   <  2:544164/135‑1 (MQ=255)
                                                                        aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTCGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACACAGCGCTTGCTAGTACCTCTTAAAAAgatg                                                   <  2:372562/135‑1 (MQ=255)
                                                                        aaCCTGCGACTTAAAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACTAATCGCTTTTTAGTACCTCTTAAAAAgatg                                                   <  1:491885/135‑1 (MQ=255)
                                                                         aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatga                                                  <  1:96766/135‑1 (MQ=255)
                                                                             gctaCTTATAAAGACGACAGACGGTCATCAGTTGTAAACCCAAAAACAAAAAAATCGAAAAAACCGAAGAACAATCTTCAGATAAAAACAAATAAGAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCt                                              <  2:429063/132‑1 (MQ=255)
                                                                              cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTa                                             >  2:465952/1‑135 (MQ=255)
                                                                                 cTTATAAAGAAGAGAGACGTTCATCAGTTGTAATTTATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATTTCAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                          <  1:525442/135‑1 (MQ=255)
                                                                                  ttATAAAGAAGCCAGACGTTCATCAGTTGTAATTTATGCAAAAAAATTATCGAAAAATCTGAAGAACAATCTTCAGATAAAAACAAATAATAACTAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCaa                                         <  1:81105/135‑1 (MQ=255)
                                                                                  ttATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCaa                                         >  1:100725/1‑135 (MQ=255)
                                                                                     taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACACAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcg                                      <  2:638918/135‑1 (MQ=255)
                                                                                      aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAGATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                     <  1:342078/135‑1 (MQ=255)
                                                                                       aggaagtCAGACGTTCATCAGTTGTAATTGAGGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                    <  2:338484/133‑1 (MQ=255)
                                                                                       acgaagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                    >  2:45911/3‑135 (MQ=255)
                                                                                            gaccgacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                               >  1:188919/1‑135 (MQ=255)
                                                                                               aggcGTTCATCAGTTGGAATTGATGCAAAAAAATTATCGAAAAAACTGAATAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCACGCGCTTTTctata                            <  2:227803/132‑1 (MQ=255)
                                                                                               agacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAATAACAATCTTCAGATACATACAAATAATAACTGTGGGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCGAGCGCTTTTctata                            <  1:129565/135‑1 (MQ=255)
                                                                                                 acGTTCATCAGTTTTACTGGATGCAAAAAAAGTATCGCAAAAACTGAAGAACAAGCTTCAGATAAAAACAAATAGTAACAGAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCCCGCGCTTTTctatact                          <  2:398185/135‑1 (MQ=255)
                                                                                                 acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTCGAATGCGCCTTTctatact                          >  2:582107/1‑135 (MQ=255)
                                                                                                  cGTTCATCAGTCGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatacta                         <  2:82224/135‑1 (MQ=255)
                                                                                                    ttCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactata                       >  1:421432/1‑135 (MQ=255)
                                                                                                      catctGTTGTAATTAATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatatc                     >  1:202630/1‑134 (MQ=255)
                                                                                                      catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatata                     >  1:434669/1‑135 (MQ=255)
                                                                                                                   tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaaa        <  1:600153/135‑1 (MQ=255)
                                                                                                                          aaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTa  <  1:515777/134‑1 (MQ=255)
                                                                                                                                ||                                                                                                                               
TTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTA  >  USA300TCH1516_ALE/2455570‑2455826

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: