Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,329,682	G→A	V117V (GTG→GTA)	sohB →	S49 peptidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,329,682	0	G	A	83.3%	74.2 / 7.5	36	V117V (GTG→GTA)	sohB	S49 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base A (13/17);  total (14/22)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAACGCCAGCGTGGCGAGTTACGGGTCAACAATCTCAGCGAACAGTATAAGGAGATGAAAGAAGAACTGGCCGCGGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCGTCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTT  >  NC_000913/1329436‑1329905                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                taaACGCCAGCGTGGCGAGTTACGGGTCAACAATCTCAGCGAACAGTATAAGGAGATGAAAGAAGAACTGGCCGCGGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACt                                                                                                                                                                                                                             >  1:59817/1‑251 (MQ=255)                 gAGTTACGGGTCAACAATCTCAGCGAACAGGATAAGGAGATGAAAGAAGAACTGGCCGCGGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCAGGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGAt                                                                                                                                                                                                             <  2:50170/251‑1 (MQ=255)                      aCGGGTCAACAATCTCAGCGAACAGTATAAGGAGATGAAAGAAGAACTGGCCGCGGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCGTCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACg                                                                                                                                                                                                         <  1:353844/250‑1 (MQ=255)                                                  aaGGAGATGAAAGAAGAACTGGCCGCGGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGa                                                                                                                                                                            <  1:113613/251‑1 (MQ=255)                                                     gagaTGAAAGAAGAACTGGCCGCTGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCa                                                                                                                                                                         <  1:605567/251‑1 (MQ=255)                                                                      ggCCGCGGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACGCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTgt                                                                                                                                                        <  2:337452/251‑2 (MQ=255)                                                                            ggcgCTGATGGACTCACCTCAGCAAAAACCGTGGCACAAAGCGCAGAAGAAACAGCATAAGCAAGAAGCCAACGCCGAACACGCGAAAGCCAAACTGGGAGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTCCGTGAAGAGATACCGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGaaatcc                                                                                                                                                  >  2:611823/1‑248 (MQ=255)                                                                            ggcgCTGATGGACTCACCTCAGCAAAAACAGTGGCACACAGCGCAGAAGAAAAAGCATCAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAACCTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGAAGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCACACCGCAGGACAAGGTTGTGATCCGTCTGGAAAGcc                                                                                                                                                  >  2:276166/1‑251 (MQ=255)                                                                               gcTGATGGCCTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGGGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCtg                                                                                                                                               >  1:327742/1‑251 (MQ=255)                                                                                   aTGGACTCACATCAGCAAAAAAAGTGGCACAACGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCtggtgg                                                                                                                                           >  1:513793/1‑251 (MQ=255)                                                                                       aCTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatg                                                                                                                                       <  1:160562/251‑1 (MQ=255)                                                                                                            ggCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATTGTGCATGGTTACGGGCTggc                                                                                                                   <  2:385381/250‑1 (MQ=255)                                                                                                            ggCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggcg                                                                                                                  >  1:226496/1‑251 (MQ=255)                                                                                                            ggCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggc                                                                                                                   <  1:100188/250‑1 (MQ=255)                                                                                                               acaAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACGGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCGTCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCt                                                                                                               <  2:594119/251‑1 (MQ=255)                                                                                                                   aGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGc                                                                                                           <  1:583676/251‑1 (MQ=255)                                                                                                                         gaagaaAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGc                                                                                                     <  2:327742/251‑1 (MQ=255)                                                                                                                             aaaaaGCATAAGCAAGAAGCGACGGCAGCAAAAGCGAAAGCCAAACTGGGAGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGACCCGTCTGGAAAGCACTGGTGGCATGGTGCATGGTTACGGGCGGGCGGCTTCGCAACTGCAGCg                                                                                                 >  2:496107/1‑251 (MQ=255)                                                                                                                                     tAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGcgcg                                                                                         >  1:152733/1‑251 (MQ=255)                                                                                                                                       aGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAATCCCTGGTGGCATGGTGCATGTTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGAt                                                                                       <  2:441468/251‑1 (MQ=255)                                                                                                                                                     gcagcaAAAGCTAAAGCCAAACTGGGCGATGTGGCAACGGACAGTAAACCCC‑GCATCTGGGTGTTGGATTTTAAAGGCAGCATGGACGCCCATGAAGGAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCGGGTGGCATGGTGCATGGTTACGGGCTGGCGGCGTCGCAACTGCAGCGGCTGCGCGATAAAAACATTCCttt                                                                         <  2:81914/251‑1 (MQ=255)                                                                                                                                                                gAAAGCCAAACTGGGCGAGGTGGCAACTGCCAGTAAACCCC‑GCGTCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTAAGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCAGGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGAt                                                                                       <  2:133209/226‑1 (MQ=255)                                                                                                                                                                gAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCGTCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGTCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCAAGTTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGAt                                                                                       >  1:133209/1‑226 (MQ=255)                                                                                                                                                                               ggAGGTGGCAAATGACAGTTAACCCC‑GCGTGTGGGTGCTGGATGGTAAAGGCAGCATGGAAGCCCATGAAGTGAACTCGCTGCGGGAAGAGATAACGGCTGTAGTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCGGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATGCCTTTAACTGTTACGGTAGACAAAGTCGCTg                                               <  2:26309/250‑1 (MQ=255)                                                                                                                                                                               cGAGGTGGCAACTGACAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTg                                               >  1:536500/1‑251 (MQ=255)                                                                                                                                                                                 aGGTGGCAACTGACAGTCCACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGcc                                             >  1:170960/1‑251 (MQ=255)                                                                                                                                                                                  ggtagCAACTGACA‑TAAACCCCCGCGTCTGGGTACTTGATTTGAAATGCAGCATGGGCGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACGCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTTCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGGAGACAAAGTCGCTGCCa                                            <  2:186262/247‑1 (MQ=255)                                                                                                                                                                                             aCAGTAAACCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCGTTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCCGCGGCGGCTAc                                 >  2:492768/1‑251 (MQ=255)                                                                                                                                                                                                  aaaCCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGcggcgggatacaagt                            >  1:612459/1‑242 (MQ=255)                                                                                                                                                                                                  aaaCCCC‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGGGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATGCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACatga                             <  2:226496/250‑1 (MQ=255)                                                                                                                                                                                                      ccc‑GCATCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGcc                        <  2:318085/251‑1 (MQ=255)                                                                                                                                                                                                           catcggggTGCTGGATTTTAAAGGCAGCATGGACGCCCAAGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCCGGTTGGGCTCCGTCGGGAAAGCCCTGGTGGCATGTTGCATGGTTACGGGCTGGCGGCTGCGCAACTGCATCGTCTGCGCGATAAAAGCATTGCTTTAACTGGTACGGGAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCtgtg                    <  2:271015/246‑1 (MQ=255)                                                                                                                                                                                                                gggTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTgg                 <  1:451785/249‑1 (MQ=255)                                                                                                                                                                                                                  gTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCGGTGTGGCGGa             <  2:14384/251‑1 (MQ=255)                                                                                                                                                                                                                      tGGATTTTAAAGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTTCCAGCGGCGGTTACaaaaaagcacatgttgggggaaa          >  1:54764/1‑228 (MQ=255)                                                                                                                                                                                                                              aaaGGCAGCATGGACGCCCATGAAGTAAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGAGAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGttt  <  1:276166/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                TAAACGCCAGCGTGGCGAGTTACGGGTCAACAATCTCAGCGAACAGTATAAGGAGATGAAAGAAGAACTGGCCGCGGCGCTGATGGACTCACATCAGCAAAAACAGTGGCACAAAGCGCAGAAGAAAAAGCATAAGCAAGAAGCGAAAGCAGCAAAAGCGAAAGCCAAACTGGGCGAGGTGGCAACTGACAGTAAACCCC‑GCGTCTGGGTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTT  >  NC_000913/1329436‑1329905

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