Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 14,133 C→T intergenic (+54/‑35) dnaK → / → dnaJ chaperone protein DnaK/chaperone protein DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091314,1330CT100.0% 219.5 / NA 61intergenic (+54/‑35)dnaK/dnaJchaperone protein DnaK/chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/29);  total (32/29)

CCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCC  >  NC_000913/13997‑14268
                                                                                                                                        |                                                                                                                                       
ccGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGcc                                                                                                                                      >  1:89930/1‑140 (MQ=255)
ccGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGcc                                                                                                                                      >  1:1298998/1‑140 (MQ=255)
ccGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGcc                                                                                                                                      >  1:1209384/1‑140 (MQ=255)
 cGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAACAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGccc                                                                                                                                     <  2:978444/140‑1 (MQ=255)
 cGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGccc                                                                                                                                     <  1:1211973/140‑1 (MQ=255)
      cTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCa                                                                                                                                <  2:202191/140‑1 (MQ=255)
      cTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCa                                                                                                                                <  2:89930/140‑1 (MQ=255)
       tGATGCTTCTGCAAACAACGCCAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCAt                                                                                                                               <  2:1298998/140‑1 (MQ=255)
        gATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCAtt                                                                                                                              <  2:878126/140‑1 (MQ=255)
        gATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCAtt                                                                                                                              <  2:539527/140‑1 (MQ=255)
              tCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTa                                                                                                                        >  2:252793/1‑140 (MQ=255)
                  cAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAgggg                                                                                                                    >  2:701038/1‑140 (MQ=255)
                         cgcgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTa                                                                                                             >  2:160417/1‑140 (MQ=255)
                              aaGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTaaaaaa                                                                                                        <  2:416228/140‑1 (MQ=255)
                                   gACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATgg                                                                                                   >  2:1178146/1‑140 (MQ=255)
                                   gACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATgg                                                                                                   >  1:932476/1‑140 (MQ=255)
                                        tgttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAg                                                                                              <  2:1155793/140‑1 (MQ=255)
                                          ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAg                                                                                                       >  1:187112/1‑129 (MQ=255)
                                          ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAg                                                                                                       <  2:187112/129‑1 (MQ=255)
                                          ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCa                                                                                            >  2:938972/1‑140 (MQ=255)
                                            gtCGACGCTGAATTTGAAGAAGTCAAAAACAAAAAATAATCGCCCTATAAACGGGAAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTATTTTAGGGGCATTTTAAAAAAGATGGCTAAGCAAg                                                                                          >  2:897382/1‑140 (MQ=255)
                                             tCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGa                                                                                         >  1:910758/1‑140 (MQ=255)
                                                  gCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGattatt                                                                                    >  2:1127944/1‑140 (MQ=255)
                                                    tGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTAc                                                                                  >  2:906506/1‑140 (MQ=255)
                                                     gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACg                                                                                 <  1:897382/140‑1 (MQ=255)
                                                                agTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGc                                                                      <  1:1127944/140‑1 (MQ=255)
                                                                agTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGc                                                                      >  1:65653/1‑140 (MQ=255)
                                                                agTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGc                                                                      <  1:252793/140‑1 (MQ=255)
                                                                          aaaaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa                                                            <  1:1427983/140‑1 (MQ=255)
                                                                           aaaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc                                                           >  1:782199/1‑140 (MQ=255)
                                                                                 aatCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGa                                                     <  1:1178146/140‑1 (MQ=255)
                                                                                        cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGt                                              <  1:160417/140‑1 (MQ=255)
                                                                                        cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGt                                              <  2:1342146/140‑1 (MQ=255)
                                                                                                 ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa                                     >  1:788505/1‑140 (MQ=255)
                                                                                                 ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa                                     >  1:1220372/1‑140 (MQ=255)
                                                                                                  gTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaa                                    <  1:803975/140‑1 (MQ=255)
                                                                                                  gTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaa                                    <  2:932476/140‑1 (MQ=255)
                                                                                                   tAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaaa                                   >  1:1007890/1‑140 (MQ=255)
                                                                                                     aTTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAgg                                 >  2:956637/1‑140 (MQ=255)
                                                                                                       tataCTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaa                                    <  2:756186/135‑1 (MQ=255)
                                                                                                       tataCTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaa                                    >  1:756186/1‑135 (MQ=255)
                                                                                                            tGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACaa                          >  2:453296/1‑140 (MQ=255)
                                                                                                            tGACACGGGCGAAGGGGAATTTCCTCTCTGCCCATGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACaa                          >  2:1222563/1‑140 (MQ=255)
                                                                                                              acacGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAAc                        >  2:533015/1‑140 (MQ=255)
                                                                                                                acGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGc                      >  2:1047253/1‑140 (MQ=255)
                                                                                                                   ggCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTg                   >  2:826926/1‑140 (MQ=255)
                                                                                                                    gCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                  <  1:906506/140‑1 (MQ=255)
                                                                                                                    gCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                  >  1:1161257/1‑140 (MQ=255)
                                                                                                                    gCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                  <  2:1007890/140‑1 (MQ=255)
                                                                                                                     cGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGc                      >  1:150219/1‑135 (MQ=255)
                                                                                                                     cGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGc                      <  2:150219/135‑1 (MQ=255)
                                                                                                                     cGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGc                 <  1:938972/140‑1 (MQ=255)
                                                                                                                        aGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                  <  2:445257/136‑1 (MQ=255)
                                                                                                                        aGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                  >  1:445257/1‑136 (MQ=255)
                                                                                                                         ggggAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg             <  2:1220372/140‑1 (MQ=255)
                                                                                                                          gggAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGa            >  1:42649/1‑140 (MQ=255)
                                                                                                                               tttCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATac       <  1:453296/140‑1 (MQ=255)
                                                                                                                                    tctcTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                  >  1:1223025/1‑124 (MQ=255)
                                                                                                                                    tctcTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                  <  2:1223025/124‑1 (MQ=255)
                                                                                                                                    tctcTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCAccc  >  2:222204/1‑140 (MQ=255)
                                                                                                                                    tctcTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCAccc  <  1:382540/140‑1 (MQ=255)
                                                                                                                                    tctcTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCAccc  >  2:143891/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCC  >  NC_000913/13997‑14268

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