Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 790,960 T→A 6.7% Q46L (CAG→CTG)  PP_0677 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947790,9600TA 6.7% 193.3 / 5.1 60Q46L (CAG→CTG) PP_0677lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base A (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GCTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAA  >  NC_002947/790817‑791105
                                                                                                                                               |                                                                                                                                                 
gcTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTggcggc                                                                                                                                             <  2:66725/150‑1 (MQ=255)
        tATTGGCGATGTCATTCACTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCActgctg                                                                                                                                     >  2:808429/1‑150 (MQ=255)
         aTTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACtgctgc                                                                                                                                    <  2:167410/150‑1 (MQ=255)
          ttGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccgggc                                                                                                                                        >  2:1297561/1‑141 (MQ=255)
          ttGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccgggc                                                                                                                                        <  1:1297561/145‑5 (MQ=255)
             gcgcTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGcc                                                                                                                                >  2:1443394/1‑150 (MQ=255)
                  gTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCt                                                                                                                           <  1:59069/150‑1 (MQ=255)
                      ttcGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGaa                                                                                                                       <  1:756341/150‑1 (MQ=255)
                      ttcGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGaa                                                                                                                       <  2:1002446/150‑1 (MQ=255)
                      ttcGCTGTTCTGCAGGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGaa                                                                                                                       <  2:1345917/150‑1 (MQ=255)
                          cTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGACCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCa                                                                                                                   >  2:971645/1‑150 (MQ=255)
                           tGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCAt                                                                                                                  >  2:1332893/1‑150 (MQ=255)
                           tGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGACCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCAt                                                                                                                  <  1:971645/150‑1 (MQ=255)
                                  cATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATa                                                                                                           <  1:1332893/150‑1 (MQ=255)
                                        gTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCagg                                                                                                     >  2:166667/1‑150 (MQ=255)
                                                   ggAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGt                                                                                          >  1:1145922/1‑150 (MQ=255)
                                                   ggAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGt                                                                                          >  2:992267/1‑150 (MQ=255)
                                                     aaGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTgg                                                                                        <  1:1557311/150‑1 (MQ=255)
                                                      aGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGt                                                                                       <  2:1598843/150‑1 (MQ=255)
                                                         cATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCActg                                                                                                                                        >  1:233084/1‑98 (MQ=255)
                                                         cATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCActg                                                                                                                                        <  2:233084/98‑1 (MQ=255)
                                                         cATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGa                                                                                    >  2:1149908/1‑150 (MQ=255)
                                                               cccAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTgg                                                                                                                         <  1:441667/107‑1 (MQ=255)
                                                               cccAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTgg                                                                                                                         >  2:441667/1‑107 (MQ=255)
                                                                 cAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGc                                                                            <  1:808429/150‑1 (MQ=255)
                                                                  aGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCt                                                                           >  1:1060809/1‑150 (MQ=255)
                                                                     cGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCgt                                                                        >  2:822142/1‑150 (MQ=255)
                                                                            gTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACTATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTgg                                                                                        <  1:1266631/127‑1 (MQ=255)
                                                                            gTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACTATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTgg                                                                                        >  2:1266631/1‑127 (MQ=255)
                                                                             tACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTggcgg                                                                                                                                              <  1:899967/72‑1 (MQ=255)
                                                                             tACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTggcgg                                                                                                                                              >  2:899967/1‑72 (MQ=255)
                                                                                 gACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGAc                                                                                   <  1:745301/127‑1 (MQ=255)
                                                                                 gACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGAc                                                                                   >  2:745301/1‑127 (MQ=255)
                                                                                    gCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATa                                                                                                           >  1:106904/1‑100 (MQ=255)
                                                                                    gCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATa                                                                                                           <  2:106904/100‑1 (MQ=255)
                                                                                    gCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAAc                                                         <  1:1149908/150‑1 (MQ=255)
                                                                                      tCGTGCAGGTAGCGCGGTACGTTCTGGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAAccc                                                       <  1:1156994/150‑1 (MQ=255)
                                                                                          gCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACa                                                   <  1:1443394/150‑1 (MQ=255)
                                                                                           cAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAg                                                  <  1:1100884/150‑1 (MQ=255)
                                                                                                        aCGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATg                                     <  2:1060809/150‑1 (MQ=255)
                                                                                                             cTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGca                                   >  1:815651/1‑147 (MQ=255)
                                                                                                             cTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGca                                   <  2:815651/147‑1 (MQ=255)
                                                                                                              ttGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGca                                   <  2:1219068/146‑1 (MQ=255)
                                                                                                              ttGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGca                                   >  1:1219068/1‑146 (MQ=255)
                                                                                                                    aaGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAAc                         <  1:166667/150‑1 (MQ=255)
                                                                                                                       ctGCCCACCATCAGCCCGGCCGCCTGGCGGCCAATGCTTCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAg                      <  2:793125/149‑1 (MQ=255)
                                                                                                                         gCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTa                    >  1:53297/1‑150 (MQ=255)
                                                                                                                           ccaccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGa                                                     >  2:180477/1‑115 (MQ=255)
                                                                                                                           ccaccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGa                                                     <  1:180477/115‑1 (MQ=255)
                                                                                                                             accaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACcgc                >  1:907839/1‑150 (MQ=255)
                                                                                                                              ccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACcgcg               <  2:1646396/150‑1 (MQ=255)
                                                                                                                               caTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGa              <  1:992267/150‑1 (MQ=255)
                                                                                                                                    gCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGggtggt                                                                 >  2:750066/1‑94 (MQ=255)
                                                                                                                                    gCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGggtggt                                                                 <  1:750066/94‑1 (MQ=255)
                                                                                                                                    gCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAg         >  1:641719/1‑150 (MQ=255)
                                                                                                                                    cagtggCCGCCAGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAg         <  1:1509744/146‑1 (MQ=255)
                                                                                                                                    cagtggCCGCCAGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAg         >  2:1509744/5‑150 (MQ=255)
                                                                                                                                         gccgccTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGt                                                                                                >  1:1556510/1‑58 (MQ=255)
                                                                                                                                         gccgccTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGt                                                                                                <  2:1556510/58‑1 (MQ=255)
                                                                                                                                           cgccTGGCGGCCACTGCTCCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACaaa  >  1:874948/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAA  >  NC_002947/790817‑791105

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: ‑