Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,638,365 G→A S96S (AGC→AGT bamB ← outer membrane protein assembly factor BamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,638,3650GA97.1% 108.4 / NA 35S96S (AGC→AGTbamBouter membrane protein assembly factor BamB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/19);  minor base T (0/1);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GCTTCGTTCAGCGCTTGTAACTGACCGTTACTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGT  >  NC_000913/2638123‑2638610
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
gCTTCGTGCAGCGCTTGTAACTGACCGTTACTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGa                                                                                                                                                                                                                                               <  1:157105/251‑1 (MQ=255)
    cGTTCAGCGCTTGTAACTGACCGTTACTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGAcc                                                                                                                                                                                                                                             >  1:129058/1‑249 (MQ=255)
      ttCAGCGCTTGTAACTGACCGTTACTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGACTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGAcc                                                                                                                                                                                                                                             >  1:133624/1‑247 (MQ=255)
                      gACCGTTACTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACGTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCatcatc                                                                                                                                                                                                                          <  2:409364/250‑1 (MQ=255)
                        ccGTTACTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGc                                                                                                                                                                                                                       >  1:87950/1‑251 (MQ=255)
                             aCTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTc                                                                                                                                                                                                                   <  1:50886/250‑1 (MQ=255)
                                ggTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAgc                                                                                                                                                                                                                <  1:183405/250‑1 (MQ=255)
                                     ggATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAGGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGTCAGCACAGATTTCTTTGCCATCATCCGCATTCAGCGCtttt                                                                                                                                                                                                          >  2:15381/1‑251 (MQ=255)
                                        ttAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACt                                                                                                                                                                                                       >  1:182995/1‑251 (MQ=255)
                                             cccggACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTCCCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACTGACCAGATGTCTTTGCCATCATCCGCATGCAGCGCTTTTACTAAAcc                                                                                                                                                                                                  <  2:77002/247‑1 (MQ=255)
                                               cAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAg                                                                                                                                                                                                <  1:273953/251‑1 (MQ=255)
                                                            accaccGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCAt                                                                                                                                                                                                                     >  2:404024/1‑217 (MQ=255)
                                                            accaccGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCAt                                                                                                                                                                                                                     <  1:404024/217‑1 (MQ=255)
                                                             ccaccGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGc                                                                                                                                                                                   >  1:166692/1‑250 (MQ=255)
                                                               accGGGCGCGACCGTGCTTCACCAGCAACTTTAGTTTGCAATACCACCGTACCATCGCTGGTATTCAGCGCGTCAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACCCCGCCAGTAAGTACTGCAGGCTCTTTCGAGAACCAGCCATCTTTCTCGGCCAGACTGCCAGACCCAATTTCTTTGCCCTCATCCGCATTCAGCGCTTTTCCTAAAcaacagcgatccgttgcat                                                                                                                                                                                 >  2:270118/1‑232 (MQ=255)
                                                                              gctccACCCGTGACGTTTGTTGGCTATGCCGCAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCTTGTTCGCTGCCAATGTAGACATGCCCTCCAGACACTGGCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGATCCAGCCATCTTTCACTGCCAGTCTGACAGACCATATTTCTTTGCCAGCATCCGCATTCAGCGCTTGTACTAAACCAGCGCGGTCCGCTGCATAGACAGCGTTGTCCg                                                                                                                                                                  <  2:239146/246‑1 (MQ=255)
                                                                               cTTCACCCGCGACTTTATTTTGCCATGCCACAGTACCAACGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTAGCTGCCAATGTAGACTTGCCCACCAGACACGGTCACACCGCCAGAAAATAATGCAAGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGCTTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAACCCAGCGCGTTCCGCTGCATAGCCAACGTTATCCGcc                                                                                                                                                                >  2:312023/1‑251 (MQ=255)
                                                                                 tCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGcc                                                                                                                                                                <  1:345971/249‑1 (MQ=255)
                                                                                         gACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGAt                                                                                                                                                      <  1:199159/251‑1 (MQ=255)
                                                                                                            acaGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGtt                                                                                                                                   <  1:131771/251‑1 (MQ=255)
                                                                                                                       gCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGCCAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGc                                                                                                                        >  2:384569/1‑251 (MQ=255)
                                                                                                                          ggTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTGGCCAATGCCGCTa                                                                                                                      <  1:88362/250‑1 (MQ=255)
                                                                                                                                        aaaCCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACAAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTGGCCAATGCCGCTACCAACGGAAGTGCTc                                                                                                       >  1:27940/1‑251 (MQ=255)
                                                                                                                                               cgcCTTTCCGCTGCCAATGAAGACATGCCCACCAGACACGGGCCCACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGATACCAACGGAAGTGCTCCACGCCg                                                                                                <  2:318904/251‑1 (MQ=255)
                                                                                                                                                    gttCGCTGCCAATTTAGACATGCCCACCAGACGCGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGAGTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTc                                                                                           <  2:129058/250‑1 (MQ=255)
                                                                                                                                                             cAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAAc                                                                                  <  1:48723/251‑1 (MQ=255)
                                                                                                                                                                       aTGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCa                                                                         >  2:120793/1‑250 (MQ=255)
                                                                                                                                                                                   cacGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATg                                                            >  1:281819/1‑251 (MQ=255)
                                                                                                                                                                                             ccGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCtt                                                  <  1:291412/251‑1 (MQ=255)
                                                                                                                                                                                                       gTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACAtct                                        <  1:123247/251‑1 (MQ=255)
                                                                                                                                                                                                          aTGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACAtcttct                                     <  1:81409/251‑1 (MQ=255)
                                                                                                                                                                                                             catgCTCTGTAGAGATCCAGCCATCTTTCTCGGCCAGACTGACATACCAGATTTCTTTGCCATCATCCGCATTCAGTGCTGTTACTAAACCAGCGCGGTCCGCTGCATAGACATCGTTGTCTGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAAGGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCg                                  <  2:1765/248‑1 (MQ=255)
                                                                                                                                                                                                             cAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCg                                  <  2:160946/251‑1 (MQ=255)
                                                                                                                                                                                                                   ctTTAGAGAACCAGCCATCTTTCTCGGCCAGACTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGtt                             >  1:348651/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ccAGACTGACAGAACAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGt  >  1:87171/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
GCTTCGTTCAGCGCTTGTAACTGACCGTTACTGGTGTGGATTAACACCAGACCGTCGCTGACCACCGGGCGCGAAAGTGCTTCACCCGCGACTTTAGTTTGCCATGCCACAGTACCATCGCTGGTATTCAGCGCGTAAACCTGCGCCTTTTCGCTGCCAATGTAGACATGCCCACCAGACACGGTCACACCGCCAGAAAGTAATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATAGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGT  >  NC_000913/2638123‑2638610

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