Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,223	G→A	R19H (CGT→CAT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,223	0	G	A	97.4%	105.3 / ‑4.3	38	R19H (CGT→CAT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (24/13);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATA  >  NC_000913/13975‑14468                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       caACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGGCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTg                                                                                                                                                                                                                                                      <  2:432586/251‑1 (MQ=255)       gCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCGATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATc                                                                                                                                                                                                                                                 <  2:667126/250‑1 (MQ=255)         ccAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGCCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTAGAAACGGGTAATTATACTGCCACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGAAGGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCaa                                                                                                                                                                                                                                               >  2:269623/1‑249 (MQ=255)              agcagACTGCCGGTGCTGGTGCTTCTGCAAACATCGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCCATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGc                                                                                                                                                                                                                                         <  1:125886/251‑1 (MQ=255)                                aTGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCAt                                                                                                                                                                                                                       >  1:66729/1‑251 (MQ=255)                                                          gACGATATTGTCGACGCTGAATTTGAAGAAGTCAACGACAAAAAATAATCGCCCTATAAACGGGTAATTAAAATGACACGGGAGAAGGGGAATTTCATCTCCGCCCGTGCATTAATCTCGGGGCAATTTAAAAAAAATGGCTAAGAAAGATTATTACGAGATTTTAGGCGTTTACAAAACAGCGGAAGAGCATGAAATCAAAAAGCACTACAACAGCCTGGCCATGaacaaacacccgaaccaaaacaagg                                                                                                                                                                                             >  2:437208/1‑228 (MQ=255)                                                              aTGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTACTTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTACGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAACTCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGAAcctacaccaggatg                                                                                                                                                                                          >  2:162009/1‑238 (MQ=255)                                                                gttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAAGCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGATATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGAc                                                                                                                                                                                        <  2:281858/250‑1 (MQ=255)                                                                       aCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCTGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGc                                                                                                                                                                                <  1:469094/251‑1 (MQ=255)                                                                             aaTTTGAAGAAGGCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAATGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCGCCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGc                                                                                                                                                                          <  2:635245/251‑1 (MQ=255)                                                                               tttGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCg                                                                                                                                                                         >  1:150789/1‑250 (MQ=255)                                                                                     gaagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAATAAGGCCTAAAAACGCCTGGCCATGAAATACCACCCGGACCGTAACAAGGGAGACAAAGAGGCCGAGGCGAAAttt                                                                                                                                                                   >  2:587434/1‑250 (MQ=255)                                                                                         tCAAAGACCAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGTGCGAAGGGGAATTTCCTCTCAGCCCGTGCTTTCATCTAGGGGCAATTTAAGAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAACCAGCGGAAGAGCATGAAATCAACAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCAGAGGCGAAATTT‑aaa                                                                                                                                                               >  2:552910/1‑249 (MQ=255)                                                                                          cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACTCGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTACTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCAGGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGTCCGAGGCGAAATTT‑TAagag                                                                                                                                                            <  2:50874/251‑1 (MQ=255)                                                                                                                      cGGGTAATTATACTGACACGGGCGAAGGGGCATTTCCTCTCCGCCAGTGCATTCATCTCGGGGCAATTTACAAAAGATGGCTAAGCCAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGTCCGTAACCAGGGGGAACAAGAGGCCGAGGCGACATTTCAAAG‑GACCATGGAAGCTTCTGAAGTGCTGACCg                                                                                                                                >  2:306945/1‑251 (MQ=255)                                                                                                                                    gACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGt                                                                                                                  >  1:417677/1‑251 (MQ=255)                                                                                                                                        cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCgg                                                                                                              >  1:600350/1‑251 (MQ=255)                                                                                                                                           gcggaGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAACAGCTGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAACAGTCCTCCAAACGCCTGGCCATGAAATACCACCCGTACCGTAACCAGGGCGACAAAGCGGCCGCGGCGAAATTT‑AAAGAGACCAAGGAAGCTTATAAAGTTATGACCGACTCGCAAACACTTGCGGCa                                                                                                            >  2:439230/5‑250 (MQ=255)                                                                                                                                             gAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGTTGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGACATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGGTCAAGGAAGCTTATGAAGTTCTGACCGACTCGCACAAACGTGCGGCATAc                                                                                                         >  2:71381/1‑251 (MQ=255)                                                                                                                                             gAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATAc                                                                                                         >  1:445905/1‑251 (MQ=255)                                                                                                                                                                           cATCTAGGGGCAATTTAAAAAAGATGGCTAAACAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGCAATCACAAAGGCCTACAAACCCCTGGCCATGCAATACCACCCGGACCGTAACCAGTGTGACAAAGCGGCCGAGGCGCAATTT‑AAAGAGATCACGGAAGCTTATGAAGTTCTGCCCGACTCGCACAAACGTGCGGCATACGATCCGTATGGTAATGCTGCGTTTGAGcc                                                                            >  1:77937/1‑249 (MQ=255)                                                                                                                                                                                aGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAg                                                                          <  1:587434/247‑1 (MQ=255)                                                                                                                                                                                 ggtgCAAGTTAAAAAAGATGGCTAAGCAAGATTAATACGAGATTTGAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTggt                                                                     <  2:261512/248‑2 (MQ=255)                                                                                                                                                                                             aaaaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGggcggcggc                                                         <  1:319034/251‑1 (MQ=255)                                                                                                                                                                                               aaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGcggcggcgg                                                       <  2:107992/251‑1 (MQ=255)                                                                                                                                                                                                aGATGGCTAAGCAAGATGATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGt                                                      <  2:322042/251‑1 (MQ=255)                                                                                                                                                                                                     gCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGCGATCAAGGCAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGGGGTATGGGCGGCGGCGGtttttg                                                 >  2:668959/1‑249 (MQ=255)                                                                                                                                                                                                        aaGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTggcgg                                              >  1:460353/1‑251 (MQ=255)                                                                                                                                                                                                           cAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAATCGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGAAGACAAAGAGGCCGAGGAGAAATTT‑AAAGAGATAAAGGACGCTTCTGCAGTTCTTACCGACTCGCAAAAACTGGCGGCTAACGCTCATAATGGTAATGCTGCGTTTGAGCAAGTAGGTATGGACGGCGGATGTTTTgtcggcgg                                           >  2:594888/1‑251 (MQ=255)                                                                                                                                                                                                            aagcttattaCGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTggcggcggg                                          >  2:501538/5‑250 (MQ=255)                                                                                                                                                                                                                         gaTTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCg                             >  1:206841/1‑251 (MQ=255)                                                                                                                                                                                                                               aGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGAGGTAttt                       >  2:146322/1‑251 (MQ=255)                                                                                                                                                                                                                                      cccAAAACAGCGAAAGAGCATGAAACCAAAAAGGCCGACAAACGTCTGGCCATGAAATACCTCCCGGACCGTAACCAGTGTGACAAAGAGGTCGAGTCGAAATTT‑AAAGCGATCAAGTAATCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGGTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCTTAGACTTCATCGATATTTTTGGTGa                <  2:398963/250‑1 (MQ=255)                                                                                                                                                                                                                                         aaaaCAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGt             >  1:616821/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCg     >  1:286495/1‑250 (MQ=255)                                                                                                                                                                                                                                                    agagCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGTCCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATTCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGTCGata  <  1:243486/251‑1 (MQ=255)                                                                                                                                                                                                                                                    agagCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGTAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGGTATTTTTGGTGACGTTTTCGGCgct   >  1:240593/1‑248 (MQ=255)                                                                                                                                                                                                                                                    agagCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGACATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGAGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGGTTTCGGCGatt  >  2:175303/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       CAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTT‑AAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATA  >  NC_000913/13975‑14468

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