Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,015,773 C→A 8.5% A5S (GCT→TCT)  PP_5679 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,015,7730CA 8.5% 169.5 / 2.8 57A5S (GCT→TCT) PP_5679hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base A (3/2);  total (31/26)
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.23e-01

GACGAGGATTCCGTTGGCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTT  >  NC_002947/5015624‑5015922
                                                                                                                                                     |                                                                                                                                                     
gacgaGGATTCCGTTGGCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGc                                                                                                                                                       <  2:227420/150‑1 (MQ=255)
gacgaGGATTCCGTTGGCATACTTGGCCTCATTAGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGc                                                                                                                                                       >  2:435966/1‑150 (MQ=255)
        ttCCGTTGGCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCgg                                                                                                                                               >  1:170381/1‑150 (MQ=255)
         tCCGTTGGCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGt                                                                                                                                              <  1:1059057/150‑1 (MQ=255)
                gCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCa                                                                                                                                                     >  1:523711/1‑136 (MQ=255)
                gCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCa                                                                                                                                                     <  2:523711/136‑1 (MQ=255)
                gCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTa                                                                                                                                       <  1:786529/150‑1 (MQ=255)
                              aTTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGa                                                                                                                         <  2:199545/150‑1 (MQ=255)
                               tttGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAg                                                                                                                        <  2:170381/150‑1 (MQ=255)
                                   ggTTGTACATGCCCGTGCGGATACTCGCCCGTAGCCGCTTGCAGAGTGACTGCTGGGGACGAGGCTTCTCGCCAGCCCTGGCCGTATATGTAGCCATTGCTGCTTTGGTCGGTGCCAGAGCGGTCATCCTGGCGCTCTCGTGTGAGTCCa                                                                                                                    <  2:351593/150‑1 (MQ=255)
                                    gTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCAc                                                                                                                   >  2:1247825/1‑150 (MQ=255)
                                    gTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCCTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCAc                                                                                                                   >  2:928468/1‑150 (MQ=255)
                                     ttGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGt                                                                                                                                              <  1:630960/122‑1 (MQ=255)
                                     ttGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGt                                                                                                                                              >  2:630960/1‑122 (MQ=255)
                                      tGAACTGGCCCGTGCGGTTACGCACCTGTAGCCGCTTGCAGAGTGGCTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTccggta                                                                                                                 <  1:685475/150‑5 (MQ=255)
                                       gAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATc                                                                                                                                          >  1:1144623/1‑124 (MQ=255)
                                       gAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATc                                                                                                                                          <  2:1144623/124‑1 (MQ=255)
                                           tGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCAGGCTGGGTCATCCTAGCGCTCTCGTTTGAGTGCACgatgaga                                                                                                            >  1:976744/1‑150 (MQ=255)
                                             gCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACgatgagatg                                                                                                          >  2:688921/1‑150 (MQ=255)
                                               ccGTGCGGATCCTCCCCCGTAGCCCCTTGCACATCGACTGCTGGGCACTAAGCTTTTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGaa                                                                                                        <  2:976744/150‑1 (MQ=255)
                                               ccGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGCATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGaa                                                                                                        >  2:1257583/1‑150 (MQ=255)
                                                cGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAAt                                                                                                       <  1:435966/150‑1 (MQ=255)
                                                 gTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATa                                                                                                      >  1:1100933/1‑150 (MQ=255)
                                                     ggATACTCACCCCTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGcc                                                                                                  >  1:29363/1‑150 (MQ=255)
                                                             aCCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGCTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCgg                                                                                          >  1:561524/1‑150 (MQ=255)
                                                             aCCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCa                                                                                                                                            >  2:1001557/1‑100 (MQ=255)
                                                             aCCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCa                                                                                                                                            <  1:1001557/100‑1 (MQ=255)
                                                                       gCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACgat                                                                                                                <  2:497740/118‑1 (MQ=255)
                                                                       gCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACgat                                                                                                                >  1:497740/1‑118 (MQ=255)
                                                                         ttGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGtc                                                                                        <  2:446350/140‑1 (MQ=255)
                                                                         ttGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGtc                                                                                        >  1:446350/1‑140 (MQ=255)
                                                                               agTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGc                                                                                               >  1:630619/1‑127 (MQ=255)
                                                                               agTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGc                                                                                               <  2:630619/127‑1 (MQ=255)
                                                                                  gACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCCTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACt                                                                     <  1:928468/150‑1 (MQ=255)
                                                                                        tgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGACAGCt                                                                                                                                                  <  2:1035781/67‑1 (MQ=255)
                                                                                        tgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGACAGCt                                                                                                                                                  >  1:1035781/1‑67 (MQ=255)
                                                                                        tgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGACAGCt                                                                                                                                                  >  1:1036222/1‑67 (MQ=255)
                                                                                        tgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGACAGCt                                                                                                                                                  <  2:1036222/67‑1 (MQ=255)
                                                                                                aaGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGa                                                       <  1:688921/150‑1 (MQ=255)
                                                                                                    ttCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATaa                                                                                                     >  2:966054/1‑100 (MQ=255)
                                                                                                    ttCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATaa                                                                                                     <  1:966054/100‑1 (MQ=255)
                                                                                                             gCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGAccc                                          >  1:1229347/1‑150 (MQ=255)
                                                                                                              cccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCg                                         >  2:334799/1‑150 (MQ=255)
                                                                                                                cTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCg                                                                <  1:590753/125‑1 (MQ=255)
                                                                                                                cTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCg                                                                >  2:590753/1‑125 (MQ=255)
                                                                                                                cTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaa                                       >  1:680807/1‑150 (MQ=255)
                                                                                                                 tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                      >  2:723448/1‑150 (MQ=255)
                                                                                                                        ataGGTAGCCACTGCTTCCTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTc                               >  2:259539/1‑150 (MQ=255)
                                                                                                                           ggTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCtt                            <  2:29363/150‑1 (MQ=255)
                                                                                                                            gTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtg                           <  2:1229347/150‑1 (MQ=255)
                                                                                                                            gTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtg                           >  2:144975/1‑150 (MQ=255)
                                                                                                                                ccACTGCTTCTTTGCTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGCGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGaa                       <  1:1226797/150‑1 (MQ=255)
                                                                                                                                ccACTGCTTCCTTGATCGGCGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGaa                       <  1:259539/150‑1 (MQ=255)
                                                                                                                                 cACTGCTTCTTTGATCGGAGCCAGCTCGGTCCTCCTCGTGATCTCGTATGCGTTCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTGTCTGACGCGAAAGTACGTCCTTGTGAAc                      >  1:543416/1‑150 (MQ=255)
                                                                                                                                  aCTGCTTCTTTGATCGGAGACAGCTCTGTCAGCCGAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACt                     >  1:670484/1‑150 (MQ=255)
                                                                                                                                     gCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgat                  <  2:1100933/150‑1 (MQ=255)
                                                                                                                                          tttGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAg             >  1:1131074/1‑150 (MQ=255)
                                                                                                                                              aTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTc         >  2:614087/1‑150 (MQ=255)
                                                                                                                                                  gAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTgg     >  2:953808/1‑150 (MQ=255)
                                                                                                                                                     ccAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCtt  <  2:1131074/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACGAGGATTCCGTTGGCATACTTGGCCTCATTTGGGTTGAACTGGCCCGTGCGGATACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTT  >  NC_002947/5015624‑5015922

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