Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,279,849 A→G intergenic (+21/‑131) ghxP → / → yjcE guanine/hypoxanthine permease, high affinity, guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,279,8490AG98.6% 248.2 / NA 69intergenic (+21/‑131)ghxP/yjcEguanine/hypoxanthine permease, high affinity, guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (35/33);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGG  >  NC_000913/4279714‑4279973
                                                                                                                                       |                                                                                                                            
cgcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                           <  2:769910/139‑1 (MQ=255)
cgcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATCCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                           <  2:1065942/139‑1 (MQ=255)
 gcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCggg                                                                                                                          >  1:771450/1‑139 (MQ=255)
             aTTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTg                                                                                                              >  2:28941/1‑139 (MQ=255)
               tGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGcc                                                                                                            >  1:1345202/1‑139 (MQ=255)
                    gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGttt                                                                                                       >  1:11705/1‑139 (MQ=255)
                           ttCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAtttt                                                                                                >  2:974914/1‑139 (MQ=255)
                            tCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTtc                                                                                               >  2:1197781/1‑139 (MQ=255)
                                gCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCg                                                                                           >  2:700705/1‑139 (MQ=255)
                                    aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcaca                                                                                       <  2:11705/139‑1 (MQ=255)
                                    aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcaca                                                                                       <  1:700705/139‑1 (MQ=255)
                                        cgcAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgt                                                                                   >  2:1373448/1‑139 (MQ=255)
                                                gAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCtt                                                                           <  1:1197781/139‑1 (MQ=255)
                                                          aCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                           <  1:1044036/81‑1 (MQ=255)
                                                          aCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                           >  2:1044036/1‑81 (MQ=255)
                                                           cGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTc                                                                <  2:771450/139‑1 (MQ=255)
                                                                gTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaaa                                                           >  2:1373280/1‑139 (MQ=255)
                                                                gTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaaa                                                           <  1:552896/139‑1 (MQ=255)
                                                                 ttATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAttt                                                                                                 >  2:117581/1‑100 (MQ=255)
                                                                 ttATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAttt                                                                                                 <  1:117581/100‑1 (MQ=255)
                                                                        cGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGca                                                                                         >  1:517410/1‑101 (MQ=255)
                                                                        cGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGca                                                                                         <  2:517410/101‑1 (MQ=255)
                                                                        cGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATgag                                                   >  2:274600/1‑139 (MQ=255)
                                                                        cGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATgag                                                   >  2:398814/1‑139 (MQ=255)
                                                                        cGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGACTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATgag                                                   >  2:1080653/1‑139 (MQ=255)
                                                                            gcgcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAtt                                               >  1:98063/1‑139 (MQ=255)
                                                                            gcgcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAtt                                               <  1:685706/139‑1 (MQ=255)
                                                                              gcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCa                                             <  1:1373448/139‑1 (MQ=255)
                                                                                  gTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGtt                                                                                                        <  2:390297/76‑1 (MQ=255)
                                                                                  gTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGtt                                                                                                        >  1:390297/1‑76 (MQ=255)
                                                                                  gTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCa                                         <  1:1080653/139‑1 (MQ=255)
                                                                                   tGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCAt                                        <  1:1103052/139‑1 (MQ=255)
                                                                                    gACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                       <  1:28941/139‑1 (MQ=255)
                                                                                            aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGcc                                                                                                            >  1:563620/1‑62 (MQ=255)
                                                                                            aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGcc                                                                                                            <  2:563620/62‑1 (MQ=255)
                                                                                            aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGcc                                                                                                            >  1:1118792/1‑62 (MQ=255)
                                                                                            aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGcc                                                                                                            <  2:1118792/62‑1 (MQ=255)
                                                                                                 ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                           >  1:868598/1‑42 (MQ=255)
                                                                                                 ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                           <  2:868598/42‑1 (MQ=255)
                                                                                                 ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAg                          >  1:770207/1‑139 (MQ=255)
                                                                                                  gTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                         >  1:1495100/1‑139 (MQ=255)
                                                                                                  gTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                         <  2:98063/139‑1 (MQ=255)
                                                                                                    ggCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACg                                    <  1:18364/126‑1 (MQ=255)
                                                                                                    ggCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACg                                    >  2:18364/1‑126 (MQ=255)
                                                                                                    ggCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                  >  2:1302978/1‑128 (MQ=255)
                                                                                                    ggCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                  <  1:1302978/128‑1 (MQ=255)
                                                                                                     gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcaca                                                                                       <  2:1126793/74‑1 (MQ=255)
                                                                                                     gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcaca                                                                                       >  1:1126793/1‑74 (MQ=255)
                                                                                                     gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                >  1:1045277/1‑81 (MQ=255)
                                                                                                     gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                <  2:1045277/81‑1 (MQ=255)
                                                                                                           cTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaaa                <  1:804076/139‑1 (MQ=255)
                                                                                                              tCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGc                                                                             >  2:1459377/1‑75 (MQ=255)
                                                                                                              tCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGc                                                                             <  1:1459377/75‑1 (MQ=255)
                                                                                                               cTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACAt            <  1:974914/139‑1 (MQ=255)
                                                                                                                 aaTCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCg          <  1:398814/139‑1 (MQ=255)
                                                                                                                   tCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                       >  1:1150101/1‑108 (MQ=255)
                                                                                                                   tCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                       <  2:1150101/108‑1 (MQ=255)
                                                                                                                   tCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCa        >  1:437551/1‑139 (MQ=255)
                                                                                                                        cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAtt                                               >  1:862423/1‑95 (MQ=255)
                                                                                                                        cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAtt                                               <  2:862423/95‑1 (MQ=255)
                                                                                                                        cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTGAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAg   >  2:674942/1‑139 (MQ=255)
                                                                                                                        cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                         <  2:173102/117‑1 (MQ=255)
                                                                                                                        cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                         >  1:173102/1‑117 (MQ=255)
                                                                                                                        cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAg   >  1:640623/1‑139 (MQ=255)
                                                                                                                         aTCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAgg  <  2:1495100/139‑1 (MQ=255)
                                                                                                                          tCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                       <  2:198599/101‑1 (MQ=255)
                                                                                                                          tCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                       >  1:198599/1‑101 (MQ=255)
                                                                                                                                tCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtctc                            <  1:1390040/106‑1 (MQ=255)
                                                                                                                                tCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtctc                            >  2:1390040/1‑106 (MQ=255)
                                                                                                                                       |                                                                                                                            
CGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGG  >  NC_000913/4279714‑4279973

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