Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 237,540 C→A 7.7% A69E (GCG→GAG) ‡ yafT → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913237,5400CA 7.7% 214.5 / 3.8 63A69E (GCG→GAG) ‡yafTlipoprotein
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base A (3/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGTATCCCG  >  NC_000913/237399‑237676
                                                                                                                                             |                                                                                                                                        
aaTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGAt                                                                                                                                   <  4:201353/149‑1 (MQ=255)
 aTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATg                                                                                                                                  <  2:441653/149‑1 (MQ=255)
            atgaAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTc                                                                                                                       <  1:223460/149‑1 (MQ=255)
               aaaaGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    >  2:402561/1‑149 (MQ=255)
               aaaaGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    >  2:321220/1‑149 (MQ=255)
                   gaagaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTc                                                                                                                >  3:216908/1‑149 (MQ=255)
                     agaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                                                                                                              <  3:73657/149‑1 (MQ=255)
                      gaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTcccc                                                                                                             >  3:80103/1‑149 (MQ=255)
                        aaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                                                                                                           >  3:55628/1‑149 (MQ=255)
                         aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                                                                                                              >  1:271954/1‑145 (MQ=255)
                         aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                                                                                                              <  2:271954/145‑1 (MQ=255)
                         aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaata                                                                                                                                 >  4:341849/1‑123 (MQ=255)
                         aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaata                                                                                                                                 <  3:341849/126‑4 (MQ=255)
                                 aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt                                                                                                  <  4:207441/149‑1 (MQ=255)
                                 aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaa                                                                                                                                   >  2:263536/1‑115 (MQ=255)
                                 aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaa                                                                                                                                   <  1:263536/116‑2 (MQ=255)
                                       cacaAAGTAATATTAATAAAAACACTCCCAAGCAACTGACAGACAAAGATTTACTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATa                                                                                            <  2:269787/149‑1 (MQ=255)
                                        acaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATAc                                                                                           <  1:321220/149‑1 (MQ=255)
                                        acaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGAGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATAc                                                                                           >  3:299442/1‑149 (MQ=255)
                                          aaaGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACtt                                                                                         >  3:63915/1‑149 (MQ=255)
                                                atatTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   >  3:456264/1‑133 (MQ=255)
                                                atatTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   <  4:456264/133‑1 (MQ=255)
                                                   ttAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    <  4:200119/113‑1 (MQ=255)
                                                   ttAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    >  3:200119/1‑113 (MQ=255)
                                                       taaAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCa                                                                            >  1:410841/1‑149 (MQ=255)
                                                         aaaaCACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                                                                                                      >  1:319102/1‑105 (MQ=255)
                                                         aaaaCACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                                                                                                      <  2:319102/105‑1 (MQ=255)
                                                          aaaCACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAAcc                                                                         >  3:11023/1‑149 (MQ=255)
                                                             cacTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACcgcg                                                                      >  3:162346/1‑149 (MQ=255)
                                                                  cccAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTaa                                                                                              >  3:219472/1‑120 (MQ=255)
                                                                  cccAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTaa                                                                                              <  4:219472/120‑1 (MQ=255)
                                                                    cAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGa                                                               >  2:209676/1‑149 (MQ=255)
                                                                         aCTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCAc                                                                   <  2:335230/140‑1 (MQ=255)
                                                                         aCTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCAc                                                                   >  1:335230/1‑140 (MQ=255)
                                                                            gacagacaAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTa                                                       <  4:162346/149‑1 (MQ=255)
                                                                             acagacaAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTAt                                                      <  4:55628/149‑1 (MQ=255)
                                                                                gacaAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCa                                                   >  4:323201/1‑149 (MQ=255)
                                                                                  caAAGATTTATTCAGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCagg                                                 <  4:200853/149‑1 (MQ=255)
                                                                                   aaaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCg                                                                               <  2:423149/118‑1 (MQ=255)
                                                                                   aaaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCg                                                                               >  1:423149/1‑118 (MQ=255)
                                                                                   aaaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCagga                                                >  4:282413/1‑149 (MQ=255)
                                                                                    aaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   >  3:25684/1‑97 (MQ=255)
                                                                                    aaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   <  4:25684/97‑1 (MQ=255)
                                                                                      gATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAggagga                                             <  4:80103/149‑1 (MQ=255)
                                                                                        tttATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGgaga                                           >  2:375225/1‑149 (MQ=255)
                                                                                         ttATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGAt                                          <  1:209676/149‑1 (MQ=255)
                                                                                         ttATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATAAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTAAGCAGGAGGAGAt                                          <  4:299442/149‑1 (MQ=255)
                                                                                            ttCGGTAATGAAACCACTCTGGCTGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGTAAAc                                                            >  1:449575/1‑128 (MQ=255)
                                                                                            ttCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCg                                       >  1:147833/1‑149 (MQ=255)
                                                                                            ttCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCCCACCGGAAAc                                                            <  2:449575/128‑1 (MQ=255)
                                                                                                    tGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTGGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAAtatt                               >  1:416439/1‑149 (MQ=255)
                                                                                                     gAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAAtatt                               >  4:22649/1‑148 (MQ=255)
                                                                                                     gAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAAtatt                               <  3:22649/148‑1 (MQ=255)
                                                                                                                 gCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCcac                  >  3:455174/1‑149 (MQ=255)
                                                                                                                     tATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATtc              >  2:94951/1‑149 (MQ=255)
                                                                                                                      aTCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATtct             >  1:214382/1‑149 (MQ=255)
                                                                                                                       tCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATAGTGTTTCCACATtctc            <  4:455174/149‑1 (MQ=255)
                                                                                                                         cGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTg          <  1:402561/149‑1 (MQ=255)
                                                                                                                          gAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTgg         >  2:310364/1‑149 (MQ=255)
                                                                                                                           aaGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGt        <  4:216908/149‑1 (MQ=255)
                                                                                                                                gATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGt                                                                                       <  4:211097/65‑1 (MQ=255)
                                                                                                                                gATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGt                                                                                       >  3:211097/1‑65 (MQ=255)
                                                                                                                                 atatTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGTATCCCg  <  1:375225/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                        
AATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGTATCCCG  >  NC_000913/237399‑237676

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