Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,896,445 G→A 8.7% intergenic (+78/‑130) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,896,4450GA 8.7% 158.1 / 4.5 57intergenic (+78/‑130)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (23/29);  new base A (2/3);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

ACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTT  >  NC_002947/4896297‑4896588
                                                                                                                                                    |                                                                                                                                               
acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACTTAGCGGGTAAATACCCACCGAAAGGGGCGGACCCTGCGCGGTGACCCCTAGGAGCggg                                                                                                                                                 >  7:118759/1‑149 (MQ=255)
   agaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGttt                                                                                                                                              <  5:125231/149‑1 (MQ=255)
     aagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTAc                                                                                                                                            <  6:278233/149‑1 (MQ=255)
          ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTa                                                                                                                                             <  1:183118/143‑1 (MQ=255)
          ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTa                                                                                                                                             >  2:183118/1‑143 (MQ=255)
              gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                                   <  8:18725/149‑1 (MQ=255)
                gCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                                 >  5:104342/1‑149 (MQ=255)
                         ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcag                                                                                                                        <  6:345117‑M2/149‑4 (MQ=255)
                            gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCa                                                                                                                     <  5:267700/149‑1 (MQ=255)
                                  ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGaa                                                                                                                  <  1:236149/146‑1 (MQ=255)
                                      cGCTGTGGGAATGGGTTTTTTCCCAAAGCGTGATAGGGCCTGCGCCGGCGTCTTCGCGGGTATTTCCGCTCCTACTGGCGCGGACCTTGCGCGGTCCCCCTTAGGAGCGGGTTTACCCGCTAAGAGGCCGGCACAGGc                                                                                                                      <  8:145332/138‑1 (MQ=38)
                                      cGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCCCGGTCCCCCGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGCAcaggc                                                                                                                      <  7:145332‑M2/138‑6 (MQ=255)
                                        cTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGCGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCa                                                                                                         <  1:47362/149‑1 (MQ=255)
                                           tgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaaa                                                                                                       <  2:246396/148‑1 (MQ=255)
                                                            ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTACCGGGTAAATACCCAACAAAAGGCGCGGACCGGGCGCGGTCCCCCCGAGGAGCCGGTTTACCCCCCAAGCGGCCGGCCCAGGCAGAAGCAAAAACAAACCATTTTTTGTTAtttt                                                                                      >  4:125716/1‑147 (MQ=255)
                                                                                 gcCGGCCTCTTCGCGGTAAAATCCGCTCCAACAGGCGCGGACCGTGCGCGGTCCCCCGTAGTAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta                                                                <  5:12405/149‑1 (MQ=255)
                                                                                       ctctTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGt                                                                                <  3:274953/127‑1 (MQ=255)
                                                                                       ctctTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGt                                                                                >  4:274953/1‑127 (MQ=255)
                                                                                       ctctTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCg                                                          >  6:137232/1‑149 (MQ=255)
                                                                                         ctTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                                                        >  2:245540/1‑149 (MQ=255)
                                                                                          ttCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTg                                                       >  2:331152/1‑149 (MQ=255)
                                                                                          ttCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTg                                                       >  2:220461/1‑149 (MQ=255)
                                                                                          ttCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGATATAGTCGCTg                                                       >  8:319439/1‑149 (MQ=255)
                                                                                            cgcgGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                                                        <  6:104342/146‑1 (MQ=255)
                                                                                             gcgGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGtt                                                                 <  3:267018/136‑1 (MQ=255)
                                                                                             gcgGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGtt                                                                 >  4:267018/1‑136 (MQ=255)
                                                                                               gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAAt                                                  >  7:161051/1‑149 (MQ=255)
                                                                                                 gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgg                                                     <  2:344619‑M2/144‑71 (MQ=255)
                                                                                                 gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgg                                                     >  1:344619‑M2/1‑74 (MQ=255)
                                                                                                 gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt                                                >  1:323962/1‑149 (MQ=255)
                                                                                                 gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt                                                <  8:161051/149‑1 (MQ=255)
                                                                                                   aaaTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttct                                              <  7:118287‑M2/149‑78 (MQ=255)
                                                                                                     aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttg                                            <  7:213233‑M2/149‑80 (MQ=255)
                                                                                                      tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaac                                                                            >  2:78013‑M2/1‑69 (MQ=255)
                                                                                                      tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaac                                                                            <  1:78013‑M2/116‑48 (MQ=255)
                                                                                                         gCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggc                                        >  8:118287‑M2/1‑66 (MQ=255)
                                                                                                         gCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATACACGCACCAGTTATAGTCGCTGGCAATCTGCCCCaa                                         >  4:41489/1‑148 (MQ=255)
                                                                                                         gCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtgg                                         <  8:269331‑M2/148‑83 (MQ=255)
                                                                                                         gCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttg                                                                                                        <  8:320195‑M2/85‑20 (MQ=255)
                                                                                                         gCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttg                                                                                                        >  7:320195‑M2/1‑66 (MQ=255)
                                                                                                            ccTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCATATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAAcccc                                     >  8:243621/1‑149 (MQ=255)
                                                                                                              tACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtc                                   <  3:91277‑M2/149‑89 (MQ=255)
                                                                                                                    cgcgGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGa                             <  5:345117/149‑1 (MQ=255)
                                                                                                                      cgGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGATATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc                           <  7:319439/149‑1 (MQ=255)
                                                                                                                       ggACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGaa                            <  2:323962/147‑1 (MQ=255)
                                                                                                                          ccGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAgg                       <  1:331152/149‑1 (MQ=255)
                                                                                                                           cGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcccc                                           >  4:209785/1‑128 (MQ=255)
                                                                                                                           cGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcccc                                           <  3:209785/128‑1 (MQ=255)
                                                                                                                           cGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGa                      >  8:196725/1‑149 (MQ=255)
                                                                                                                            gTGCGCTGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGAc                     <  3:125716/149‑1 (MQ=255)
                                                                                                                                gcgGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAACGCATTTTTTGTTGTTTAACCCGTATAAGCCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGACCAAGGACGCac                 >  4:95297/1‑148 (MQ=255)
                                                                                                                                      ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGa            >  1:137206/1‑148 (MQ=255)
                                                                                                                                      ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTAATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGa            >  8:174664/1‑148 (MQ=255)
                                                                                                                                       cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaa           >  7:24269/1‑148 (MQ=255)
                                                                                                                                        ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAg                        >  7:319039/1‑134 (MQ=255)
                                                                                                                                        ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAg                        <  8:319039/134‑1 (MQ=255)
                                                                                                                                          gTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaaaa       <  5:137232/149‑1 (MQ=255)
                                                                                                                                               aGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCtt  <  8:99371/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTT  >  NC_002947/4896297‑4896588

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: