Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,783,086 G→A G192G (GGC→GGT gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,783,0860GA91.4% 186.2 / 8.8 70G192G (GGC→GGTgpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (35/29);  total (38/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 = 1.42e-01

CAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGT  >  NC_000913/3782950‑3783212
                                                                                                                                        |                                                                                                                              
cAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGcacc                                                                                                                              >  2:536954/1‑139 (MQ=255)
           cccAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCa                                                                                                                   <  2:251961/139‑1 (MQ=255)
                cgccgcACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCa                                                                                                              <  1:406329/139‑1 (MQ=255)
                   cgcACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATg                                                                                                           <  1:536954/139‑1 (MQ=255)
                          aGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATGCCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCg                                                                                                    >  2:1183117/1‑139 (MQ=255)
                            aCGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGa                                                                                                  <  1:1629248/139‑1 (MQ=255)
                                    tttCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTa                                                                                          <  2:805520/139‑1 (MQ=255)
                                      tCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTaaa                                                                                        >  1:1051399/1‑139 (MQ=255)
                                              cccACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGaaa                                                                                >  1:471945/1‑139 (MQ=255)
                                              cccACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGaaa                                                                                >  1:317440/1‑139 (MQ=255)
                                                        aTCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTg                                                                                                                     >  1:1440006/1‑92 (MQ=255)
                                                        aTCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTg                                                                                                                     <  2:1440006/92‑1 (MQ=255)
                                                        aTCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACg                                                                                                                 >  1:118167/1‑96 (MQ=255)
                                                        aTCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACg                                                                                                                 <  2:118167/96‑1 (MQ=255)
                                                             cgcCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgc                                                                 >  2:832455/1‑139 (MQ=255)
                                                               ccGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcag                                                               >  2:597178/1‑139 (MQ=255)
                                                                  tACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg                                                                       <  2:389592/128‑1 (MQ=255)
                                                                  tACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg                                                                       >  1:389592/1‑128 (MQ=255)
                                                                  tACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGcagcag                                                            >  1:510655/1‑139 (MQ=255)
                                                                        cATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAgctgc                                                        >  2:498535/1‑137 (MQ=255)
                                                                        cATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAgctgc                                                        <  1:498535/137‑1 (MQ=255)
                                                                           tCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGgag                                                   >  2:114222/1‑139 (MQ=255)
                                                                           tCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGgag                                                   <  1:114222/139‑1 (MQ=255)
                                                                               aCCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGatca                                               >  2:1492495/1‑139 (MQ=255)
                                                                                     aCCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCa                                         <  2:1051399/139‑1 (MQ=255)
                                                                                               tCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGAt                               <  2:317440/139‑1 (MQ=255)
                                                                                                 ggACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTg                                                                                                                     <  1:1117531/51‑1 (MQ=255)
                                                                                                 ggAAATCCCCGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTg                                                                                                                     >  2:1117531/1‑51 (MQ=38)
                                                                                                       ccccGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAg                       >  1:435356/1‑139 (MQ=255)
                                                                                                       ccccGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAg                       >  2:959036/1‑139 (MQ=255)
                                                                                                        cccGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAgg                      >  2:335168/1‑139 (MQ=255)
                                                                                                         ccGCACCAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGc                     <  1:832455/139‑1 (MQ=255)
                                                                                                            cACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAg                  <  1:871901/139‑1 (MQ=255)
                                                                                                               cAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg                                                                       <  2:968273/83‑1 (MQ=255)
                                                                                                               cAATGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg                                                                       >  1:968273/1‑83 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCgg                                                                                                   <  1:183039/53‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCgg                                                                                                   >  2:183039/1‑53 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGatca                                               <  1:1486677/105‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGatca                                               >  2:1486677/1‑105 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGc                                                                         >  2:598883/1‑79 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGc                                                                         <  1:598883/79‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg                                                                       >  2:1029130/1‑81 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg                                                                       <  1:1029130/81‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcagc                                                              >  1:341352/1‑90 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcagc                                                              <  2:341352/90‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGatcat                                              >  2:1586732/1‑106 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGatcat                                              <  1:1586732/106‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCa                                         <  1:1522850/111‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCa                                         >  2:1522850/1‑111 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGAt                               >  1:587662/1‑121 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGAt                               <  2:587662/121‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGAt                               >  2:81138/1‑121 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGAt                               <  1:81138/121‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCgg                            <  1:256343/124‑1 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCgg                            >  2:256343/1‑124 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGaaa             >  1:1215298/1‑139 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGaaa             >  2:1622072/1‑139 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGaaa             >  1:1267623/1‑139 (MQ=255)
                                                                                                                 aTGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGaaa             >  1:828346/1‑139 (MQ=255)
                                                                                                                  tGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGa                                                                                                  <  1:1082324/53‑1 (MQ=255)
                                                                                                                  tGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGa                                                                                                  >  2:1082324/1‑53 (MQ=255)
                                                                                                                  tGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAAt            >  2:279105/1‑139 (MQ=255)
                                                                                                                  tGGCAATAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAAt            >  1:1272553/1‑139 (MQ=255)
                                                                                                                      aaTAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGc                                                                      >  2:884642/1‑77 (MQ=255)
                                                                                                                      aaTAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGc                                                                      <  1:884642/77‑1 (MQ=255)
                                                                                                                      aaTAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCt        <  2:510655/139‑1 (MQ=255)
                                                                                                                      aaTAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCt        <  1:597178/139‑1 (MQ=255)
                                                                                                                        tAACGTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGt      >  1:461068/1‑139 (MQ=255)
                                                                                                                            gTTTTTCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGt  <  2:435356/139‑1 (MQ=255)
                                                                                                                               tttCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCa                                         <  2:494299/97‑1 (MQ=255)
                                                                                                                               tttCACCGCACCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCa                                         >  1:494299/1‑97 (MQ=255)
                                                                                                                                        |                                                                                                                              
CAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGT  >  NC_000913/3782950‑3783212

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