Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,963 T→C 6.3% I43V (ATT→GTT)  yfbT ← sugar phosphatas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,411,9630TC 6.3% 213.0 / 3.9 63I43V (ATT→GTT) yfbTsugar phosphatas
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (2/2);  total (32/31)
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.71e-01

CGGAAGCGCGGTAATACCTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCTTTGCACCGCACGCGGGCCTCCTCAGGCATTG  >  NC_000913/2411832‑2412109
                                                                                                                                   |                                                                                                                                                  
cGGAAGCGCGGTAATACCTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCa                                                                                                                                         <  4:263322/143‑1 (MQ=255)
cGGAAGCGCGGTAATACCTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCa                                                                                                                                         >  3:263322/1‑143 (MQ=255)
                ccTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAg                                                                                                                                               <  4:614991/121‑1 (MQ=255)
                ccTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAg                                                                                                                                               >  3:614991/1‑121 (MQ=255)
                 cTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAACGAAAGCCAGCACCTCTTCCGGCGCTa                                                                                                                          <  1:171793/141‑1 (MQ=255)
                 cTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAACGAAAGCCAGCACCTCTTCCGGCGCTa                                                                                                                          >  2:171793/1‑141 (MQ=255)
                  ttCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGtgaatg                                                                                                                                                   <  3:1074042/115‑1 (MQ=255)
                  ttCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGtgaatg                                                                                                                                                   >  4:1074042/1‑115 (MQ=255)
                   tCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGc                                                                                                                                              <  2:798179/119‑1 (MQ=255)
                   tCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGc                                                                                                                                              >  1:798179/1‑119 (MQ=255)
                                  tCGTTGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGc                                                                                                                                              <  2:287983/104‑1 (MQ=255)
                                  tCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGc                                                                                                                                              >  1:287983/1‑104 (MQ=255)
                                  tCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTc                                                                                                              >  2:122825/1‑136 (MQ=255)
                                  tCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTc                                                                                                              <  1:122825/136‑1 (MQ=255)
                                  tCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGc                                                                                                  >  3:715318/1‑148 (MQ=255)
                                    gATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTc                                                                                                <  4:715318/148‑1 (MQ=255)
                                            ccAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACctc                                                                                                                                     <  1:557117/103‑1 (MQ=255)
                                            ccAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACctc                                                                                                                                     >  2:557117/1‑103 (MQ=255)
                                                   cgTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCt                                                                                                             >  1:493812/1‑120 (MQ=255)
                                                   cgTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCt                                                                                                             <  2:493812/120‑1 (MQ=255)
                                                     tAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAAccc                                                                                                                      >  2:57735/1‑109 (MQ=255)
                                                     tAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAAccc                                                                                                                      <  1:57735/109‑1 (MQ=255)
                                                     tAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGtt                                                                                      >  1:722033/1‑141 (MQ=255)
                                                     tAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGtt                                                                                      <  2:722033/141‑1 (MQ=255)
                                                      aaaCTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTc                                                                                                                                  <  3:775709/96‑1 (MQ=255)
                                                      aaaCTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTc                                                                                                                                  >  4:775709/1‑96 (MQ=255)
                                                           cggcggCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTg                                                                                                            <  3:4546/113‑1 (MQ=255)
                                                           cggcggCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTg                                                                                                            >  4:4546/1‑113 (MQ=255)
                                                                gcAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAAt                                                                    >  1:1022724/1‑148 (MQ=255)
                                                                        aGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGc                                                            >  2:1064316/1‑148 (MQ=255)
                                                                             cGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTaa                                                                                                                         <  2:22731/82‑1 (MQ=255)
                                                                             cGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTaa                                                                                                                         >  1:22731/1‑82 (MQ=255)
                                                                             cGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACg                                                                                                                >  4:356997/1‑91 (MQ=255)
                                                                             cGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACg                                                                                                                <  3:356997/91‑1 (MQ=255)
                                                                                 tttGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACctct                                                                                                                                    <  2:146139/67‑1 (MQ=255)
                                                                                 tttGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACctct                                                                                                                                    >  1:146139/1‑67 (MQ=255)
                                                                                    gCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGa                                                <  2:1022724/148‑1 (MQ=255)
                                                                                      ccGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATc                                              >  2:613008/1‑148 (MQ=255)
                                                                                        gCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCaa                                            <  1:613008/148‑1 (MQ=255)
                                                                                        gCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCaa                                            <  1:1064316/148‑1 (MQ=255)
                                                                                        gCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAACGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGc                                                                                                  <  4:138656/94‑1 (MQ=255)
                                                                                        gCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAACGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGc                                                                                                  >  3:138656/1‑94 (MQ=255)
                                                                                         ccATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTa                                                                                   >  4:259476/1‑108 (MQ=255)
                                                                                         ccATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTa                                                                                   <  3:259476/108‑1 (MQ=255)
                                                                                          cattttaTGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGtt                                                         <  4:816023/127‑1 (MQ=255)
                                                                                          cattttaTGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGtt                                                         >  3:816023/7‑133 (MQ=255)
                                                                                            tAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTa                                                                                   <  3:392352/105‑1 (MQ=255)
                                                                                            tAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTa                                                                                   >  4:392352/1‑105 (MQ=255)
                                                                                                         agagGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCgcaggca                                                                          >  3:954575/1‑101 (MQ=255)
                                                                                                         agagGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCgcaggca                                                                          <  4:954575/101‑1 (MQ=255)
                                                                                                                aTCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCa                                                   >  1:413746/1‑117 (MQ=255)
                                                                                                                aTCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCa                                                   <  2:413746/117‑1 (MQ=255)
                                                                                                                aTCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCTTTGCACCGCAcgc                    >  3:125715/1‑148 (MQ=255)
                                                                                                                  cGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGa                                                                      >  3:1062422/1‑96 (MQ=255)
                                                                                                                  cGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGa                                                                      <  4:1062422/96‑1 (MQ=255)
                                                                                                                  cGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCaa                                            >  2:218008/1‑122 (MQ=255)
                                                                                                                  cGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCaa                                            <  1:218008/122‑1 (MQ=255)
                                                                                                                    ccTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCtt                                >  2:297061/1‑132 (MQ=255)
                                                                                                                    ccTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCtt                                <  1:297061/132‑1 (MQ=255)
                                                                                                                        tttACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAg                                                             <  4:1141966/99‑1 (MQ=255)
                                                                                                                        tttACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAg                                                             >  3:1141966/1‑99 (MQ=255)
                                                                                                                             cGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGaaa                                     <  3:976695/118‑1 (MQ=255)
                                                                                                                             cGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGaaa                                     >  4:976695/1‑118 (MQ=255)
                                                                                                                              gTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATcca                                                                 <  2:1016078/89‑1 (MQ=255)
                                                                                                                              gTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATcca                                                                 >  1:1016078/1‑89 (MQ=255)
                                                                                                                                  atgaaAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCTTTGCACCGCACGCGGGCCTCCTCAGGCATTg  >  2:803565/1‑148 (MQ=255)
                                                                                                                                   tgaaAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCTTTGCAc                           <  1:214592/122‑1 (MQ=255)
                                                                                                                                   tgaaAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCTTTGCAc                           >  2:214592/1‑122 (MQ=255)
                                                                                                                                   |                                                                                                                                                  
CGGAAGCGCGGTAATACCTTCGGTTTCCGTGGCCTCGATGTGCTCCAGACGCGTAAACTCGGCGGCAATATCAGCCTCGGATTTGCCCGCCATAAAATGGCGCAGAGAGGTGATCGCCTGTTTACCGTGAATGAAAGCCAGCACCTCTTCCGGCGCTAACCCATGACGTCTGGCCCAGTTGCTCCACGCCCGTTCTACCGCAGGCAGGGAATCCACCAGCGTTCCATCAAGATCAAACAGAAAACCTTTGCACCGCACGCGGGCCTCCTCAGGCATTG  >  NC_000913/2411832‑2412109

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