Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,712,518	(C)6→7	intergenic (+360/‑251)	nth → / → dtpA	DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,712,512	1	.	C	100.0%	214.7 / NA	59	intergenic (+354/‑257)	nth/dtpA	DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (34/25);  total (34/25)

ATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712376‑1712645                                                                                                                                          |                                                                                                                                       aTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑cc                                                                                                                                     <  2:315231/139‑1 (MQ=255)  tCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑ccc                                                                                                                                    <  1:402333/139‑1 (MQ=255)   cAGGTTAACATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑cccc                                                                                                                                   <  2:643144/139‑1 (MQ=255)    aGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑ccccc                                                                                                                                  <  2:357470/139‑1 (MQ=255)                                         ttGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttt                                                                                              >  1:993049/1‑139 (MQ=255)                                         ttGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttt                                                                                              >  1:467035/1‑139 (MQ=255)                                           gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  1:483335/1‑139 (MQ=255)                                           gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  2:737565/1‑139 (MQ=255)                                             ccTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                                                                          >  2:179114/1‑139 (MQ=255)                                             ccTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                                                                          >  1:242101/1‑139 (MQ=255)                                                aCAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaa                                                                                       >  1:884807/1‑139 (MQ=255)                                                         tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                              >  1:871044/1‑139 (MQ=255)                                                         tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                              >  1:666006/1‑139 (MQ=255)                                                         tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                              >  1:501802/1‑139 (MQ=255)                                                          acacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                                                             <  1:176480/139‑1 (MQ=255)                                                             cTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                          <  1:163976/139‑1 (MQ=255)                                                              tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGAt                                                                         >  2:473324/1‑139 (MQ=255)                                                              tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑ccc                                                                                                                                    <  2:536693/79‑1 (MQ=255)                                                              tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑ccc                                                                                                                                    >  1:536693/1‑79 (MQ=255)                                                                gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTaa                                                                                 >  2:544023/1‑129 (MQ=255)                                                                gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTaa                                                                                 <  1:544023/129‑1 (MQ=255)                                                                    tttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAtt                                                                   <  1:179114/139‑1 (MQ=255)                                                                          aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑ccccccctct                                                                                                                             >  2:757769/1‑70 (MQ=255)                                                                          aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑ccccccctct                                                                                                                             <  1:757769/74‑5 (MQ=255)                                                                            aaaCTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           <  2:666006/139‑1 (MQ=255)                                                                                        ctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcc                                               >  2:721302/1‑139 (MQ=255)                                                                                         tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  2:298772/1‑139 (MQ=255)                                                                                         tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:501802/139‑1 (MQ=255)                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  1:75611/139‑1 (MQ=255)                                                                                             tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  1:298772/139‑1 (MQ=255)                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           >  2:195235/1‑120 (MQ=255)                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           <  1:195235/120‑1 (MQ=255)                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:609689/131‑1 (MQ=255)                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  1:609689/1‑131 (MQ=255)                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  1:567193/1‑139 (MQ=255)                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACACTTCCCGGAATGGTTGAATTCCCCGc                                           >  1:1002465/1‑134 (MQ=255)                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACACTTCCCGGAATGGTTGAATTCCCCGc                                           <  2:1002465/134‑1 (MQ=255)                                                                                                   aaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    <  1:163313/139‑1 (MQ=255)                                                                                                     aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  1:129478/1‑130 (MQ=255)                                                                                                     aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  2:129478/130‑1 (MQ=255)                                                                                                         tAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGATTTCCCCGCCTCAGTTATATGt                              >  1:621018/1‑139 (MQ=255)                                                                                                         tAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  1:388453/123‑1 (MQ=255)                                                                                                         tAAATAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  2:388453/1‑123 (MQ=255)                                                                                                          aaaTAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTa                             >  1:775942/1‑139 (MQ=255)                                                                                                           aaTAGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:823330/139‑1 (MQ=255)                                                                                                              aGCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAg                         >  2:96540/1‑139 (MQ=255)                                                                                                               gCAGAACATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGa                        <  1:507270/139‑1 (MQ=255)                                                                                                                   aaCATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  1:485662/1‑118 (MQ=255)                                                                                                                   aaCATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:485662/118‑1 (MQ=255)                                                                                                                     cATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    >  2:57582/1‑121 (MQ=255)                                                                                                                     cATATCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    <  1:57582/121‑1 (MQ=255)                                                                                                                        atCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               >  2:348826/1‑75 (MQ=255)                                                                                                                        atCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  1:348826/75‑1 (MQ=255)                                                                                                                        atCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            >  1:1030595/1‑126 (MQ=255)                                                                                                                        atCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  2:1030595/126‑1 (MQ=255)                                                                                                                        atCGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaa               >  2:989675/1‑139 (MQ=255)                                                                                                                          cGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  2:389325/73‑1 (MQ=255)                                                                                                                          cGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               >  1:389325/1‑73 (MQ=255)                                                                                                                          cGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           <  2:342599/93‑1 (MQ=255)                                                                                                                          cGCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           >  1:342599/1‑93 (MQ=255)                                                                                                                           gCCGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAgg            <  2:884807/139‑1 (MQ=255)                                                                                                                             cGTTAAGCCATCACCCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAACGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATCATTACAAAGGa           >  1:120337/1‑139 (MQ=255)                                                                                                                             cGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAAt                                                                                                             <  1:763835/40‑1 (MQ=255)                                                                                                                             cGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAAt                                                                                                             >  2:763835/1‑40 (MQ=255)                                                                                                                             cGTTAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc          >  1:38801/1‑139 (MQ=255)                                                                                                                               ttAAGCCATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACtt        >  1:556000/1‑139 (MQ=255)                                                                                                                                     cATCAC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGa  <  2:993049/139‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       ATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712376‑1712645

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: ‑