Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,401	C→T	G78G (GGC→GGT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,401	0	C	T	91.3%	108.9 / 1.0	48	G78G (GGC→GGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (29/13);  total (32/16)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCA  >  NC_000913/14155‑14643                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ccaTTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAAGACTACAAACGCATAGCCATGAAATACCAACCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAAGTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGAATCGCAAAAACGTGCGGCATACGATCATGCTGGTCATGCGGCGTTTGAGCAAG‑GTGGTATgg                                                                                                                                                                                                                                                >  2:139663/3‑251 (MQ=255)       aaaaaaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTAACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGCATGggcggcg                                                                                                                                                                                                                                          >  1:140002/1‑251 (MQ=255)                     aGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTCTGAAGTTCTGACCGACTCGCAAAAACGAGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTggcggc                                                                                                                                                                                                                            >  1:357470/1‑251 (MQ=255)                        aaGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTggcggcggg                                                                                                                                                                                                                         >  1:177935/1‑250 (MQ=255)                           attattACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGgcgc                                                                                                                                                                                                                       >  1:108414/1‑250 (MQ=255)                                       ttttAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGat                                                                                                                                                                                                          >  1:321790/1‑251 (MQ=255)                                          tAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATAtt                                                                                                                                                                                                       >  1:347539/1‑251 (MQ=255)                                           aGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATAttt                                                                                                                                                                                                      >  1:341772/1‑251 (MQ=255)                                             gCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATAttttt                                                                                                                                                                                                    <  1:47654/251‑1 (MQ=255)                                               gTTTCCAAAACAGCGGAAGAGCATGAAATCACAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAATGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAATTCTGACCGACTCGCAAAAACGTGCGCCATGCGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTTTGGGCGGCGGCTGTTTTggcggcggagccaacttcctcggtattttttg                                                                                                                                                                                                  >  2:28957/1‑227 (MQ=255)                                                         cAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTc                                                                                                                                                                                        >  1:46451/1‑251 (MQ=255)                                                             ggAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCg                                                                                                                                                                                    >  1:175680/1‑251 (MQ=255)                                                                  agcatGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGGGAAAACGAGGCCGAGGAGAAATTTAAAGAGATCAAGGAAGCTTAAGAAGTTCTGACCGACTCGCATAAACGTGCGGCAGACGGTCAGTCTGGTCATGCTGCGTTTGAG‑AAGTGAGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGGCTTCAGCGATATTTTTAGTGGCGTTTTAGGCGacttt                                                                                                                                                                               >  2:303326/5‑247 (MQ=255)                                                                       taaaTCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCCCCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGTCTCTCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCTGTTTTGTCGGCGGCGCATACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTg                                                                                                                                                                           <  1:283773/249‑1 (MQ=255)                                                                        aaaTCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGCAAAAGAGACCGAGGCGCAATTTAAAGAGATCGAGGACGCTTATGAAGTCCTGACCGACTCGCAAAAACGTGCGGCATACGATAAGTATGGTCATGCTGCGTTTGCGCACG‑GTGGAATGGGCGGCGGCGGTTTTGGCGGCGGCGCCGACTTCAGCGATCTTTTTGGTGACGTTTTCGGCGATATTTTTggc                                                                                                                                                                         >  1:375083/1‑251 (MQ=255)                                                                            caaaaaGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGgcggcg                                                                                                                                                                     >  1:243670/4‑251 (MQ=255)                                                                             aaaaaagCCTACCAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGcggcgg                                                                                                                                                                    >  1:76097/7‑251 (MQ=255)                                                                              aaaaGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GGGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGGGACGTTTTCGGCGATATTTTTGGcggcgg                                                                                                                                                                    >  1:306000/2‑250 (MQ=255)                                                                                        cAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGAGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTATGACCGACTCGCAAAAACGAGCGGCATACGATCAGTATGGTAATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGAGCAGGCGTTAGCGTGATTTTTGGTGACGTTTTCGGCGGTGTTTTTtgcggagggaggtggtgt                                                                                                                                                          >  2:152248/1‑237 (MQ=255)                                                                                                            accaccCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGGTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCTCGGTGCt                                                                                                                                     >  1:375834/1‑251 (MQ=255)                                                                                                                cccGGACCGTAACCAGGGTGACAAAGATGCCGATGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTGCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTGGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCATCGATATTTTTGGTGACGTTTTCGGCGATACTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGAtt                                                                                                                                 <  1:161677/251‑1 (MQ=255)                                                                                                                cccGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGTGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGAtt                                                                                                                                 <  1:92435/251‑1 (MQ=255)                                                                                                                 ccGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTCAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCCGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCGGAttt                                                                                                                                >  1:339477/1‑251 (MQ=255)                                                                                                                              aGGGTGACAAAGAGGACGAGGCGAAATTCAAAGAGATCACGGCAGCTTATGAAGTACTGACCGACTCGCACAAACGCGCGGCATACGATCAGTATGGTCATGCTGCGTTTGCGCAAG‑GTGGTATGGGAGGCGGCGGTTTTGGCAGCGGCGTAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCAGACGTGGTCGTTAACGTGCATAGTGCGGTGCTGATTTAAGCCATAACATg                                                                                                                   >  1:146740/1‑251 (MQ=255)                                                                                                                                  ggACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTGCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGc                                                                                                               <  1:28957/250‑1 (MQ=255)                                                                                                                                       aaGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGGGGTCGTCAACGTACGGCGCGCGGTGCAGATTTACGCTATAACAGGGAGCTCAcc                                                                                                          >  1:138100/1‑251 (MQ=255)                                                                                                                                         gagGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCAACTCGCAAAAACGTGCGGCATACGATAAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCt                                                                                                        <  1:177621/251‑1 (MQ=255)                                                                                                                                                        tttAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTg                                                                                                                                                               <  1:204442/181‑1 (MQ=255)                                                                                                                                                        tttAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GGGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGAGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGaaggg                                                                                                                                                               >  2:204442/1‑177 (MQ=255)                                                                                                                                                                       gAAGCTTATGAAGTTCTTACCGACACGCAAAAACGTGCGGCATACGATCAGTATGGTCAGGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAaga                                                                          <  2:177935/251‑1 (MQ=255)                                                                                                                                                                        aaGCGTATGAATTACTGACCGAATCGCAAAAATGTGCGGTCTAAGATGAGTATGGTCATGCTGCGTTTGAGCAAG‑TTGGTATGGTCGGCGGCGTTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTTACGGGTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCACCGTGCGGCGCGCGGTGCTGATTTACGCTCAAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAagag                                                                         <  2:105491/251‑1 (MQ=255)                                                                                                                                                                           cTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAag                                                                           >  1:6378/1‑246 (MQ=255)                                                                                                                                                                            gttggAAGTGCTTACCGACGTGCAAAAACGGGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGCATGTGCGGCGGCGGTTGTGGCGGCGGCGCAGACGTCAGCGATATTTTTTGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGTTTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATcc                                                                     <  2:140002/247‑1 (MQ=255)                                                                                                                                                                                   gTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTcc                                                              <  2:6378/251‑1 (MQ=255)                                                                                                                                                                                       tGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGAct                                                          >  1:187073/1‑251 (MQ=255)                                                                                                                                                                                           cGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACAAAAGAGATCCGCATTCCGACTCTgg                                                      <  1:377326/251‑1 (MQ=255)                                                                                                                                                                                                               cATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCAc                                  <  1:139663/251‑1 (MQ=255)                                                                                                                                                                                                                    gATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGt                               >  1:29117/1‑249 (MQ=255)                                                                                                                                                                                                                     aTCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCACCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGACGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGc                            >  1:366732/1‑251 (MQ=255)                                                                                                                                                                                                                                   gctgcGTTGGAGCAAG‑GTGGTATGGGCGGCGGCTGTTTTGGCGGCGGCGCATACTTTAGCGATATTGTTGGGGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCGGGAAGAGTGTGACGTTTGCCACGGTAGCGTTGCAAAACCAgg              <  2:306000/251‑1 (MQ=255)                                                                                                                                                                                                                                    ctgcGTTTGAGCAAG‑GTGGTCTGGGCGGCGGCGGTTTTGGCGTCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATCTTTTTTGCGGCGGACGTTGTCGTCAACATGCGGCGCGCGGGGCTGCTTTACGCTATAACATGGAGCACCACCTCTAAGAAACTATACTGGGGGTGACCATAGAGATCCTCATTACGCCTCTGGAAAAGTGGAACTTTTGCCACAGTAGCgggtaacaccacgat             >  2:112218/1‑238 (MQ=255)                                                                                                                                                                                                                                         tttGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAAAATGGAGCTCACCCTCGAAGGAGATGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGAGGGACGTTTGCCACGGTAGCGGTGCAAAACCAGGGacac         >  1:27098/1‑250 (MQ=255)                                                                                                                                                                                                                                          ttGAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTAcaca        <  1:185087/250‑1 (MQ=255)                                                                                                                                                                                                                                          gtGAGCAAG‑GTGGGATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTAcaca        <  1:182964/249‑1 (MQ=255)                                                                                                                                                                                                                                            gAGCAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGcc     >  1:178360/1‑251 (MQ=255)                                                                                                                                                                                                                                               cAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGGCGTTGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGATGTACGTGGCGCGACACAAGGGAACCGCATACCGACACTGGACGAGAATGAACTTTGCCACGGTAACGGTGCAGAACCCGGTacaacacccgg  >  2:285128/1‑243 (MQ=255)                                                                                                                                                                                                                                               cAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGc   >  1:248395/1‑250 (MQ=255)                                                                                                                                                                                                                                               cAAG‑GTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGGGACGTTTTCGGCGATATTTTTGGCGGCGGACGGGGGCGTCAACGTGCGGCGCGCGGTGCTGCTTTACGCTATAACATGGCGCTCACCCTCGGAGAAGCTGGACGTGGCGTGACAAAAGAGATCCGCAGTCCGACTCAGGCAGGGTGTGACGTTTGCCCCGGTAGAGGTGCaaagcccaggaccagtccgcc  >  2:49571/1‑233 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAG‑GTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCA  >  NC_000913/14155‑14643

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: ‑