Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,379,710	C→A	P21P (CCG→CCT)	frdC ←	fumarate reductase membrane protein FrdC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,379,710	0	C	A	100.0%	355.0 / NA	100	P21P (CCG→CCT)	frdC	fumarate reductase membrane protein FrdC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (52/48);  total (52/48)

TTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAA  >  NC_000913/4379576‑4379844                                                                                                                                          |                                                                                                                                       ttttGTAAAAAGTCGACGAATCCCGCCCAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCaa                                                                                                                                   >  1:1871926/1‑140 (MQ=255)    tGTAAAAAGTCGACGAATCCCGCCCAGGCTTCC‑GGGCCATTTTT‑CAGGTCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAAttt                                                                                                                                >  1:1204382/1‑140 (MQ=255)     aaaaaaaGTGGACG‑ATCCCCCCCCGGGTTCCGGGGCTATTTTG‑AAGGAAAAC‑CACCCCAAAATTCAGTCAATTGCTGAACCACACAGGGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAAtttt                                                                                                                               <  2:1151016/138‑1 (MQ=255)      aaaaaaGTCGACGAACCCCCCCCGGGCTTCG‑GGGCCTTTTTT‑AGGGGCAAA‑CAGCCCAAAAATCATTTCAATGCGAAACCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAAttttt                                                                                                                              <  2:1354055/139‑1 (MQ=255)       aaaaaGTCGACGAATCCCGCCCAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAAtttttt                                                                                                                             >  1:892762/1‑140 (MQ=255)          agcgcgacaaACCCCGCCCAGCCTCCC‑GGGCCATTTTC‑CAGGGCAAA‑CAGCCGGAAAATAAGTTCAATGCTAACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTcca                                                                                                                          <  2:1449695/136‑1 (MQ=255)          aaGTCGACGAATCCCGCCCAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTcca                                                                                                                          >  2:953683/1‑140 (MQ=255)           aGTCGACGAACCCCCCCCAGGCTCCC‑GGGCTTTTTTC‑GAGGGCAAA‑CACCCAAAAAATCAGTTAAAGGCAGACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTccac                                                                                                                         <  2:2398723/140‑1 (MQ=255)                  aaacccccccCAGGCTTCG‑GGGCTTTTTTT‑CAGGGAAAA‑CAGCCCAAAAATCAGTTAAATGCTAAACCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTgg                                                                                                                  <  2:872801/136‑1 (MQ=255)                  aaaccccGCCCAGGTTTCG‑GGGCTTTTTTT‑CGGGGCAAA‑CAGCCAAAAAATCAGTTCAATGCTAACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTgg                                                                                                                  <  2:28103/136‑1 (MQ=255)                          cccAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCAtt                                                                                                          >  2:117894/1‑140 (MQ=255)                          cccAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCAtt                                                                                                          >  1:1911065/1‑140 (MQ=255)                          cccAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGACAATTTTTTCCACCAGGTGGACGTCAtt                                                                                                          >  1:1834104/1‑140 (MQ=255)                            cAGGCTTCC‑GGGCTTTTTTT‑CAGGGCAAA‑CAGCCAAAAAATGAGTTCAATGCTAAACCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTgg                                                                                                        <  2:1345956/140‑1 (MQ=255)                             gggCTTCC‑GGGCCTTTTTT‑CAGGGAAAA‑CACCCAAAAAATCATTTCAATGCTAAACCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGc                                                                                                       <  2:1440878/139‑1 (MQ=255)                              gggttCC‑GGGCCTTTTTTCCAGGACAAC‑CACCCCAAAAATAATTAAAATGCAAACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGc                                                                                                       <  2:600820/137‑1 (MQ=255)                              ggCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGcc                                                                                                      >  1:2052535/1‑140 (MQ=255)                               gCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAAATTTTTCCACCAGGTGGACGTCATTGGCCg                                                                                                     >  2:2387145/1‑140 (MQ=255)                                   cc‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACa                                                                                                 >  1:374299/1‑140 (MQ=255)                                     ggggCCTTTTTT‑CGGGACAAA‑CACCCAAAAAATCAGTTAAATGCTAACCCACACAGCGGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACAt                                                                                                <  2:1004132/139‑1 (MQ=255)                                      gggcttttttc‑CAGGGCAAA‑CGGCCAAAAAATCAGTTAAATGCTAACCCAAACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATa                                                                                               <  2:1746112/134‑1 (MQ=255)                                      gggCCTTTTTT‑CGGGGCAAA‑CAGCCAAAAAATCATTTAAATGCAAACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATa                                                                                               <  2:1911065/140‑1 (MQ=255)                                      gggCCTTTTTC‑CGGGGCAAA‑CAGCCGAAAAAGAAGTTAAATGCTAACCCCCGCGGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGACAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATa                                                                                               <  2:1834104/140‑1 (MQ=255)                                          ttttttt‑CGGGCAAAA‑CACCCAAAAAATCATTTAAATGCTAAACCCCACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGt                                                                                           <  2:1152996/138‑1 (MQ=255)                                          ctttttt‑CAGGGAAAA‑CAGCCAAAAAATCATTTAAATGCTAAACCCCACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGt                                                                                           <  1:1998685/138‑1 (MQ=255)                                           aTTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGtt                                                                                          >  1:1552800/1‑140 (MQ=255)                                            ttttt‑CAGGGCTAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGttt                                                                                         >  1:662997/1‑140 (MQ=255)                                            ttttt‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGttt                                                                                         >  1:1020747/1‑140 (MQ=255)                                            ttttt‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGttt                                                                                         >  2:24206/1‑140 (MQ=255)                                                   aGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGttt                                                                                   >  1:1191260/1‑140 (MQ=255)                                                   aGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGtct                                                                                   >  1:1004132/1‑138 (MQ=255)                                                     ggCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAg                                                                                 >  1:611241/1‑140 (MQ=255)                                                      gCAAACCAGCCGAAAAACCAGTTAAATGGTGACCCACCCGGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAg                                                                                 <  2:1191260/140‑1 (MQ=255)                                                        aaa‑CAGCCCAAAAATCAGTTCAATGCTAACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcg                                                                              <  2:174774/140‑1 (MQ=255)                                                        aaa‑CAGCCAAAAAATCATTTAAATGCTAACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcg                                                                              <  2:662997/140‑1 (MQ=255)                                                        aaa‑CAGCCAAAAAAAAATTTCATTGCAAAACCACACCGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcg                                                                              <  2:1096359/140‑1 (MQ=255)                                                        aaa‑CACCCCAAAACCCATTCCAATGCGAACCCACACACCCGAAGCCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcg                                                                              <  2:564213/140‑1 (MQ=255)                                                          a‑CAGCCCGAAAATCATTTAAATGCTGACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcgtc                                                                            <  1:1503440/140‑1 (MQ=255)                                                          a‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcgtc                                                                            >  1:174774/1‑140 (MQ=255)                                                               cccGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGt                                                                        >  1:1096359/1‑140 (MQ=255)                                                                ccGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGtt                                                                       >  1:371141/1‑140 (MQ=255)                                                                 cGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTg                                                                      >  1:308720/1‑140 (MQ=255)                                                                  aaaaaTCAGTTCAATGCAGACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGc                                                                     <  2:2052535/139‑1 (MQ=255)                                                                   aaatTCAGTCAAATGGTGACCCCCACAGCGGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCa                                                                    <  1:1360306/140‑1 (MQ=255)                                                                   aaaaTCATTTAAATGCTGACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCa                                                                    <  1:959171/140‑1 (MQ=255)                                                                   aaaaTCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCa                                                                    <  1:24206/140‑1 (MQ=255)                                                                             aaaTGCTAACCCCCACAGCAGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCg                                                          <  2:2238065/139‑1 (MQ=255)                                                                             aaaTGCTAACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCg                                                          <  2:892762/139‑1 (MQ=255)                                                                                 gCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGt                                                      >  1:2194510/1‑140 (MQ=255)                                                                                     aaCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGc                                                           >  1:2395285/1‑131 (MQ=255)                                                                                     aaCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGc                                                           <  2:2395285/131‑1 (MQ=255)                                                                                     aaCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGttt                                                    >  1:882371/1‑138 (MQ=255)                                                                                     aaCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGttt                                                    <  2:882371/138‑1 (MQ=255)                                                                                     aaCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCa                                                  >  1:1021878/1‑140 (MQ=255)                                                                                     aaCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAAATTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCa                                                  <  1:2387145/140‑1 (MQ=255)                                                                                        cacacaGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAggg                                               <  2:1728212/140‑1 (MQ=255)                                                                                         acacaGCGGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGt                                              <  2:2383466/140‑1 (MQ=255)                                                                                         acacaGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGt                                              <  2:619579/140‑1 (MQ=255)                                                                                          cacaGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTc                                             >  1:1725831/1‑140 (MQ=255)                                                                                          cacaGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTc                                             >  1:1951867/1‑140 (MQ=255)                                                                                             aGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGt                                              >  1:1650660/1‑136 (MQ=255)                                                                                             aGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGt                                              <  2:1650660/136‑1 (MQ=255)                                                                                               ccGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTCTCAGGGTCGCGAt                                        >  2:1049113/1‑140 (MQ=255)                                                                                                    aCCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGaa                                   <  2:2194510/140‑1 (MQ=255)                                                                                                         cgTGCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCt                              >  2:471094/1‑140 (MQ=255)                                                                                                            gCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTc                           <  1:267856/140‑1 (MQ=255)                                                                                                            gCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTc                           <  2:1552800/140‑1 (MQ=255)                                                                                                            gCCTTCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTc                           <  1:1340851/140‑1 (MQ=255)                                                                                                               ttCGCGCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGaa                        >  1:969130/1‑140 (MQ=255)                                                                                                                 cgcgCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCg                          >  1:192734/1‑136 (MQ=255)                                                                                                                 cgcgCAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCg                          <  2:192734/136‑1 (MQ=255)                                                                                                                   cgcAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACttt                    <  2:2393582/140‑1 (MQ=255)                                                                                                                   cgcAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACttt                    <  1:117894/140‑1 (MQ=255)                                                                                                                   cgcAGCATGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTCTCAGGGTCGCGATAAGAAAGTCTTTCGAACttt                    <  1:1049113/140‑1 (MQ=255)                                                                                                                      agcaTGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTa                 >  2:1894645/1‑140 (MQ=255)                                                                                                                        caTGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGa                         <  1:1295030/130‑1 (MQ=255)                                                                                                                        caTGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGa                         >  2:1295030/1‑130 (MQ=255)                                                                                                                        caTGTAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACt               <  2:1082179/140‑1 (MQ=255)                                                                                                                            tAAAAGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGc           >  2:109342/1‑140 (MQ=255)                                                                                                                             aaaaGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGcc          <  2:1020747/140‑1 (MQ=255)                                                                                                                                aGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGccc         <  1:1448687/138‑1 (MQ=255)                                                                                                                                aGCGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGccc         >  2:1448687/1‑138 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACttt                    >  1:2314927/1‑125 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACttt                    <  2:2314927/125‑1 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTg            <  2:340006/133‑1 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTg            >  1:340006/1‑133 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct     >  2:1551026/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct     >  2:437715/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct     >  2:461386/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct     >  1:1356537/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct     >  2:1564500/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct     >  2:157484/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct     >  2:794372/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTAATTTGCCctgct     >  2:2116586/1‑140 (MQ=255)                                                                                                                                  cGATAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTAATTTGCCctgct     >  2:2116587/1‑140 (MQ=255)                                                                                                                                    aTAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGa   >  2:845813/1‑140 (MQ=255)                                                                                                                                    aTAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGa   >  1:2144407/1‑140 (MQ=255)                                                                                                                                    aTAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGa   >  2:1561938/1‑140 (MQ=255)                                                                                                                                    aTAAAAAGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGa   >  2:2463650/1‑140 (MQ=255)                                                                                                                                     tAAAAAGGCATTTTTTTCCACCAGGGGGAGGTCTTGGGCGGAAAATAGGGTTTAGGTTTAGCCGCCATGTTGCACTCCTTAGCGTGGTTTCAGGGCGGCGTTAGGAAAGTCTTTCGAACTTTTTACTTTGCCCTGCTGaa  >  2:710459/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       TTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCC‑GGGCCATTTTT‑CAGGGCAAA‑CAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAA  >  NC_000913/4379576‑4379844

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