Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,367,368	G→T	17.7%	D100Y (GAT→TAT)	arnC →	undecaprenyl phosphate‑L‑Ara4FN transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,367,368	0	G	T	17.7%	137.1 / 3.2	58	D100Y (GAT→TAT)	arnC	undecaprenyl phosphate‑L‑Ara4FN transferase
Reads supporting (aligned to +/- strand):  ref base G (10/37);  new base T (9/2);  total (19/39)
Fisher's exact test for biased strand distribution p-value = 3.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCTTACCGGAATTAATCAGGCGCACCACCACAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGA  >  NC_000913/2367133‑2367606                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    aaGCTTACCGGAATAAACCGGGCGCACACCCACAGCCTGGGAATCGTTGGGGAAAGAGTCTGAGAGCCTGCTGATTGATGACGCCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATg                                                                                                                                                                                                                                                  <  2:1158672/239‑1 (MQ=255)         cGGAATTAATCAGGCGCACCACCACAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGGTTGTGACACGTTACTGGAGACTTTATTATTACACTT‑GATGCCggtgt                                                                                                                                                                                                                                           >  2:845536/1‑234 (MQ=255)                   cAGGCGCACCACCACAGCCTGTGAATCGTTGGGGAAAGAGGAGGAGACCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGc                                                                                                                                                                                                                                <  2:533077/239‑1 (MQ=255)                       cgcACCGCCACAGCCTGTGAATCGATGGGGAAAGAGTAGGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGcc                                                                                                                                                                                                                               <  2:672769/236‑1 (MQ=255)                                       gtgAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCtg                                                                                                                                                                                                          <  1:1010417/239‑1 (MQ=255)                                         gAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTTAGTCACCGTTTTTGCGGCTTAATTATTTTCCTT‑TATGCCCCTCTCCAGAATTCCCCCG‑AAAAAATC‑CCCCCCCcggt                                                                                                                                                                                                        >  1:444266/1‑239 (MQ=255)                                             cGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCa                                                                                                                                                                                                    <  1:75532/239‑1 (MQ=255)                                                      aaGAGTATGAGTCCCTGCTGTTTGTTGCCGCAAG‑TAGCGATAATCCC‑GCCCATAGGCTGGTCGAAGCCTCACAGGGGGAGACC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGCCACGTTACTGGCGACTTAATTATTCCCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGAt                                                                                                                                                                                           <  2:47417/239‑1 (MQ=255)                                                         ggTATGAGAGCCGGCTTATTGACGCCGGCAGTTAGC‑AAAATTCCGGCGCATATGTTGTCCAAACCCCCAAAGGGGAAAAAC‑AGCCATTTTGTTTCTATTTTGTTAAACCGAAATTCCGGGCAACATTCAGGGATT‑T‑GGGGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATg                                                                                                                                                                                          <  2:11831/236‑1 (MQ=255)                                                         agTATGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTTATTCACGTTACTGGGGACCTTATTATTATCCTT‑TGTGCCGCTCCTCCCAATCCGCCCA‑AAGAAAAC‑CCCCCCCCGGTGGCTAAAAcccctgga                                                                                                                                                                                        >  2:1139454/1‑233 (MQ=255)                                                                tacccgcttgattaatGCCCGCAG‑GTGAGATATTCCC‑GCGAATATCCTGTGCTAAGCCTCTAAAGCGGAGACC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGAAACATTCAGCGATTAT‑GGGGGGATTCAGTAACGTTACTGGCGACTTAATTATTACCCTG‑GATGCCGAGCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacg                                                                                                                                                                                 <  2:43185/231‑1 (MQ=255)                                                                  tcttgctgATTGAGACCGGCAG‑TGGCGATAATCCC‑GCGCATATGCTGGTCAAAGCCCCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                               <  2:439051/236‑1 (MQ=255)                                                                  tCCTGCTGATTGAGGACGGCAG‑TAGCAAAAATCCC‑GCGAATATCCTGGTGGAAGCCCCACAAGCGGAGACC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                               <  2:306655/239‑1 (MQ=255)                                                                  tCCTGCGTATGGAGGACGGCGG‑TGGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACGAGCGGAGAAC‑AGCCATATTGTGTTTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTAACGTTACTGGCGACTTAATTATTGCCCTT‑GATGCCGATCTCCAGAAACCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                               <  2:589613/239‑1 (MQ=255)                                                                  gccTGCTGATTGATGACGGCAG‑TAGCAATAATTCC‑GCGCATATGCTGGTCAAAGCCCCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                               <  2:778735/238‑1 (MQ=255)                                                                   ccTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGCTTTAGTCACGTTATTGGCGGCTTAATTATTAACACT‑TATGCCGCTCTCCAGAAACAGCCCA‑AAGAAAAA‑CCCCCCCCGTGGGGAAAAACCCACGAAAGTTacgccg                                                                                                                                                                              >  1:1140137/1‑239 (MQ=255)                                                                     tgctgATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGTGTTTAATCACGTTACTTGCGACCTAATTTTTTCCCTT‑TTTGCCGGTCTCCCGAAACCACCCGAAAAAAAAC‑CCCCCCCTGGGGGCAAAA‑CCGATGAAAGGTacgacggg                                                                                                                                                                            >  2:918768/1‑237 (MQ=255)                                                                     tgctgATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGGTACTGGCGACTTCATTATTACCCCT‑GATTGCCACCTCCAGGATCCCCCAC‑AAGAAATC‑CCCCCCCTGGTGGCGAAAGCCGACGAAGGGTacaacggc                                                                                                                                                                            >  1:255111/1‑237 (MQ=255)                                                                        tgattgatGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGt                                                                                                                                                                          <  1:82400/238‑1 (MQ=255)                                                                          attgatGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTTAGTCACGTTACTGGCGACTTAATTTTTAACCTT‑TTTGCCGATCTCCCGAATCCGCCAG‑AAGAAAAC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGGTacccactgggagg                                                                                                                                                                       >  1:59281/1‑228 (MQ=255)                                                                           ttgatGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTCCCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAgg                                                                                                                                                                       <  1:94127/238‑1 (MQ=255)                                                                            tgatgaCGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTTAGTCACGTTTCTGGCGGCTTTATTTTTAACCTT‑TATGGCGATCTTCAGAAAACGCCCG‑AAGAAATC‑CACCCCCTGGTGGCGAAAAACGATGAAGGTTACGACGGGGTAgggg                                                                                                                                                                     >  1:399787/1‑238 (MQ=255)                                                                            tgatgaCGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGa                                                                                                                                                                     <  1:512218/239‑1 (MQ=255)                                                                            tgaggaCGGCAG‑TGGCGATAATTCC‑GCGCATATGCTGGTCGAGGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGa                                                                                                                                                                     <  2:618158/239‑1 (MQ=255)                                                                                   ggCAG‑TAGCGATAATCCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACg                                                                                                                                                              <  1:641864/239‑1 (MQ=255)                                                                                                                                aaGGGGAGAACAAGCATTTTTGTGTCTATTTGGCTTAACCCAAATTACGTGCAAGATTCAGCGATTAT‑GGCGGGTTTCAGCCACGTTGCTGGCAACTTAATTTTTCCCCTT‑GATGCCGTACTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGaa                                                                                                                  <  2:930752/239‑1 (MQ=255)                                                                                                                                        acc‑AGCAATTTTGTTTCTATTTTGCTTAACCGCAATTCCGGGAAACATTCAGCGATTAG‑GGGGGGATTCAGCCACGTTACTGGCGATTTAATTATTCCCCTG‑GATGCCAATCTCCAGAATCCGCCAG‑AAGAAACC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTaa                                                                                                         <  2:59281/237‑1 (MQ=255)                                                                                                                                         accatctattttGTGTCTATTTTGCTTAACCGAAATTCCGGGCAACATTCAGCGATTAT‑GGGGGGATTCAGTCACGTTACTGGCGACTTAATTATTCCCCTT‑GATGCCGATCCCCAGAATCCGCCAG‑AAGAAACC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTaa                                                                                                         <  1:918768/229‑1 (MQ=255)                                                                                                                                              cctttttGGGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGCCCCGTTACTGGCAACTTAATAATTCCCCTT‑GATGCCAATCCCAAGAACCCGCCAA‑AAGAAA‑C‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGc                                                                                                    <  2:508671/233‑1 (MQ=255)                                                                                                                                                 ttttGTGTGTTTTTTGCTTAACCCAAATACCGGGCACCATTCAGCGATTAG‑GGCGGAATTCAGTCACGTTCCTGGCGACTTAATTATCCCCCTG‑GAGGCCGATCCCCAGAATCCCCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTGGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTa                                                                                                 <  2:399787/237‑1 (MQ=255)                                                                                                                                                   ttGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGGTTTAGTCACGGTATAAGGGAATCAATTATTTATCTT‑TAAGCCGAGCTCCACAAACCCCCAG‑AAGGAATCACCCCCCCCGGGGGTAAAA‑ACGAAGAAGGGTAAGCAGGGGGATGGACTGCACGCCAG‑AAACGCCAGGAAAGCGGGTTTCGTAACAACGCATTTGAAgaggataatcggcgttc                                                                                                 >  2:185516/1‑227 (MQ=255)                                                                                                                                                        tCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGGTTCAGTCACGTTACTTGCGGCTTTATTATTACCCTT‑TATGCCGATCTCCAGAATCCCCCAG‑AAGAAATC‑CCCCCCCTGGTGGCAAAAACCCATGAAAGTTACGACGTGGGAGGGACTGTACGCCCG‑AACCGCCCGGACAGCCGGTTTCGTAAAACCCC‑TTCGAAGATGATTAACCGGCTTATTCatc                                                                                           >  1:589919/1‑236 (MQ=255)                                                                                                                                                        tCTATTTTCCTAACCCCCAATCGC‑GGCAACAATTAGCGGTTTTGGGAGGAATTCAGTCCCGTTTGGGTCGCCTTACTTCTTACCCTG‑TATCCCGTCCCCAAAAACCCCCCAG‑AAGAAACC‑CCCCCCCTGTTGGCAAAAGCAGATGAAGTTAACGACGTGTTAGGACCTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑gcg                                                                                         <  2:298486/239‑1 (MQ=255)                                                                                                                                                        gtgatttgacTTACCCGCAATAACGGGCAAGATACAGCATTTAG‑GGCGGGATTCAGTCACGTTACGGACGACTAATTGTTTCCCCTT‑GAGGCCAAGCCCCAGATTCCGCAAG‑AAGAAACC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTGCGACGTGGAAGGGACTGTACGCCAG‑ATCCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑gcg                                                                                         <  2:255111/230‑1 (MQ=255)                                                                                                                                                                  ttAAACGAAAATACGGGCAACATACATCGATTAT‑GCGGGGATTCGTTCCCTTTGCTGGTGCTTTAATTTTTCCCGTG‑GAGTCCGATCTCCAGAACCCGCGAA‑AAAAAAGC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGTTTACGATGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGc                                                                               <  2:444266/239‑1 (MQ=255)                                                                                                                                                                    aaCCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGTGTTTAGTCCCGTTACTTGCGACCTAATTATTATCCTT‑TTTGCCGATCTCCCGAATCCGCCAG‑CAGAAAAC‑CCCCCCCTGGTGGCCCAAGCCGGCGAAAGTTACGACGTGGGTGGGACGGTTCGGCAGAAACCCCCAGG‑CCGCCGGGTTCCTTAAACCCC‑TCCCGAGAAGAGTAAACCGCTTATTCATGCGcacaacggcac                                                                              >  2:333737/1‑234 (MQ=255)                                                                                                                                                                        gCAATTCCGGGCACATTTCAGGGATTAT‑GGCGGGATCCAGCAACGTTCCTGGCAACTAAATTATTCCCCTT‑GAGGCCGATCCCCAAAACCCGCCAG‑AAAAATCC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTAGGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCg                                                                         <  2:522021/239‑1 (MQ=255)                                                                                                                                                                            ttACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATgg                                                                     <  1:583008/239‑1 (MQ=255)                                                                                                                                                                                  caacaagtCAGCGTTTAG‑GGGGGGATCCAGTCAGTTTACTGGCGACTTAATTATTCCCCTG‑GAGGCCGATCCCCAGAACCCCCCAG‑AAAAAACC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGAt                                                                <  2:538376/231‑1 (MQ=255)                                                                                                                                                                                   aaaCAGTCAGGGATTAG‑GGGGGAAGCCAGCCGCGTTGTGGGCGACATAATTATTCCCCTG‑GTTGCCGAGCTCAGGAACCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGAAAAAGCCGAGGAAGGTTACGACGTGGTAGGGCCGGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTa                                                              <  2:987187/238‑1 (MQ=255)                                                                                                                                                                                     aaaTTCAGCGTTTAT‑GGGGGGTTTCAGCAACGTTCCTGGGGATTAATTTTTTCCCCTG‑GAGGCCGATCCCAAGAATCCCCAGG‑AAGAATCC‑CCCCGCCGGGGGGCAAAACCCGATGAGGGTTACGAGGTGGTAGGGACTGTCCGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACg                                                            <  2:589919/237‑1 (MQ=255)                                                                                                                                                                                              ggttAG‑GGCGGAATTCAGCCGTGTTCGT‑GCGGCTTAAATTATACCCTTAGAGGCCGATCTCCAAAACCCTCCAA‑AAGAAATC‑CCCCGCCGGGTGGAAAAAGCCGATGAAGGTTACGACGTGGTGGGGCCTGTACGCCAG‑ACCCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTa                                                      <  2:1206403/234‑1 (MQ=255)                                                                                                                                                                                                 tAT‑GGCGGGATTCAGTCACGTTCCTGGCGACTTAATTTTTCCCCTT‑GATCCCGATCTCCAGAAACCGCCAG‑AAGAAACC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTGGGGATTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTAtgctgc                                                <  1:928314/239‑1 (MQ=255)                                                                                                                                                                                                           aTTCAGTCACGTGACGGGAGATTTAATTATGCCCCTG‑TGTGCAGATGTCCAGAACCCGCCAG‑AAGAAATC‑CCCCGCCTGGGGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTAt                                        <  2:1070730/238‑1 (MQ=255)                                                                                                                                                                                                                    cGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAACCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCata                              <  1:1139454/239‑1 (MQ=255)                                                                                                                                                                                                                     gTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTa                                                              <  1:720222/206‑1 (MQ=255)                                                                                                                                                                                                                     gTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTa                                                              >  2:720222/1‑206 (MQ=255)                                                                                                                                                                                                                         cTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTc                         >  1:297502/1‑239 (MQ=255)                                                                                                                                                                                                                          tGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcg                        >  2:460712/1‑239 (MQ=255)                                                                                                                                                                                                                           ggCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AGGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGGAGGGACGGTACGCCAG‑AACCGCCAGGACAGGTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCTCACCACTGGCAAAGCGATGGGTGATTACGGTGGTATTCTTCGCGCCTCTCGCCGTCATATTGTcgt                       >  2:236727/1‑238 (MQ=255)                                                                                                                                                                                                                              gACTTAATTATTCCCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatg                     <  1:845536/238‑1 (MQ=255)                                                                                                                                                                                                                               aCTTATTTATTCCCCTG‑GATGCCGATCCCCAAAAGCGGCAGA‑AAGAAACC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGTTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcg                   <  2:297560/239‑1 (MQ=255)                                                                                                                                                                                                                               aCTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCCTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcg                   >  1:872031/1‑239 (MQ=255)                                                                                                                                                                                                                                 ttAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcgat                 >  1:475296/1‑239 (MQ=255)                                                                                                                                                                                                                                       attcCCCTT‑GATGCCGATCTCCAAAATCCCCGAG‑AAGAAATC‑CCCCGCCTGGTGCCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACGGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCa           <  1:185516/239‑1 (MQ=255)                                                                                                                                                                                                                                         tccccGT‑GATGCCAATCGCCGAAACCCGCCAA‑AAAAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt         <  1:518576/237‑1 (MQ=255)                                                                                                                                                                                                                                          aCCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt         >  1:158310/1‑238 (MQ=255)                                                                                                                                                                                                                                             cTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCa     >  1:896004/1‑239 (MQ=255)                                                                                                                                                                                                                                                 gATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa  >  2:1038544/1‑239 (MQ=255)                                                                                                                                                                                                                                                 gATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa  >  1:651384/1‑239 (MQ=255)                                                                                                                                                                                                                                                 gATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa  >  1:651234/1‑239 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    AAGCTTACCGGAATTAATCAGGCGCACCACCACAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAG‑TAGCGATAATTCC‑GCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAAC‑AGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTAT‑GGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTT‑GATGCCGATCTCCAGAATCCGCCAG‑AAGAAATC‑CCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAG‑AACCGCCAGGACAGCTGGTTTCGTAAAACCGC‑TTCGAAGATGATTAACCGGCTTATTCA‑GCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGA  >  NC_000913/2367133‑2367606

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