Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,416,597 A→T 8.5% intergenic (+87/+28) pdxB → / ← PP_2118 erythronate‑4‑phosphate dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,416,5970AT 8.5% 175.3 / 7.9 59intergenic (+87/+28)pdxB/PP_2118erythronate‑4‑phosphate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base T (2/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

ACGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATA  >  NC_002947/2416452‑2416740
                                                                                                                                                 |                                                                                                                                               
aCGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAgcg                                                                                                                                             >  1:1213038/1‑150 (MQ=255)
aCGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGACAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAgcg                                                                                                                                             >  2:720779/1‑150 (MQ=255)
       ggCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGtt                                                                                                                                      <  2:792080/150‑1 (MQ=255)
           ccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCag                                                                                                                                  >  2:758316/1‑150 (MQ=255)
              cagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAgaga                                                                                                                               >  1:1246229/1‑150 (MQ=255)
                 ttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGACCTCTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAAtt                                                                                                                            >  1:601302/1‑150 (MQ=255)
                     cgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGcc                                                                                                                        <  2:778849/150‑1 (MQ=255)
                                  tGCCCTGGGTTGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGCTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCCCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGc                                                                                                           <  1:1227213/150‑1 (MQ=255)
                                  tGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTCGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCa                                                                                                                                   >  1:322742/1‑126 (MQ=255)
                                  tGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTCGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCa                                                                                                                                   <  2:322742/126‑1 (MQ=255)
                                      cTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTAc                                                                                                       >  2:792194/1‑150 (MQ=255)
                                        gggtggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGTTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTAcgc                                                                                                     <  1:410755/150‑1 (MQ=255)
                                                 gtTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGacac                                                                                            <  1:758316/150‑1 (MQ=255)
                                                    ggTCTAAGCACGCCTTGTCTTGACCTTGCAAGTAAATGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGCCACGGc                                                                                         >  2:157052/1‑150 (MQ=255)
                                                     gTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGaa                                                                                                                              >  2:320838/1‑112 (MQ=255)
                                                     gTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGaa                                                                                                                              <  1:320838/112‑1 (MQ=255)
                                                     gTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCt                                                                                        >  1:1024465/1‑150 (MQ=255)
                                                        tAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttc                                                                                     <  1:423361/150‑1 (MQ=255)
                                                        tAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttc                                                                                     <  1:423368/150‑1 (MQ=255)
                                                             aCGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCa                                                                                <  2:620644/150‑1 (MQ=255)
                                                              cGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGACAAGCCCACTGGCTGGGCGGGTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAg                                                                               <  1:720779/150‑1 (MQ=255)
                                                               gCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTCGCGTACGCGCTCGacac                                                                                            <  1:1184723/136‑1 (MQ=255)
                                                               gCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTCGCGTACGCGCTCGacac                                                                                            >  2:1184723/1‑136 (MQ=255)
                                                                    gTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCCCTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCa                                                                                                                                   <  2:1277688/92‑1 (MQ=255)
                                                                    gTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCa                                                                                                                                   >  1:1277688/1‑92 (MQ=255)
                                                                     tCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTg                                                                        >  1:99294/1‑150 (MQ=255)
                                                                      cTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGGCGGGTTGCAGACAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGc                                                                       >  1:664483/1‑150 (MQ=255)
                                                                       ttGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGcc                                                                                                                        <  2:755460/100‑1 (MQ=255)
                                                                       ttGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGcc                                                                                                                        >  1:755460/1‑100 (MQ=255)
                                                                          aCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTg                                                                        <  1:773814/145‑1 (MQ=255)
                                                                          aCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTg                                                                        >  2:773814/1‑145 (MQ=255)
                                                                            cTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCct                                                                 >  2:153550/1‑150 (MQ=255)
                                                                                    taggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTTGatc                                                         <  1:114074/150‑1 (MQ=255)
                                                                                     aggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatc                                                         <  2:890827/149‑1 (MQ=255)
                                                                                     aggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatc                                                         >  1:890827/1‑149 (MQ=255)
                                                                                        tgtgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAggg                                                                             >  1:318001/1‑126 (MQ=255)
                                                                                        tgtgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAggg                                                                             <  2:318001/126‑1 (MQ=255)
                                                                                           ggTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCt                                                  >  2:820499/1‑150 (MQ=255)
                                                                                                  tGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGAt                                           >  2:1182663/1‑150 (MQ=255)
                                                                                                    ttGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCg                                         >  1:1169585/1‑150 (MQ=255)
                                                                                                     tGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCgg                                        <  2:1169585/150‑1 (MQ=255)
                                                                                                         gcggcCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGa                                                           >  2:1282869/1‑127 (MQ=255)
                                                                                                         gcggcCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGa                                                           <  1:1282869/127‑1 (MQ=255)
                                                                                                             ccAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGa                                                           >  2:15872/1‑123 (MQ=255)
                                                                                                             ccAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGa                                                           <  1:15872/123‑1 (MQ=255)
                                                                                                              cAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCCGCTGGATCATCTGCTCGGTGATCGGCACTTcgcg                               >  1:664412/1‑150 (MQ=255)
                                                                                                                  cccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcccc                           >  1:188475/1‑150 (MQ=255)
                                                                                                                   ccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCt                          >  2:1229210/1‑150 (MQ=255)
                                                                                                                         gctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGt                    <  1:820499/150‑1 (MQ=255)
                                                                                                                               ccAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTc                   >  2:49548/1‑145 (MQ=255)
                                                                                                                               ccAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTc                   <  1:49548/145‑1 (MQ=255)
                                                                                                                                  tctgcaaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCt                                                     >  1:986181/6‑108 (MQ=39)
                                                                                                                                  tctgcaaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCt                                                     <  2:986181/103‑1 (MQ=39)
                                                                                                                                   tAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCg                                         <  2:709068/119‑1 (MQ=255)
                                                                                                                                   tAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCg                                         >  1:709068/1‑119 (MQ=255)
                                                                                                                                    aaaaaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGaa         <  2:99294/150‑1 (MQ=255)
                                                                                                                                         cccGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttcttc                                                                                  >  2:61348/1‑72 (MQ=39)
                                                                                                                                         cccGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttcttc                                                                                  <  1:61348/72‑1 (MQ=39)
                                                                                                                                           cGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATa  <  2:1246229/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATA  >  NC_002947/2416452‑2416740

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