Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	15,226	C→T	G353G (GGC→GGT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,226	0	C	T	95.1%	166.9 / ‑2.0	65	G353G (GGC→GGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (23/35);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

ATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAG  >  NC_000913/14984‑15472                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         aTGGCGGCGCTGGGGGGCGAAATCGAAGTAACGAACCTTGATGGTCGCGTCAAACTGAAAGTGCCTGACGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTATGTCCGCGGGGGCGCACA‑GGGAGAATTGCTGTGCCGCGTTGTCGTCGAAACACCGGGAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGAAAGAAAGCTTCGGGGGCCCAAACGGTGAGccca                                                                                                                                                                                                                                                   >  2:444427/1‑247 (MQ=255)   ggcgtgggtgtgTGGCGCAACCGAAGTAACGACCCTTGATGTGCGCGTCAAGCTGATCGTGCCGGGCGAAACCCAGCCCGGTAAGCTGTTCCGTATGCGCGGTAAAGTC‑GGCAAGTCTGTACGCGGTGGCGCCCT‑GGGTGATTTGCTTTGCCGCGTTGTCGGCGAAACACCGGTAGGCCTGAACGAGAATCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTTGCCCATCAGGTGAGCacaaca                                                                                                                                                                                                                                                <  2:43044/240‑1 (MQ=255)                ggCGAAATCGAAGGACCGCCCCTTGATGGTCGCTTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAGGGC‑GTCAAGTCTGTCCGCGGTGGTGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGGCGAAACACCTGTAGACCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGTTTAGCACAACAGCCCGCGCTCaa                                                                                                                                                                                                                                    <  2:386906/250‑1 (MQ=255)                        cGAATTACCGACCCCTGATGGTCGCGTCAAACTGAAAGCGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGTTAAAGGC‑GTCAAGTCTGGCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑Gcttc                                                                                                                                                                                                                           <  2:573966/250‑1 (MQ=255)                        cGAAGTACCGACCCCTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGGCCGGTAAGCTTTGCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTGGCTGTGCCGCGGTGGCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑Gcttct                                                                                                                                                                                                                          <  2:532365/251‑1 (MQ=255)                           aGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTg                                                                                                                                                                                                                       <  1:199035/251‑1 (MQ=255)                           aGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTg                                                                                                                                                                                                                       <  1:238270/251‑1 (MQ=255)                               ccGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCCAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCCTTGATgg                                                                                                                                                                                                                   >  2:392563/1‑251 (MQ=255)                                cGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAACGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGC‑GCTTCTTTGCTggg                                                                                                                                                                                                                  >  2:53315/1‑250 (MQ=255)                                cGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGgt                                                                                                                                                                                                                  >  1:429318/1‑251 (MQ=255)                                 gACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCGGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGgtg                                                                                                                                                                                                                 >  1:140686/1‑251 (MQ=255)                                 gACCCTTGATGGACGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTATGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGGACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGAGAGCACAACAGCGCGCGCTCAAAGC‑GCTTCTTTGAGGgtg                                                                                                                                                                                                                 >  2:469343/1‑251 (MQ=255)                                  aCCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGgtgt                                                                                                                                                                                                                <  1:457498/251‑1 (MQ=255)                                       tGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTCGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTgaa                                                                                                                                                                                                             <  1:450551/249‑1 (MQ=255)                                           ggTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGtt                                                                                                                                                                                                       <  2:157971/251‑1 (MQ=255)                                             gcgcgTCAAACTGAAAGTGTCTGGTGAAATCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GGCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGATCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGttt                                                                                                                                                                                                      <  2:44961/249‑1 (MQ=255)                                              cgggTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGttttt                                                                                                                                                                                                    <  2:187374/248‑1 (MQ=255)                                              cgcgTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGttttt                                                                                                                                                                                                    >  1:276008/1‑251 (MQ=255)                                                cgTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAAAAGCCCGCGCTCAAAGA‑GCTTCGTTGATGGTGTGAAGAAGTTTTTTg                                                                                                                                                                                                  >  2:432761/1‑251 (MQ=255)                                                      accgAAAGTGCATGGCGAAACCCAAACCGGTAAGCTTTTCCGAATGCTCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACC‑GGGTGATTTGCTGTGCCGCGTTGACGCCGCAACACAGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGCGCTGCACGTAAGATTCGGTGGACAAACGGGTGAGAACTACAGACCGTGCACAACGA‑GCTTATTTGATGGGGATAAGGAGTTATTTGACGAcc                                                                                                                                                                                            >  2:517157/4‑251 (MQ=255)                                                        tGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCaaa                                                                                                                                                                                                                                   >  1:465769/1‑211 (MQ=255)                                                              tGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTACAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAACCAGCTGCTGAAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGGGGGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGCTGGGGTGGAGCAGGTTTTTGAAGACCTTACCCGc                                                                                                                                                                                    >  2:134636/1‑251 (MQ=255)                                                               gCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCt                                                                                                                                                                                   >  1:673448/1‑251 (MQ=255)                                                                   ggCGAAACCCAGACCGGTAAGCGTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTaa                                                                                                                                                                                 <  1:310865/249‑1 (MQ=255)                                                                    gCGAAACCCAGACCGTTAAGCTTTTCCGTATGGGCGGTAAAGGC‑GTCAAGTCGGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCt                                                                                                                                                                              <  2:190867/251‑1 (MQ=255)                                                                      gtaccccggaccggttggccttttcGTATGCGCGGTTCAGGCTGT‑AAGTCTTTCCGCTGTGGCGCCCA‑GGGTGATTTGCTGTGCAGAGTTGTCTTCGAAACACCGGTATGCCTGAACGAGCAGTAGAAACAGCTGTTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCAAACCGGTGAGCACAACCGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGATCCGCTAACCTcc                                                                                                                                                                            <  2:268430/227‑1 (MQ=255)                                                                       aaaCCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAGGA‑GCTTCTTTGATGGCGTGAAGAAGTTTTTTGACGCCCTGACCCGCTAACCTccc                                                                                                                                                                           >  2:146516/1‑251 (MQ=255)                                                                                     aatcttttacGTATGTGCGGTAAAGGC‑GGCAAGTCTGTCCGCGGTGTCGCACA‑GGGTGATTTGCTGTGCCGCGTGGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGGTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCg                                                                                                                                                             <  2:591745/242‑1 (MQ=255)                                                                                       tcatttCCGTATGCGCGGTAAAGGC‑GTCAAGTATGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCAGCGTTGTCGTCGAACCACCGGGAGGCCTGAAAGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGATGGCCCAACAGATGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGAGGGTGTGCAGAAGTTTTTTGACGACCTGCCCCGCTACCCTCCCCCAAAGCATGCCAGTg                                                                                                                                                           >  2:538321/5‑251 (MQ=255)                                                                                          tttCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTgg                                                                                                                                                          <  1:159333/248‑1 (MQ=255)                                                                                           ttCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGc                                                                                                                                                        >  1:400429/1‑250 (MQ=255)                                                                                             ccGTATGCGCGGTAAAGGT‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGTGTTGTCTTCGAAACACCGGTAGGCCTGAACGATAAGCAGAAACCGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGCAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCATCAGCCTTCCCGTTGGCAg                                                                                                                                                      <  2:400429/250‑1 (MQ=255)                                                                                              cGTATGCGAGGTAAAGGC‑GTAAAGTCTGACCGCGGAGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGAAGGCCTGAACGAGAAGAAGAAACAGCTGCTGAAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGAGAGCACAAAAGAACTCGCT‑AAACACGCTTCTTTGATGGTGTGTAGAAGTTTTTTGAAGCACTGACCCGCTCACCTATCCATAAGCCTGTCCggaggaaggc                                                                                                                                                    >  2:372615/1‑242 (MQ=255)                                                                                               gTATGCGCGGGAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGct                                                                                                                                                   <  1:163389/251‑2 (MQ=255)                                                                                                gaTGCGCGGTAAAGGC‑GGCAAGTCTGCCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTTTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACATCCCGCGCTCAAAGA‑GCTTCTGTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTACCCTCCCCAAAAGCCTGCCCGTGTGCAGGctt                                                                                                                                                  <  2:118754/250‑3 (MQ=255)                                                                                                     gcgGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGAGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGctggggt                                                                                                                                              >  1:481007/1‑244 (MQ=255)                                                                                                      gcaggataggggTCAAGTCTGGCCGCGGTGTCGC‑CACGGGGGCTGTGCTGGGCCGCGTTGGCGTCGATCAACCGGTAGGCCTGACTGCGAAGAAGAAACAGCTGCTGCAAGTTCTTCAAGAAAGCTTCTGTGGCACAACCGGTGTGCACAACAGCTAGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACTTTACCCGCTAACCTCCCCCATAGCCTGCCCGTGTGCAgttttgggta                                                                                                                                             <  2:532170/240‑10 (MQ=255)                                                                                                          aaaGGC‑GTCAAGGCAGTCCGCGGTGGCGCACA‑GGGTGATGTGCTGTGCCGAGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAATAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCGGCCCGTGGGCAGGCTTGGTTaaaaa                                                                                                                                         <  2:26773/250‑1 (MQ=255)                                                                                                               c‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTaaacttagggg                                                                                                                                   >  2:639051/1‑243 (MQ=255)                                                                                                                     aGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGtt                                                                                                                              <  1:89197/251‑1 (MQ=255)                                                                                                                     aGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGtt                                                                                                                              >  1:376911/1‑251 (MQ=255)                                                                                                                                              gTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTgaaa                                                                                                        >  1:342802/1‑246 (MQ=255)                                                                                                                                                 aTTTGCTGGGCCGTGTTGTCGTCGAAACTCCGGTCGGCCTGAACGAGAATCAGAAACAGCTGCTGCAAGAGCTGCAAGAAATCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTGTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCATGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTgg                                                                                                   <  2:429318/251‑1 (MQ=255)                                                                                                                                                    tGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCg                                                                                                 <  1:582369/250‑1 (MQ=255)                                                                                                                                                      cTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGCAAAGTGGCgggg                                                                                              >  1:516477/1‑251 (MQ=255)                                                                                                                                                           ccGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGagat                                                                                          >  1:59642/1‑247 (MQ=255)                                                                                                                                                              ggTTGCCGTCGAAACACCGGTAGGCCTGAATGAGAAGCAGAAACAGCTGTTGCAAGATCTGCAAGAAAGCTTCGGTGGACCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTGTGACGACCTGACTCGCTAACCTCCCCAAAAGCCTGCCCGTTGGCAGGATTGGGTAAAATTAGGGTTCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTcc                                                                                      <  1:455746/250‑1 (MQ=255)                                                                                                                                                                  gtcgtcGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                                                                   <  1:345227/250‑1 (MQ=255)                                                                                                                                                                    cgtcgAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAg                                                                                <  1:218810/251‑1 (MQ=255)                                                                                                                                                                        gAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGaacaaacccctcagacgc                                                                            >  1:229314/1‑234 (MQ=255)                                                                                                                                                                          aaccccGGTAGGCCTGTACAAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCATCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCGGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAATTGGCGGGGATCACTCTCATAAGCGCTaa                                                                          <  2:17864/247‑1 (MQ=255)                                                                                                                                                                          aaCACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAcaacccacaaggaca                                                                           >  2:574735/1‑236 (MQ=255)                                                                                                                                                                           acacCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAAc                                                                         <  1:451476/251‑1 (MQ=255)                                                                                                                                                                            cacCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                                        <  1:469343/251‑1 (MQ=255)                                                                                                                                                                                 gTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAgg                                                                  <  1:236775/251‑1 (MQ=255)                                                                                                                                                                                              agaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTAc                                                     <  1:538321/247‑1 (MQ=255)                                                                                                                                                                                                   cAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGACACTACCCCAAGCGCTAAATATAAGGGtgtgggatttgtcccg                                                 >  1:172517/1‑236 (MQ=255)                                                                                                                                                                                                     gAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTGACTCCCCTAAGCGCTAACT‑Taaagggtggggggttatgcctgtt                                              >  1:232317/1‑229 (MQ=255)                                                                                                                                                                                                              gctgcaagagcggcaagaAAGCTGCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGAGGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTa                                      <  2:376911/250‑1 (MQ=255)                                                                                                                                                                                                                       gctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                            <  1:53315/251‑1 (MQ=255)                                                                                                                                                                                                                           caagaAAGCTGCGTTGGCCCAACCGGTGAGCACACCAGCCCGCGATGTAAGA‑GCTTCTTTGATGGTGAGAAGTAGTTTTTTGACGGTCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTGGGGTAAAAGTAGGGTGCGTTGAAGATATGCGAGAACCTGAAAAGTGGCGGGGATCTCTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                        <  2:152633/251‑1 (MQ=37)                                                                                                                                                                                                                                     tcgctGGCCCAAACGGGGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGAACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAGGAATTACTCTCAc               <  2:129217/246‑1 (MQ=38)                                                                                                                                                                                                                                       ggTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGCTGGTGTGACGACGTTTTTTGACGACCTGACCCGCTAAACTCCCCAAAGGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTACTCCCAAAAACGCATACA‑TAAAGGGTGTGGCATTACGCCCCATCTGAATTTATGTGCCGATGAATTAATCTCACgct            >  1:286245/1‑249 (MQ=255)                                                                                                                                                                                                                                                aCCGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa   <  1:668826/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ccGGTGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg  <  1:432761/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         ATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGC‑GTCAAGTCTGTCCGCGGTGGCGCACA‑GGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGA‑GCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAG  >  NC_000913/14984‑15472

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