Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,354,272 A→C 12.8% S371A (TCC→GCC) ‡ nrdA ← ribonucleoside‑diphosphate reductase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,354,2720AC12.8% 134.5 / 12.5 56S371A (TCC→GCC) ‡nrdAribonucleoside‑diphosphate reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base C (4/3);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GTGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCATGA  >  NC_002947/1354127‑1354415
                                                                                                                                                 |                                                                                                                                               
gTGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGAcccc                                                                                                                                             <  2:715403/150‑4 (MQ=255)
 tGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCt                                                                                                                                            <  1:1133125/150‑1 (MQ=255)
   cccTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGc                                                                                                                                          >  1:464276/1‑150 (MQ=255)
    ccTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCg                                                                                                                                         <  1:572612/150‑1 (MQ=255)
      ttGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCTGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGTCCAGCTGctac                                                                                                                                       <  2:1213116/150‑4 (MQ=255)
      ttGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCggg                                                                                                                                       >  1:208535/1‑150 (MQ=255)
                agCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGc                                                                                                                             >  1:221352/1‑150 (MQ=255)
                 gCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGcc                                                                                                                            <  2:475234/150‑1 (MQ=255)
                  cccAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGa                                                                                                                                                 >  2:632771/1‑128 (MQ=255)
                  cccAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGa                                                                                                                                                 <  1:632771/128‑1 (MQ=255)
                    cAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCggg                                                                                                                                       >  2:976591/1‑136 (MQ=255)
                    cAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCggg                                                                                                                                       <  1:976591/136‑1 (MQ=255)
                      gTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGt                                                                                                                       <  2:192096/150‑1 (MQ=255)
                                aGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAgggt                                                                                                             <  2:221352/150‑1 (MQ=255)
                                        aGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAgggtcggggtcg                                                                                                     <  1:1383159/150‑1 (MQ=255)
                                         gTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGa                                                                                                    >  2:84028/1‑150 (MQ=255)
                                                 ccAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGt                                                                                            >  1:391240/1‑150 (MQ=255)
                                                       cACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGa                                                                                      >  2:887211/1‑150 (MQ=255)
                                                       cACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGa                                                                                      >  1:118394/1‑150 (MQ=255)
                                                          cGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCTATGTAGTCGAacg                                                                                   >  2:33512/1‑150 (MQ=255)
                                                          cGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGCGGTCGAGGCCATGTAGTCGAacg                                                                                   >  2:277247/1‑150 (MQ=255)
                                                           ggCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGaa                                                                                     <  2:929714/147‑1 (MQ=255)
                                                           ggCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGaa                                                                                     >  1:929714/1‑147 (MQ=255)
                                                           ggCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAacga                                                                                  <  2:208535/150‑1 (MQ=255)
                                                            gCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAacgac                                                                                 <  2:391240/150‑1 (MQ=255)
                                                              gAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGAGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAAcgacga                                                                               >  1:169027/1‑150 (MQ=255)
                                                                       cAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGtt                                                                      <  2:464276/150‑1 (MQ=255)
                                                                          cATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTa                                                                   <  1:229955/150‑1 (MQ=255)
                                                                          cATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTa                                                                   <  1:84028/150‑1 (MQ=255)
                                                                              gCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCTATGTAGTCGAACGACGACAACAGGTTGTAGAAc                                                               <  1:33512/150‑1 (MQ=255)
                                                                              gCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAAc                                                               <  1:58900/150‑1 (MQ=255)
                                                                                    tCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAg                                                                                                                >  2:1417327/1‑95 (MQ=255)
                                                                                    tCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAg                                                                                                                <  1:1417327/95‑1 (MQ=255)
                                                                                    tCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATc                                                         >  2:854769/1‑150 (MQ=255)
                                                                                          aTCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACgg                                                   >  1:503994/1‑150 (MQ=255)
                                                                                           tCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAg                                                                         <  1:121883/127‑1 (MQ=255)
                                                                                           tCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAg                                                                         >  2:121883/1‑127 (MQ=255)
                                                                                           tCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGATGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACggg                                                  <  1:442679/150‑1 (MQ=255)
                                                                                             gcgTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAAc                                                               <  1:65430/135‑1 (MQ=255)
                                                                                             gcgTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAAc                                                               >  2:65430/1‑135 (MQ=255)
                                                                                                   tAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTg                                          >  2:1176838/1‑150 (MQ=255)
                                                                                                     gATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTc                                        >  2:345851/1‑150 (MQ=255)
                                                                                                     gATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGATGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTc                                        >  1:1358045/1‑150 (MQ=255)
                                                                                                      aTGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCt                                       <  1:854769/150‑1 (MQ=255)
                                                                                                           cGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTTTCTTTttc                                  >  2:546501/1‑150 (MQ=255)
                                                                                                                  tcgtcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAgc                           <  1:345851/150‑1 (MQ=255)
                                                                                                                     tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAAc                                                                                                                  >  1:609981/1‑60 (MQ=255)
                                                                                                                     tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAAc                                                                                                                  <  2:609981/60‑1 (MQ=255)
                                                                                                                     tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCttttt                                   <  1:631130/139‑1 (MQ=255)
                                                                                                                     tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCttttt                                   >  2:631130/1‑139 (MQ=255)
                                                                                                                                     ggTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCacac        <  1:1176838/150‑1 (MQ=255)
                                                                                                                                       ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTg                                                                     <  2:1139687/85‑1 (MQ=255)
                                                                                                                                       ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTg                                                                     >  1:1139687/3‑87 (MQ=255)
                                                                                                                                       ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATgg              <  1:1126195/140‑1 (MQ=255)
                                                                                                                                       ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATgg              >  2:1126195/3‑142 (MQ=255)
                                                                                                                                       ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGc      >  1:506267/3‑150 (MQ=255)
                                                                                                                                         gcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCAt    >  2:296328/1‑150 (MQ=255)
                                                                                                                                           agcTGGGCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCATga  >  2:708609/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GTGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCATGA  >  NC_002947/1354127‑1354415

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: