Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,182,071 G→A 7.1% intergenic (‑36/‑119) PP_1935 ← / → PP_5487 Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,182,0710GA 7.1% 162.7 / 2.0 57intergenic (‑36/‑119)PP_1935/PP_5487Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base A (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CGGCAGGATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATT  >  NC_002947/2181924‑2182216
                                                                                                                                                   |                                                                                                                                                  
cggcAGGATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTAtgc                                                                                                                                                   <  4:3946/149‑1 (MQ=255)
    aGGATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGTGTCGAGGCACGAACAGGCAGTATGCTGcc                                                                                                                                               <  8:264632/149‑1 (MQ=255)
      gATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCgg                                                                                                                                             >  3:29580/1‑149 (MQ=255)
       aTCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGCATACTGCCTGt                                                                                                                                            <  1:42553/149‑1 (MQ=255)
           aaaGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACg                                                                                                                                        >  8:55528/1‑149 (MQ=255)
             aGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACTCGGCCTAAGGCATTTTTGAAATTCATCACGCGTACGCGTCCAAAAGGAAGCAGTCAACCTATTGCGTCGAGGCAGGAACAGGCAGTATGCTGCCGGTTACGTc                                                                                                                                      <  1:183470/149‑1 (MQ=255)
                    tGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCt                                                                                                                              >  2:230663/1‑149 (MQ=255)
                      cGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGACACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCt                                                                                                                              <  1:61497/147‑1 (MQ=255)
                         aaGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCAtttt                                                                                                                                 >  5:20077/1‑142 (MQ=255)
                         aaGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCAtttt                                                                                                                                 <  6:20077/142‑1 (MQ=255)
                          aGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGttt                                                                                                                        <  7:55528/149‑1 (MQ=255)
                            gCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTc                                                                                                                      >  7:197314/1‑149 (MQ=255)
                                cGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCt                                                                                                                  >  3:304167/1‑149 (MQ=255)
                                 gCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCtt                                                                                                                 >  4:212643/1‑149 (MQ=255)
                                     aaCACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTcc                                                                                                              >  6:278526/1‑148 (MQ=255)
                                     aaCACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTc                                                                                                               <  3:212643/147‑1 (MQ=255)
                                      acacGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCg                                                                                                             <  4:29580/148‑1 (MQ=255)
                                       cacGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGa                                                                                                            >  8:292840/1‑148 (MQ=255)
                                        acGCAGTTGGTGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGaa                                                                                                           >  8:7380/1‑148 (MQ=255)
                                          gCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCAAGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATc                                                                                                        >  7:300739/1‑149 (MQ=255)
                                            aGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCTGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCt                                                                                                       >  4:194938/1‑148 (MQ=255)
                                            aGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCAAGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCt                                                                                                       >  1:63536/1‑148 (MQ=255)
                                             gTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGa                                                                                                     <  8:15047/149‑1 (MQ=255)
                                                   ggAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTc                                                                                               >  8:296022/1‑149 (MQ=255)
                                                            ggCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTa                                                                                      <  1:110672/149‑1 (MQ=255)
                                                              ccAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCtt                                                                                        >  3:299119/1‑145 (MQ=255)
                                                              ccAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTaga                                                                                    <  7:7380/149‑1 (MQ=255)
                                                              ccAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTaga                                                                                    <  5:235591/149‑1 (MQ=255)
                                                                 aatgCAATTTTGTAATTCATCAACGGTAAGGGTCCAAAATGTATCATTGAACCTATTGCGTCGCGGCACGAACATGCACTATGCTTCCGGTTAAGTCATTTTCACGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagaga                                                                                  <  1:230663/146‑1 (MQ=255)
                                                                 aaGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagagag                                                                                 >  8:203144/1‑149 (MQ=255)
                                                                 aaGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTACCGCCTTTagagag                                                                                 >  5:240353/1‑149 (MQ=255)
                                                                                   atcaCGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTc                                                                <  7:296022/148‑1 (MQ=255)
                                                                                      acgcgtacgcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTg                                                             <  4:287637/148‑1 (MQ=255)
                                                                                       cgagtacgcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCACCTTTAGAGAGTCCTCGGGAACGACTTCCTGCa                                                           <  8:300739/146‑1 (MQ=255)
                                                                                              gcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagagag                                                                                 <  3:250365/120‑1 (MQ=255)
                                                                                              gcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagagag                                                                                 >  4:250365/1‑120 (MQ=255)
                                                                                               cgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCTGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCaaaaa                                                   <  3:194938/149‑1 (MQ=255)
                                                                                                 tCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACt                                                                  >  1:25740/1‑132 (MQ=255)
                                                                                                 tCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACt                                                                  <  2:25740/132‑1 (MQ=255)
                                                                                                   cAAAAGGAAGCATTCATCCTATTGCGTCGAGGCACGAACATGCAGTATTCTGCCGGTTACGTCATTTTCTTCGTTTTCGCCTTTCCGAAATCTGATTGCTCCCGCCTTTAGATAGTCCTCGGGTACTACTTCCTGCACCCAAAAATTAg                                               <  7:203144/149‑1 (MQ=255)
                                                                                                    aaaaGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAg                                               >  5:18515/1‑148 (MQ=255)
                                                                                                          aaGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagaga                                                                                  <  7:24288/107‑1 (MQ=255)
                                                                                                          aaGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagaga                                                                                  >  8:24288/1‑107 (MQ=255)
                                                                                                          aaGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCa                                        <  4:304167/149‑1 (MQ=255)
                                                                                                                  aaCCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAgtgt                                 <  2:75927/148‑1 (MQ=255)
                                                                                                                       aTTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAatat                           >  2:165720/1‑149 (MQ=255)
                                                                                                                          gCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGa                         <  1:165720/148‑1 (MQ=255)
                                                                                                                            gTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATa                       >  8:40938/1‑148 (MQ=255)
                                                                                                                                 tgacgtaaCCGGCAGCATACTGCCTGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAAcgc                  >  8:45398/7‑148 (MQ=255)
                                                                                                                                  gCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAAcgc                  >  3:75999/1‑147 (MQ=255)
                                                                                                                                      gAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCgagtga            >  7:227630/1‑149 (MQ=255)
                                                                                                                                      gAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCgagt              >  6:151621/1‑147 (MQ=255)
                                                                                                                                       aaCCGGCAGCATACTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAatat                           >  6:204665/1‑133 (MQ=255)
                                                                                                                                       aaCCGGCAGCATACTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAatat                           <  5:204665/133‑1 (MQ=255)
                                                                                                                                        aCAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAg                                               <  2:231834/112‑1 (MQ=255)
                                                                                                                                        aCAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAg                                               >  1:231834/1‑112 (MQ=255)
                                                                                                                                            gCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACtt      >  1:292317/1‑149 (MQ=255)
                                                                                                                                                tATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTttatt  <  7:166639/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGCAGGATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTCTGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATT  >  NC_002947/2181924‑2182216

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