Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,551,573	A→G	N79S (AAC→AGC)	fdnI →	formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,551,573	0	A	G	100.0%	296.7 / NA	81	N79S (AAC→AGC)	fdnI	formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/43);  total (38/43)

CGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAACATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGAT  >  NC_000913/1551438‑1551708                                                                                                                                        |                                                                                                                                        cGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCAt                                                                                                                                      >  2:995568/1‑139 (MQ=255)  gTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATATTCGTCGCACTGATGTTTATGTTTGTGCGTGTTGTGCATCACAGCATc                                                                                                                                     <  1:497438/139‑1 (MQ=255)        cttcCCGACGCTGCAGTGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGAt                                                                                                                               >  2:1132827/1‑139 (MQ=255)        cttcCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGAt                                                                                                                               <  2:106762/139‑1 (MQ=255)               aCGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATaagaaag                                                                                                                        <  2:780657/139‑1 (MQ=255)                  cTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGAATCCCGGATAAGAAAGata                                                                                                                     >  2:1468069/1‑139 (MQ=255)                   tGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGatat                                                                                                                    >  1:522491/1‑139 (MQ=255)                        tGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               <  1:417952/139‑1 (MQ=255)                        tGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCACTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTTCGTTTCGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               <  1:1132827/139‑1 (MQ=255)                           cTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATaa                                                                                                                             <  1:261121/122‑1 (MQ=255)                           cTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATaa                                                                                                                             >  2:261121/1‑122 (MQ=255)                             gACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               <  2:1181343/134‑1 (MQ=255)                             gACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               >  1:1181343/1‑134 (MQ=255)                             gACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTg                                                                                                          >  2:277826/1‑139 (MQ=255)                             gACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATACCGTGGCTg                                                                                                          <  2:522491/139‑1 (MQ=255)                               cGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATaaga                                                                                                                           >  1:1673720/1‑120 (MQ=255)                               cGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATaaga                                                                                                                           <  2:1673720/120‑1 (MQ=255)                               cGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             <  2:929346/134‑1 (MQ=255)                               cGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             >  1:929346/1‑134 (MQ=255)                               cGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGtt                                                                                                        <  1:137577/139‑1 (MQ=255)                                    aCCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGata                                                                                                                     <  2:948903/121‑1 (MQ=255)                                    aCCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGata                                                                                                                     >  1:948903/1‑121 (MQ=255)                                         cGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATaaga                                                                                                                           <  2:1679953/110‑1 (MQ=255)                                         cGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATaaga                                                                                                                           >  1:1679953/1‑110 (MQ=255)                                          ggTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCg                                                                                             <  1:552009/139‑1 (MQ=255)                                                 cgcAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTAtt                                                                                      >  1:552492/1‑139 (MQ=255)                                                   cAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCg                                                                                             >  2:1299625/1‑130 (MQ=255)                                                   cAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCg                                                                                             <  1:1299625/130‑1 (MQ=255)                                                     gATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGaaa                                                                                  >  1:194329/1‑139 (MQ=255)                                                         ggACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTg                                                                                                          >  1:862494/1‑111 (MQ=255)                                                         ggACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTg                                                                                                          <  2:862494/111‑1 (MQ=255)                                                           aCGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGAt                                                                                                                               <  2:1278520/88‑1 (MQ=255)                                                           aCGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGAt                                                                                                                               >  1:1278520/1‑88 (MQ=255)                                                           aCGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAAt                                                                            <  1:1424103/139‑1 (MQ=255)                                                                         cGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             <  1:905745/92‑1 (MQ=255)                                                                         cGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             >  2:905745/1‑92 (MQ=255)                                                                         cGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAAt                                                                            <  1:698790/125‑1 (MQ=255)                                                                         cGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAAt                                                                            >  2:698790/1‑125 (MQ=255)                                                                           ttcttcGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCgg                                                            >  2:1585807/1‑139 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGAt                                                                                                                               <  2:976504/61‑1 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGAt                                                                                                                               >  1:976504/1‑61 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               >  2:1633117/1‑77 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               <  1:1633117/77‑1 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             <  1:571947/79‑1 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             >  2:571947/1‑79 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGt                                                                                               <  2:1510982/93‑1 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGt                                                                                               >  1:1510982/1‑93 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCa                                                                              >  1:1109928/1‑110 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCa                                                                              <  2:1109928/110‑1 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGt                                                    <  2:930049/136‑1 (MQ=255)                                                                                      tGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGt                                                    >  1:930049/1‑136 (MQ=255)                                                                                             ttCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACg                                          <  2:1006500/139‑1 (MQ=255)                                                                                                             ttatgtttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAgatgat                          <  1:277826/139‑1 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               <  1:1286595/48‑1 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGt                                                                                                               >  2:1286595/1‑48 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             <  2:1383308/50‑1 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             >  1:1383308/1‑50 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAAt                                                                            >  2:777888/1‑83 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAAt                                                                            <  1:777888/83‑1 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGtt                       <  1:618591/136‑1 (MQ=255)                                                                                                                   ttGTGCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGtt                       >  2:618591/1‑136 (MQ=255)                                                                                                                    tgtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCg                                                                                             <  2:344816/65‑1 (MQ=255)                                                                                                                    tgtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCg                                                                                             >  1:344816/1‑65 (MQ=255)                                                                                                                     gtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGt                                                    >  1:73407/1‑105 (MQ=255)                                                                                                                     gtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGt                                                    <  2:73407/105‑1 (MQ=255)                                                                                                                     gtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGt                                                    >  1:915138/1‑105 (MQ=255)                                                                                                                     gtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGt                                                    <  2:915138/105‑1 (MQ=255)                                                                                                                     gtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGt                  <  2:194329/139‑1 (MQ=255)                                                                                                                     gtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGt                  <  2:552492/139‑1 (MQ=255)                                                                                                                     gtgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGt                  <  1:1585807/139‑1 (MQ=255)                                                                                                                      tgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTg                                                                                                          <  2:412102/50‑1 (MQ=255)                                                                                                                      tgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTg                                                                                                          >  1:412102/1‑50 (MQ=255)                                                                                                                      tgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCg                                                                                             >  1:1052104/1‑63 (MQ=255)                                                                                                                      tgCGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCg                                                                                             <  2:1052104/63‑1 (MQ=255)                                                                                                                        cGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             >  1:1639718/1‑45 (MQ=255)                                                                                                                        cGTTTTGTGCATCACAGCATCCCGGATAAGAAAGATATTCCGTgg                                                                                                             <  2:1639718/45‑1 (MQ=255)                                                                                                                             tgtgCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCg                                                      >  2:849406/1‑95 (MQ=255)                                                                                                                             tgtgCATCACAGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCg                                                      <  1:849406/95‑1 (MQ=255)                                                                                                                                    cacaGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGt                                                                                               >  1:1501524/1‑47 (MQ=255)                                                                                                                                    cacaGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGt                                                                                               <  2:1501524/47‑1 (MQ=255)                                                                                                                                    cacaGCATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGT‑AGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGAt  >  2:1390942/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        CGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGCACCCGTTCTTCGGCATTGCGATTTTCGTCGCACTGATGTTTATGTTTGTGCGTTTTGTGCATCACAACATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGAT  >  NC_000913/1551438‑1551708

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