Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,513,698 C→T 7.3% N294N (AAC→AAT ydcT → putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,513,6980CT 7.3% 160.6 / 2.9 55N294N (AAC→AATydcTputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (20/30);  new base T (2/2);  total (22/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 = 3.06e-01

TGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGA  >  NC_000913/1513550‑1513840
                                                                                                                                                    |                                                                                                                                              
tGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAAc                                                                                                                                                <  5:77773/149‑1 (MQ=255)
 ggCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTTAAGTGCAGGCCACTGTCACGATCCACGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGTAc                                                                                                                                                <  8:384960/148‑1 (MQ=255)
  gCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAAcg                                                                                                                                               <  1:298650/148‑1 (MQ=255)
              tgtgGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAAc                                                                                                                                   <  7:101711/148‑1 (MQ=255)
                    aTGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGGAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTtgt                                                                                                                            <  5:88762/149‑1 (MQ=255)
                    aTGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGATATTGAACGGCGGTGAAAAACTGCTtgt                                                                                                                            <  5:344367/149‑1 (MQ=255)
                      gACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTtgt                                                                                                                            >  8:291659/1‑147 (MQ=255)
                      gACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTtgt                                                                                                                            >  3:372746/1‑147 (MQ=255)
                       aCGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGa                                                                                                                          <  5:166356/148‑1 (MQ=255)
                       aCGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGa                                                                                                                          <  3:87722/148‑1 (MQ=255)
                       aCGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGa                                                                                                                          <  3:62151/148‑1 (MQ=255)
                          ggAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCa                                                                                                                      <  7:291659/149‑1 (MQ=255)
                           gAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGCACGGCGGTGAAAAACTGCTGGTGAGTCAg                                                                                                                     >  3:218756/1‑149 (MQ=255)
                           gAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAg                                                                                                                     >  7:247063/1‑149 (MQ=255)
                             aGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGACCTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGc                                                                                                                   <  3:93764/149‑1 (MQ=255)
                             aGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGc                                                                                                                   <  8:377365/149‑1 (MQ=255)
                               cTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCa                                                                                                                 <  8:223236/149‑1 (MQ=255)
                               cTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCATGGTGAACTGCAGCCAAATGGCCCGATACAGGAGGTGCAATATCAGCGCGCGGCAACTCGTTTTGAACTGAAAAAGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGcc                                                                                                                  >  1:170420/1‑148 (MQ=255)
                                  cGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGGGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAat                                                                                                               <  6:347432/148‑1 (MQ=255)
                                  cGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCAAGGCGGTGCAATATCAGGGCGCGGCCTCTCGTTTTGAAATGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAata                                                                                                              >  2:191717/1‑149 (MQ=255)
                                    cccTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAatat                                                                                                             >  3:119132/1‑148 (MQ=255)
                                        gCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCtt                                                                                                                              <  7:406446/127‑1 (MQ=255)
                                         cGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAg                                                                                                        <  1:139765/148‑1 (MQ=255)
                                           aCCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAgg                                                                                                                    >  8:76393/1‑134 (MQ=255)
                                           aCCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAgg                                                                                                                    <  7:76393/134‑1 (MQ=255)
                                            ccGGAACATATCCGCCTCAACACCCCTGGTGCACTGCAGGCCAATGGCACGATCCAGGCGGAGCAATATCAGGGCGCGGCAACTCGTTGTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCg                                                                                                     >  1:306205/1‑148 (MQ=255)
                                             cGGAACATATCCGCCTCAACACCCCTGGTTAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGTCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCgaa                                                                                                   <  7:243007/149‑1 (MQ=255)
                                                  cATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATTTCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAAct                                                                                              <  5:336884/149‑1 (MQ=255)
                                                          ccTCAACTCCCCGGGCGACCTGCAGGCCAAGCGCACGAGCCAGGCGGAGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGCTGCAAAACTGCTTGTGAGTCAGGCGAATATGACAGGCGAAGCACTGCCTGCCa                                                                                      <  2:72742/149‑1 (MQ=255)
                                                                   cccTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAAATGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGc                                                                             >  4:140968/1‑149 (MQ=255)
                                                                    ccTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGcc                                                                            >  7:68000/1‑149 (MQ=255)
                                                                     cTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGccc                                                                           >  4:270114/1‑149 (MQ=255)
                                                                        gTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCg                                                                          <  6:128319/147‑1 (MQ=255)
                                                                        gTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGa                                                                        <  8:264096/149‑1 (MQ=255)
                                                                           aaCTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGacaa                                                                     >  3:84472/1‑149 (MQ=255)
                                                                            aCTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGacaac                                                                    <  4:84472/149‑1 (MQ=255)
                                                                                cAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGTCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAgg                                                                 <  5:151545/148‑1 (MQ=255)
                                                                                cAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGt                                                                >  6:88623/1‑149 (MQ=255)
                                                                                 aGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAActgcc                                                                                           >  4:343075/1‑121 (MQ=255)
                                                                                        tGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAATGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGttt                                                        >  7:45656/1‑149 (MQ=255)
                                                                                        tGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAATGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGttt                                                        >  7:51855/1‑149 (MQ=255)
                                                                                          gCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGATATTGAACGTCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGAGAAGAAGTGCCTGCCTCGCTCACGCCCGGACAACAGGTGATGGTTTc                                                       <  2:170420/148‑1 (MQ=255)
                                                                                          gCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAATGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTc                                                       <  8:45656/148‑1 (MQ=255)
                                                                                          gCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAATGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTc                                                       <  8:51855/148‑1 (MQ=255)
                                                                                           cACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCt                                                     >  5:168302/1‑149 (MQ=255)
                                                                                               aTCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGCCAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACCGGTGATGGTTTCCTGGTc                                                 >  4:301118/1‑149 (MQ=255)
                                                                                                   aGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTg                                                    <  5:101841/142‑1 (MQ=255)
                                                                                                   aGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTg                                                    >  6:101841/1‑142 (MQ=255)
                                                                                                            aaTATCAGGGCGCGGCAACGCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGtgatgtgatg                                    <  3:37210/149‑1 (MQ=255)
                                                                                                                tCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCTGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAATAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTg                                 <  2:20120/148‑1 (MQ=255)
                                                                                                                   gggCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGc                             <  6:168302/149‑1 (MQ=255)
                                                                                                                     gcgcGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTg                           >  3:195924/1‑149 (MQ=255)
                                                                                                                      cgcgGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTgg                          <  3:143272/149‑1 (MQ=255)
                                                                                                                           cAACTCGTTTTGTACTGAAATTGAAAGGAGGTGAAAAACTGCTTTTGAGTCAGGCCAATATGACAGGCGAATAACTGCCTGCCAAGCTCACGCCCTGACAACAGGGGATGGTTTCCTGGTCGCGTGATGTGAGGGTGCCGCTGGTTgag                     <  2:306205/149‑1 (MQ=255)
                                                                                                                                  ttttGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGGTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGt              <  8:247063/149‑1 (MQ=255)
                                                                                                                                        aCTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATgg        <  8:68000/149‑1 (MQ=255)
                                                                                                                                         cTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGc       >  2:33572/1‑149 (MQ=255)
                                                                                                                                              aTTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGa  <  4:372746/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGA  >  NC_000913/1513550‑1513840

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