Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,685,542 (C)6→5 coding (1047/1404 nt) fumC ← fumarase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,5370C.100.0% 254.5 / NA 57coding (1052/1404 nt)fumCfumarase C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (29/28);  total (29/28)

CGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685399‑1685668
                                                                                                                                          |                                                                                                                                   
cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                    <  1:267706/140‑1 (MQ=255)
   tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccc                                                                                                                                 >  1:1085626/1‑140 (MQ=255)
   tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccc                                                                                                                                 >  1:1014583/1‑140 (MQ=255)
         ttCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTg                                                                                                                          <  2:483026/140‑1 (MQ=255)
                  aataccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGGCGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                                 >  1:275456/5‑140 (MQ=255)
                  aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                                 <  2:293112/140‑1 (MQ=255)
                  aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                                 <  2:732049/140‑1 (MQ=255)
                        aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTc                                                                                                           >  1:682609/1‑140 (MQ=255)
                        aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTc                                                                                                           >  2:153972/1‑140 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         >  1:1593727/1‑140 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         >  2:807756/1‑140 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         >  1:664895/1‑140 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         >  2:861193/1‑140 (MQ=255)
                            cgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                       >  1:247897/1‑140 (MQ=255)
                            cgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                       >  2:1077543/1‑140 (MQ=255)
                               aGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                    <  1:1030049/140‑1 (MQ=255)
                               aGTGGTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                    <  2:517920/140‑1 (MQ=255)
                                 tgtttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCAcc                                                                                                  <  1:366251/140‑1 (MQ=255)
                                    ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                               <  2:543885/140‑1 (MQ=255)
                                    ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATATTGATCGCAACGTCGTTCCCCATCACCTGa                                                                                               <  2:891501/140‑1 (MQ=255)
                                         aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTTAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                                                                          <  1:1182894/140‑1 (MQ=255)
                                               ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                    <  2:1188319/140‑1 (MQ=255)
                                                tCCATGCTATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                                                                   <  2:103313/140‑1 (MQ=255)
                                                      ccATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                                                             <  1:1258830/140‑1 (MQ=255)
                                                         tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                          <  2:468482/140‑1 (MQ=255)
                                                          cTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                         <  2:667952/140‑1 (MQ=255)
                                                           tGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                        <  2:833990/140‑1 (MQ=255)
                                                             ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                      <  2:677631/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:343794/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  2:1014583/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:693491/140‑1 (MQ=255)
                                                                     cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                              >  1:1557558/1‑140 (MQ=255)
                                                                       cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                            >  2:1591775/1‑140 (MQ=255)
                                                                        aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTggg                                                           <  2:945692/140‑1 (MQ=255)
                                                                            aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCa                                                       <  1:394702/140‑1 (MQ=255)
                                                                              tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACt                                                     >  1:272835/1‑140 (MQ=255)
                                                                                      aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                             >  2:1402601/1‑140 (MQ=255)
                                                                                      aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                             <  2:1557558/140‑1 (MQ=255)
                                                                                      aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgc                                             >  2:941807/1‑139 (MQ=255)
                                                                                        tgtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCa                                           >  2:99367/1‑140 (MQ=255)
                                                                                         gtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCAt                                          <  2:1160977/140‑1 (MQ=255)
                                                                                           ggATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        <  2:1038752/140‑1 (MQ=255)
                                                                                               cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGa                                    >  2:342606/1‑140 (MQ=255)
                                                                                                     tGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                              <  1:443250/140‑1 (MQ=255)
                                                                                                       gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCg                            >  2:1208438/1‑140 (MQ=255)
                                                                                                         cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGc                          >  2:1273259/1‑140 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATGGCCACGCCGTTCCCAACCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGGGTTGGGTTCCCTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       >  2:1473951/1‑140 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       >  2:1136546/1‑140 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       >  1:148048/1‑140 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       <  1:66150/140‑1 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       >  2:348936/1‑140 (MQ=255)
                                                                                                              gACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCAtt                     >  2:725546/1‑140 (MQ=255)
                                                                                                                  ttcagttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTcc                 >  2:456680/1‑140 (MQ=255)
                                                                                                                   tcagttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                >  2:752055/1‑140 (MQ=255)
                                                                                                                   tcagttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                <  2:803500/140‑1 (MQ=255)
                                                                                                                   tcagttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                <  2:272835/140‑1 (MQ=255)
                                                                                                                      gttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGa             >  1:794131/1‑140 (MQ=255)
                                                                                                                       ttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAt            >  2:630127/1‑140 (MQ=255)
                                                                                                                        tcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTGCCCCATCGCCTGGCAGCAGAGCATGGTTAATGCCTCGCACTGTGTTGGGTTCACTTCCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAtt           <  2:1389642/140‑1 (MQ=255)
                                                                                                                          aaaGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga         <  2:74536/140‑1 (MQ=255)
                                                                                                                                 ccGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCa  >  1:570644/1‑140 (MQ=255)
                                                                                                                                 ccGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCa  >  1:1492819/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                   
CGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685399‑1685668

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