Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,515,726	T→C	Y63C (TAT→TGT)	dam ←	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,515,726	0	T	C	100.0%	266.3 / NA	73	Y63C (TAT→TGT)	dam	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/40);  total (33/40)

CGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCA  >  NC_000913/3515590‑3515860                                                                                                                                         |                                                                                                                                       cGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACag                                                                                                                                      <  2:524679/139‑1 (MQ=255)  gCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACaga                                                                                                                                     <  1:187281/139‑1 (MQ=255)  gCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACaga                                                                                                                                     >  2:1276675/1‑139 (MQ=255)      gACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACt                                                                                                                                 >  2:2047542/1‑139 (MQ=255)           aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCa                                                                                                                            >  2:588510/1‑139 (MQ=255)           aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCa                                                                                                                            <  2:617553/139‑1 (MQ=255)              ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         <  2:642994/139‑1 (MQ=255)              ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         <  1:1904591/139‑1 (MQ=255)              ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         <  2:146958/139‑1 (MQ=255)                  ccTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCa                                                                                                                            >  1:739044/1‑132 (MQ=255)                  ccTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCa                                                                                                                            <  2:739044/132‑1 (MQ=255)                     ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTAtt                                                                                                                  >  2:1569368/1‑139 (MQ=255)                          ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         >  1:423153/1‑127 (MQ=255)                          ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         <  2:423153/127‑1 (MQ=255)                          ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGatat                                                                                                             >  1:1525613/1‑139 (MQ=255)                              tgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                         <  2:1525613/139‑1 (MQ=255)                              tgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                         <  1:1276675/139‑1 (MQ=255)                               gAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                        <  1:588510/139‑1 (MQ=255)                                     ttCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTa                                                                                                                    <  1:1696863/121‑1 (MQ=255)                                     ttCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTa                                                                                                                    >  2:1696863/1‑121 (MQ=255)                                       cgcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGAtata                                                                                                >  2:907558/1‑139 (MQ=255)                                                 aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCg                                                                                                                       >  2:475154/1‑106 (MQ=255)                                                 aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCg                                                                                                                       <  1:475154/106‑1 (MQ=255)                                                 aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAg                                                                                      >  1:706674/1‑139 (MQ=255)                                                  tAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                     <  1:907558/139‑1 (MQ=255)                                                    gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCg                                                                                   <  1:1988047/139‑1 (MQ=255)                                                      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                              >  2:1733919/1‑110 (MQ=255)                                                      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                              <  1:1733919/110‑1 (MQ=255)                                                      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                        >  1:1421316/1‑116 (MQ=255)                                                      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                        <  2:1421316/116‑1 (MQ=255)                                                      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgt                                                                                 >  2:2051143/1‑139 (MQ=255)                                                       aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                                <  2:604893/139‑1 (MQ=255)                                                       aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                                <  1:1569368/139‑1 (MQ=255)                                                        aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCt                                                                                                                     >  1:2028106/1‑101 (MQ=255)                                                        aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCt                                                                                                                     <  2:2028106/101‑1 (MQ=255)                                                        aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTGTTGATATCGGAAAGGATATAACGAGAAAAGTCGgtgt                                                                               >  1:1309340/1‑139 (MQ=255)                                                            cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         <  2:538283/93‑1 (MQ=255)                                                            cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         >  1:538283/1‑93 (MQ=255)                                                            cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                        >  2:756260/1‑110 (MQ=255)                                                            cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                        <  1:756260/110‑1 (MQ=255)                                                            cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                                <  1:1955261/134‑1 (MQ=255)                                                            cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                                >  2:1955261/1‑134 (MQ=255)                                                             ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                              <  1:569777/103‑1 (MQ=255)                                                             ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                              >  2:569777/1‑103 (MQ=255)                                                                        tttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                               <  2:431215/139‑1 (MQ=255)                                                                          tCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                                                             <  1:2047542/139‑1 (MQ=255)                                                                              ggacaaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCg                                                                                                                       <  2:1240082/75‑1 (MQ=255)                                                                              gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCg                                                                                                                       >  1:1240082/1‑77 (MQ=255)                                                                                   aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaa                                                    <  1:755112/139‑1 (MQ=255)                                                                                       gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGAt                                                                                                                              >  2:1266198/1‑61 (MQ=255)                                                                                       gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGAt                                                                                                                              <  1:1266198/61‑1 (MQ=255)                                                                                       gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                >  2:754635/1‑139 (MQ=255)                                                                                        cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                      >  1:100404/1‑116 (MQ=255)                                                                                        cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                      <  2:100404/116‑1 (MQ=255)                                                                                           gcgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         >  1:90149/1‑62 (MQ=255)                                                                                           gcgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                         <  2:90149/62‑1 (MQ=255)                                                                                                cGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAg                                       >  1:1417013/1‑139 (MQ=255)                                                                                                       tacgtacTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                >  1:988217/1‑139 (MQ=255)                                                                                                            acTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                               >  1:63345/1‑55 (MQ=255)                                                                                                            acTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                               <  2:63345/55‑1 (MQ=255)                                                                                                                 tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCgg                                                                                                          >  2:542870/1‑55 (MQ=255)                                                                                                                 tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCgg                                                                                                          <  1:542870/55‑1 (MQ=255)                                                                                                                 tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                         <  1:343954/56‑1 (MQ=255)                                                                                                                 tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                         >  2:343954/1‑56 (MQ=255)                                                                                                                 tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                      >  2:566503/1‑139 (MQ=255)                                                                                                                   aGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaaa                    <  1:2051143/139‑1 (MQ=255)                                                                                                                        gCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCg               <  1:566503/139‑1 (MQ=255)                                                                                                                             ttCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTaa          <  1:754635/139‑1 (MQ=255)                                                                                                                               cacaATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAacac                                                                     >  2:651443/1‑78 (MQ=255)                                                                                                                               cacaATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAacac                                                                     <  1:651443/78‑1 (MQ=255)                                                                                                                                    tGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg                  <  1:21847/124‑1 (MQ=255)                                                                                                                                    tGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg                  >  2:21847/1‑124 (MQ=255)                                                                                                                                     gTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca  <  2:1417013/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCA  >  NC_000913/3515590‑3515860

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