Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,345,084	A→G	Y123C (TAC→TGC)	rpoN →	RNA polymerase, sigma 54 (sigma N) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,345,084	0	A	G	100.0%	232.4 / NA	64	Y123C (TAC→TGC)	rpoN	RNA polymerase, sigma 54 (sigma N) factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/26);  total (38/26)

ATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGAAAGTATAGGCGATG  >  NC_000913/3344948‑3345218                                                                                                                                         |                                                                                                                                       aTGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCt                                                                                                                                     >  2:119958/1‑140 (MQ=255)  tGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTg                                                                                                                                    <  1:33437/140‑1 (MQ=255)             cTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAggg                                                                                                                         >  2:646723/1‑139 (MQ=255)                  gCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGc                                                                                                                    >  1:274975/1‑140 (MQ=255)                          ccAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCg                                                                                                            <  1:666652/140‑1 (MQ=255)                          ccAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCg                                                                                                            <  2:274975/140‑1 (MQ=255)                                  gACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGa                                                                                                    >  1:209559/1‑140 (MQ=255)                                         tttACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGAc                                                                                             >  2:621905/1‑140 (MQ=255)                                         tttACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGAc                                                                                             >  2:235886/1‑140 (MQ=255)                                         tttACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGAc                                                                                             <  2:485311/140‑1 (MQ=255)                                                 gCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTAt                                                                                     >  2:565393/1‑140 (MQ=255)                                                 gCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTAt                                                                                     >  1:506341/1‑140 (MQ=255)                                                 gCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTAt                                                                                     >  1:954662/1‑140 (MQ=255)                                                  cTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTAtt                                                                                    <  1:49968/140‑1 (MQ=255)                                                    ggTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGc                                                                                  <  2:954662/140‑1 (MQ=255)                                                     gTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCt                                                                                 <  1:565393/140‑1 (MQ=255)                                                            aTCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTa                                                                          >  1:750156/1‑140 (MQ=255)                                                            aTCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTa                                                                          >  2:461006/1‑140 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGCCACTGACCGCGCTATTGCTACCTCTAt                                                                         >  1:719712/1‑140 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACcgcg                                                                                         <  1:62441/124‑1 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACcgcg                                                                                         >  2:62441/1‑124 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACcgc                                                                                          <  1:653906/123‑1 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACcgc                                                                                          >  2:653906/1‑123 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACcg                                                                                           <  2:83460/122‑1 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACcg                                                                                           >  1:83460/1‑122 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTAt                                                                         >  1:530168/1‑140 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTAt                                                                         >  1:67489/1‑140 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTAt                                                                         >  2:1074018/1‑140 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTAt                                                                         >  1:1077936/1‑140 (MQ=255)                                                             tCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTAt                                                                         >  2:844749/1‑140 (MQ=255)                                                                ggCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTAtcgt                                                                      >  1:920431/1‑140 (MQ=255)                                                                ggCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTAtcgt                                                                      >  1:790609/1‑140 (MQ=255)                                                                   aCCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGa                                                                   <  1:846697/140‑1 (MQ=255)                                                                     cAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATg                                                                 <  1:621905/140‑1 (MQ=255)                                                                     cAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATg                                                                 >  2:15696/1‑140 (MQ=255)                                                                        cGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCg                                                              >  1:192610/1‑140 (MQ=255)                                                                             aCTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTtgat                                                         >  1:194415/1‑140 (MQ=255)                                                                                 cATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGaaa                                                     <  2:790609/140‑1 (MQ=255)                                                                                  aTTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTTCTACCTCTATCGTCGATGCCGTTGATGAAAc                                                    >  1:325063/1‑140 (MQ=255)                                                                                    tGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCg                                                  <  2:1077936/140‑1 (MQ=255)                                                                                     gacgacGAGCTGCCGGTCTACCAGGGCGAAACGACGCGGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCGACCTCTATCGTCGATGCCGTTGATGAAACCgg                                                 >  2:957842/1‑140 (MQ=255)                                                                                     gacgacGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGACAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAAcc                                                   >  1:650040/1‑138 (MQ=255)                                                                                     gacgacGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGACAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAAcc                                                   <  2:650040/138‑1 (MQ=255)                                                                                               tGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACt                                       <  1:235886/140‑1 (MQ=255)                                                                                                   ggTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTcc                                   <  1:15696/140‑1 (MQ=255)                                                                                                     tCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCg                                 <  1:461006/140‑1 (MQ=255)                                                                                                        aCCAGGGCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGGTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTg                              >  2:843724/1‑140 (MQ=255)                                                                                                              gCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGat                        >  1:263561/1‑140 (MQ=255)                                                                                                              gCGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGat                        >  1:655979/1‑140 (MQ=255)                                                                                                               cGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGata                       >  2:569722/1‑140 (MQ=255)                                                                                                               cGAAACGACGCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGata                       >  2:587739/1‑140 (MQ=255)                                                                                                                   acgacgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATtct                   >  1:913530/1‑140 (MQ=255)                                                                                                                   acgacgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATtct                   >  2:611696/1‑140 (MQ=255)                                                                                                                     gacgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCg                 <  1:986116/140‑1 (MQ=255)                                                                                                                     gacgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCg                 <  1:611696/140‑1 (MQ=255)                                                                                                                     gacgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCg                 <  2:913530/140‑1 (MQ=255)                                                                                                                      acgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGa                <  2:209559/140‑1 (MQ=255)                                                                                                                       cgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGaa               <  2:67489/140‑1 (MQ=255)                                                                                                                       cgCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGaa               >  2:38279/1‑140 (MQ=255)                                                                                                                        gCAGACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGaaa              >  2:713601/1‑140 (MQ=255)                                                                                                                           gACCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGAAAGta           <  2:263561/140‑1 (MQ=255)                                                                                                                            aCCTTGCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGAAAGtat          <  1:613061/140‑1 (MQ=255)                                                                                                                                 gCAGGATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGAAAGTATAGGCg     >  2:435016/1‑140 (MQ=255)                                                                                                                                    ggATTGCCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGAAAGTATAGGCGATg  <  2:194415/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       ATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTACCTCTATCGTCGATGCCGTTGATGAAACCGGTTATCTGACTGTCCCGCTGGAAGATATTCTCGAAAGTATAGGCGATG  >  NC_000913/3344948‑3345218

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