Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,658,893 G→A *176* (TAG→TAA gadE → gad regulon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,658,8930GA98.7% 241.7 / NA 76*176* (TAG→TAAgadEgad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/43);  minor base C (1/0);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAG  >  NC_000913/3658757‑3659021
                                                                                                                                        |                                                                                                                                
aaaatGAAATCAATATCACGTCGACGCTCAATAGTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACATGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATc                                                                                                                                <  1:1549668/139‑1 (MQ=255)
                cACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTg                                                                                                                <  2:553853/139‑1 (MQ=255)
                       aCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGa                                                                                                         <  2:1341920/139‑1 (MQ=255)
                        cGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATaa                                                                                                                       >  1:18571/1‑124 (MQ=255)
                        cGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATaa                                                                                                                       <  2:18571/124‑1 (MQ=255)
                         gCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtg                                                                                                       >  1:97768/1‑139 (MQ=255)
                          cTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                      <  2:1362090/139‑1 (MQ=255)
                          cTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                      <  2:1528599/139‑1 (MQ=255)
                               tatTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctggc                                                                                                 <  1:362165/139‑1 (MQ=255)
                                     gCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           <  2:97768/139‑1 (MQ=255)
                                     gCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           <  2:147705/139‑1 (MQ=255)
                                         caaacaTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAggg                                                                                       >  1:790364/1‑139 (MQ=255)
                                              aTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCa                                                                                  >  1:1500002/1‑139 (MQ=255)
                                                    gATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCTTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCt                                                                            <  2:856748/139‑1 (MQ=255)
                                                      ttCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACGTAAGCAACTGAAATTGAtggctggc                                                                                                 <  1:147885/116‑1 (MQ=255)
                                                      ttCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACGTAAGCAACTGAAATTGAtggctggc                                                                                                 >  2:147885/1‑116 (MQ=255)
                                                       tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                         >  1:107729/1‑139 (MQ=255)
                                                       tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                         <  2:790364/139‑1 (MQ=255)
                                                       tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                         >  2:183796/1‑139 (MQ=255)
                                                                 aaCATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAccc                                                               >  2:712665/1‑139 (MQ=255)
                                                                  aCATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                              <  1:183796/139‑1 (MQ=255)
                                                                  aCATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                              <  1:1134179/139‑1 (MQ=255)
                                                                       aTGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                         <  2:1469252/139‑1 (MQ=255)
                                                                         gTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                       <  1:392086/139‑1 (MQ=255)
                                                                         gTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                       <  1:712665/139‑1 (MQ=255)
                                                                             aaGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                         >  1:1407301/1‑133 (MQ=255)
                                                                             aaGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                         <  2:1407301/133‑1 (MQ=255)
                                                                               gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                 >  1:1029318/1‑139 (MQ=255)
                                                                                cTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           <  2:736804/96‑1 (MQ=255)
                                                                                cTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           >  1:736804/1‑96 (MQ=255)
                                                                                   aaaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                             <  2:107729/139‑1 (MQ=255)
                                                                                    aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTACTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGgag                                            >  1:931774/1‑139 (MQ=255)
                                                                                    aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATaa                                                                                                                       >  1:920769/1‑64 (MQ=255)
                                                                                    aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATaa                                                                                                                       <  2:920769/64‑1 (MQ=255)
                                                                                    aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGgag                                            >  2:1546183/1‑139 (MQ=255)
                                                                                    aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGACGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGgag                                            >  1:965946/1‑139 (MQ=255)
                                                                                         aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGa                                                                                      <  1:1554879/92‑1 (MQ=255)
                                                                                         aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGa                                                                                      >  2:1554879/1‑92 (MQ=255)
                                                                                          gACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCTGGTg                                                         >  2:1475295/1‑120 (MQ=255)
                                                                                          gACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCTGGTg                                                         <  1:1475295/120‑1 (MQ=255)
                                                                                             gTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAgg                                                                                        <  2:1081977/86‑1 (MQ=255)
                                                                                             gTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAgg                                                                                        >  1:1081977/1‑86 (MQ=255)
                                                                                              tATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAAc                                  >  1:739285/1‑139 (MQ=255)
                                                                                                     gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                           <  2:931774/139‑1 (MQ=255)
                                                                                                     gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                           <  2:739285/139‑1 (MQ=255)
                                                                                                     gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                           <  2:1500002/139‑1 (MQ=255)
                                                                                                       cATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCaa                                                                                                                   >  2:872152/1‑49 (MQ=255)
                                                                                                       cATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCaa                                                                                                                   <  1:872152/49‑1 (MQ=255)
                                                                                                       cATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTAtt                         <  2:1367719/139‑1 (MQ=255)
                                                                                                        aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCt                                                                            <  2:1085968/87‑1 (MQ=255)
                                                                                                        aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCt                                                                            >  1:1085968/1‑87 (MQ=255)
                                                                                                        aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACtt                             >  2:749348/1‑134 (MQ=255)
                                                                                                        aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACtt                             <  1:749348/134‑1 (MQ=255)
                                                                                                          cGTCACCCTTGGTATAACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTccc                      <  1:1546183/139‑1 (MQ=255)
                                                                                                            tCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                         >  1:1335625/1‑86 (MQ=255)
                                                                                                            tCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                         >  1:1411311/1‑86 (MQ=255)
                                                                                                            tCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                         <  2:1411311/86‑1 (MQ=255)
                                                                                                            tCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                         <  2:1335625/86‑1 (MQ=255)
                                                                                                                 cTGGGTATCACATCTTATTTTTAATCAGTACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAt               <  2:965946/139‑1 (MQ=255)
                                                                                                                  tGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           >  2:929826/1‑62 (MQ=255)
                                                                                                                  tGGGTATCACATCTTAGTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           <  1:929826/62‑1 (MQ=255)
                                                                                                                    ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctggc                                                                                                 <  1:1360924/54‑1 (MQ=255)
                                                                                                                    ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctggc                                                                                                 >  2:1360924/1‑54 (MQ=255)
                                                                                                                    ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCAtt                                                                   <  1:919595/84‑1 (MQ=255)
                                                                                                                    ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCAtt                                                                   >  2:919595/1‑84 (MQ=255)
                                                                                                                    ggTATCACATCTTATTTTTAACCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAtgtg            >  2:142887/1‑139 (MQ=255)
                                                                                                                        tCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATc        <  2:1029318/139‑1 (MQ=255)
                                                                                                                          acaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           >  1:1558528/1‑54 (MQ=255)
                                                                                                                          acaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                           <  2:1558528/54‑1 (MQ=255)
                                                                                                                          acaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAccc                                                               >  2:650095/1‑82 (MQ=255)
                                                                                                                          acaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAccc                                                               <  1:650095/82‑1 (MQ=255)
                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATg                                                                                    >  2:1077022/1‑57 (MQ=255)
                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAg  >  2:1219/1‑139 (MQ=255)
                                                                                                                              cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTCGCATGACGAGGGATg                                                                                    <  1:1077022/57‑1 (MQ=255)
                                                                                                                                  tttttAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTg                   <  1:1344379/118‑1 (MQ=255)
                                                                                                                                  tttttAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTg                   >  2:1344379/1‑118 (MQ=255)
                                                                                                                                        |                                                                                                                                
AAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAG  >  NC_000913/3658757‑3659021

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: