Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,232,752 G→A A161V (GCG→GTG)  lysC ← lysine‑sensitive aspartokinase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,232,7520GA100.0% 227.9 / NA 69A161V (GCG→GTG) lysClysine‑sensitive aspartokinase 3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/38);  total (31/38)

GCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGC  >  NC_000913/4232617‑4232881
                                                                                                                                       |                                                                                                                                 
gCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACgg                                                                                                                                <  2:168002/139‑1 (MQ=255)
       tataATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATc                                                                                                                         >  2:1662995/1‑139 (MQ=255)
               cTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTg                                                                                                                 <  1:1119353/139‑1 (MQ=255)
                tGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGc                                                                                                                <  1:415792/139‑1 (MQ=255)
                 gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCa                                                                                                               >  2:1628556/1‑139 (MQ=255)
                 gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCa                                                                                                               >  2:1615652/1‑139 (MQ=255)
                 gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCa                                                                                                               <  1:1615652/139‑1 (MQ=255)
                 gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCa                                                                                                               <  1:1662995/139‑1 (MQ=255)
                         ggCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCaaa                                                                                                       <  1:25969/139‑1 (MQ=255)
                                     ttGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGt                                                                                               <  2:17404/135‑1 (MQ=255)
                                     ttGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGt                                                                                               >  1:17404/1‑135 (MQ=255)
                                      tGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTa                                                                                          <  2:1621908/139‑1 (MQ=255)
                                        tACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACg                                                                                        <  1:1628556/139‑1 (MQ=255)
                                         aCGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGc                                                                                       <  2:1711136/139‑1 (MQ=255)
                                          cGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCa                                                                                      >  1:343076/1‑139 (MQ=255)
                                           gACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCAt                                                                                     <  2:1124805/139‑1 (MQ=255)
                                             ccTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCa                                                                                   >  2:570090/1‑139 (MQ=255)
                                                ttATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACtt                                                                                >  2:214366/1‑139 (MQ=255)
                                                 tATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCa                                                                                      <  2:931731/132‑1 (MQ=255)
                                                 tATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCa                                                                                      >  1:931731/1‑132 (MQ=255)
                                                    tttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACg                                                                            <  1:805324/139‑1 (MQ=255)
                                                    tttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACg                                                                            <  1:1677154/139‑1 (MQ=255)
                                                     ttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGt                                                                           >  2:554913/1‑139 (MQ=255)
                                                       cGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTAc                                                                         >  2:1482923/1‑139 (MQ=255)
                                                                 aaaTCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCAc                                                               <  2:343076/139‑1 (MQ=255)
                                                                 aaaTCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCAc                                                               <  1:214366/139‑1 (MQ=255)
                                                                      ccTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGc                                                          >  2:1695789/1‑139 (MQ=255)
                                                                       cTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACg                                                                                        >  2:669458/1‑108 (MQ=255)
                                                                       cTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACg                                                                                        <  1:669458/108‑1 (MQ=255)
                                                                          ggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGa                                                      >  2:221898/1‑139 (MQ=255)
                                                                              aTCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTgg                                                                                            <  1:696749/97‑1 (MQ=255)
                                                                              aTCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTgg                                                                                            >  2:696749/1‑97 (MQ=255)
                                                                              aTCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACAt                                                  >  2:881224/1‑139 (MQ=255)
                                                                                   tAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGt                                             <  1:167719/139‑1 (MQ=255)
                                                                                    aaGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGtt                                            <  1:221898/139‑1 (MQ=255)
                                                                                        cTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTcgcg                                        <  1:918506/139‑1 (MQ=255)
                                                                                         ttCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCgcgc                                       <  1:570090/139‑1 (MQ=255)
                                                                                          tCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGt                                                                                                  <  2:1419814/79‑1 (MQ=255)
                                                                                          tCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGt                                                                                                  >  1:1419814/1‑79 (MQ=255)
                                                                                           cATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAg                                     <  1:1482923/139‑1 (MQ=255)
                                                                                           cATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAg                                     >  1:1082535/1‑139 (MQ=255)
                                                                                             ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGa                                   >  1:529199/1‑139 (MQ=255)
                                                                                                  aCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCgcgc                                       >  2:133532/1‑130 (MQ=255)
                                                                                                  aCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCgcgc                                       <  1:133532/130‑1 (MQ=255)
                                                                                                     tGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCAt                                                                                     >  2:980376/1‑81 (MQ=255)
                                                                                                     tGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCAt                                                                                     <  1:980376/81‑1 (MQ=255)
                                                                                                     tGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTcg                                          <  1:528871/124‑1 (MQ=255)
                                                                                                     tGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTcg                                          >  2:528871/1‑124 (MQ=255)
                                                                                                      gggAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCGAACCACTGTGCCTGAACATCGCGTTCGCGCa                                      <  1:565752/127‑1 (MQ=255)
                                                                                                      gggAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCa                                      >  2:565752/1‑127 (MQ=255)
                                                                                                            agcTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTca                              <  1:83135/129‑1 (MQ=255)
                                                                                                            agcTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTca                              >  2:83135/1‑129 (MQ=255)
                                                                                                              cTGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAgg                  <  1:554913/139‑1 (MQ=255)
                                                                                                               tGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAAt                                                                                                      >  2:78384/1‑54 (MQ=255)
                                                                                                               tGCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAAt                                                                                                      <  1:78384/54‑1 (MQ=255)
                                                                                                                gCAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGt                <  2:1082535/139‑1 (MQ=255)
                                                                                                                 cAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAg                                     >  2:536814/1‑117 (MQ=255)
                                                                                                                 cAGCGCGGCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAg                                     <  1:536814/117‑1 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAAt                                                                                                      <  2:280360/45‑1 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAAt                                                                                                      >  1:280360/1‑45 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCAc                                                               <  1:1616144/84‑1 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCAc                                                               >  2:1616144/1‑84 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCgcgc                                       <  1:636606/108‑1 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCgcgc                                       >  2:636606/1‑108 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCa        >  1:1172858/1‑139 (MQ=255)
                                                                                                                        gCCAGTTCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCa        >  2:1522210/1‑139 (MQ=255)
                                                                                                                              tCCGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGc  <  2:1172858/139‑1 (MQ=255)
                                                                                                                                cGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtg                                                             <  2:229983/78‑1 (MQ=255)
                                                                                                                                cGCCAGCACGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtg                                                             >  1:229983/1‑78 (MQ=255)
                                                                                                                                       |                                                                                                                                 
GCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGC  >  NC_000913/4232617‑4232881

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