Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,450,057 T→C E59E (GAA→GAG rplV ← 50S ribosomal subunit protein L22

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,450,0570TC100.0% 208.4 / NA 59E59E (GAA→GAGrplV50S ribosomal subunit protein L22
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/30);  total (29/30)

ATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGAACCTTCT  >  NC_000913/3449926‑3450190
                                                                                                                                   |                                                                                                                                     
aTGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGtt                                                                                                                               >  1:481140/1‑140 (MQ=255)
 tgtgGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTa                                                                                                                              >  1:883057/1‑140 (MQ=255)
  gtgGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAg                                                                                                                             >  1:1164999/1‑140 (MQ=255)
   tgGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGc                                                                                                                            >  2:739987/1‑140 (MQ=255)
    ggCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCa                                                                                                                           <  2:394613/140‑1 (MQ=255)
    ggCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCa                                                                                                                           <  1:902609/140‑1 (MQ=255)
        ggTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATgg                                                                                                                       >  2:959809/1‑140 (MQ=255)
        ggTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATgg                                                                                                                       >  2:387626/1‑140 (MQ=255)
          tGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCa                                                                                                                     <  2:481140/140‑1 (MQ=255)
          tGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCa                                                                                                                     <  1:589235/140‑1 (MQ=255)
               ttCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTc                                                                                                                >  1:838749/1‑140 (MQ=255)
                    gATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGaa                                                                                                           >  2:123285/1‑140 (MQ=255)
                    gATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGaa                                                                                                           <  2:794635/140‑1 (MQ=255)
                    gATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGAGCTCAGCGTTAGCAATGGCAGATTCCAGaa                                                                                                           >  2:344092/1‑140 (MQ=255)
                      tGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACt                                                                                                         >  1:451033/1‑140 (MQ=255)
                         gATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTc                                                                                                      >  1:175257/1‑140 (MQ=255)
                         gATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTc                                                                                                      <  2:524914/140‑1 (MQ=255)
                         gATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTc                                                                                                      <  2:509951/140‑1 (MQ=255)
                         gATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTc                                                                                                      <  2:883057/140‑1 (MQ=255)
                              gCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGAc                                                                                                 <  2:271074/140‑1 (MQ=255)
                               cACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGAcc                                                                                                <  1:689582/140‑1 (MQ=255)
                                      ttttGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAg                                                                                         <  2:477025/140‑1 (MQ=255)
                                              gcgGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTc                                                                                 >  1:382078/1‑140 (MQ=255)
                                              gcgGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTc                                                                                 >  1:759899/1‑140 (MQ=255)
                                              gcgGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTc                                                                                 >  2:874776/1‑140 (MQ=255)
                                              gcagCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTc                                                                                 >  1:633756/4‑140 (MQ=255)
                                                  cATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCttt                                                                                  >  1:118748/1‑135 (MQ=255)
                                                  cATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCttt                                                                                  <  2:118748/135‑1 (MQ=255)
                                                  cATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCttgt                                                                             <  2:796159/140‑1 (MQ=255)
                                                  cATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCttgt                                                                             <  1:1001355/140‑1 (MQ=255)
                                                       tGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGgtg                                                                        <  1:476931/140‑1 (MQ=255)
                                                       tGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGgtg                                                                        <  2:1164999/140‑1 (MQ=255)
                                                        gcgcTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGgtgt                                                                       >  1:712776/1‑140 (MQ=255)
                                                           cTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAgg                                                                    <  2:759899/140‑1 (MQ=255)
                                                             tCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTc                                                                  >  1:131720/1‑140 (MQ=255)
                                                                 gCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCaaaa                                                              <  2:84051/140‑1 (MQ=255)
                                                                     gggCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATc                                                          <  1:739987/140‑1 (MQ=255)
                                                                          ttCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCaaaa                                                              >  2:1164109/1‑131 (MQ=255)
                                                                          ttCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCaaaa                                                              <  1:1164109/131‑1 (MQ=255)
                                                                          ttCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCagag                                                     >  1:30865/1‑140 (MQ=255)
                                                                               cTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGc                                                >  1:184435/1‑140 (MQ=255)
                                                                                 aCGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGa                                              <  2:279470/140‑1 (MQ=255)
                                                                                 aCGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGa                                              <  2:712776/140‑1 (MQ=255)
                                                                                 aCGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGa                                              <  2:451033/140‑1 (MQ=255)
                                                                                 aCGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAAAGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGa                                              >  2:111833/1‑140 (MQ=255)
                                                                                               aaCTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACcgc                                >  1:533563/1‑140 (MQ=255)
                                                                                                   ttCAGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAAt                            <  2:175257/140‑1 (MQ=255)
                                                                                                      aGATCGTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAg                         <  2:30865/140‑1 (MQ=255)
                                                                                                           gTCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAg                    >  2:82842/1‑140 (MQ=255)
                                                                                                            tCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGc                   >  2:592874/1‑140 (MQ=255)
                                                                                                            tCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGc                   >  2:645739/1‑140 (MQ=255)
                                                                                                            tCAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGc                   >  2:659096/1‑140 (MQ=255)
                                                                                                             cAATGTCAGCGCCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGca                  <  1:310790/140‑1 (MQ=255)
                                                                                                                      cgcCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAgg            >  2:40722/1‑137 (MQ=255)
                                                                                                                      cgcCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAgg            <  1:40722/137‑1 (MQ=255)
                                                                                                                      cgcCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGa         <  1:82842/140‑1 (MQ=255)
                                                                                                                      cgcCATCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGa         <  1:959809/140‑1 (MQ=255)
                                                                                                                          aTCGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGAACCt     >  1:753058/1‑140 (MQ=255)
                                                                                                                            cGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGAACCTTc   >  2:127619/1‑140 (MQ=255)
                                                                                                                            cGTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGAACCTTc   >  1:283206/1‑140 (MQ=255)
                                                                                                                             gTTGTGCTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGAACCTTCt  <  1:645739/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
ATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTTTCTTGTTGGTGTAGGTCAAAATATCCAGAGCCTGCGACACTTTCTTACCGCGAATCAGGTCAGCAACAAGGCGAACCTTCT  >  NC_000913/3449926‑3450190

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