Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,956	G→C	L263L (CTG→CTC)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,956	0	G	C	90.0%	107.7 / 6.7	40	L263L (CTG→CTC)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base C (13/23);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

CTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCT  >  NC_000913/14713‑15198                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    cTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCg                                                                                                                                                                                                                                              >  1:126426/1‑251 (MQ=255)            gTCAGGGCCGCGGTACTTTGAGCAAAGACCCGTGCAACAAATGTCATGGTCATTGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGAGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGTGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCTGATTTTCTAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGa                                                                                                                                                                                                                                     <  1:217819/249‑1 (MQ=255)                          acgcttgtCAAAGACCCGGGCAACAAATGTCATGGTGAGGGTCGTGTGGAGCGCAGCAAAACGCTGTCCGGTAAACTCCCGGCAGGGGTGTGCTCTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCCTGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCACTAACCTCTATTGCGAAGTCCCGATCAACTTCGCTaggg                                                                                                                                                                                                                      <  2:248395/243‑4 (MQ=255)                                    aaGATCCGTGCAACAAATGTCATGGTCATGGGCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGTAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGt                                                                                                                                                                                                           <  1:9785/251‑1 (MQ=255)                                     aGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTg                                                                                                                                                                                                          <  1:222517/251‑1 (MQ=255)                                           gtggaacaaATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAGTCCCGGCAGGGGTGGACACTGGAGACCGCGTCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGaaa                                                                                                                                                                                                    <  1:178885/247‑1 (MQ=255)                                             gCAACAAATGGCAGGGTCATGTTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACTTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACTTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAAt                                                                                                                                                                                                   <  1:330311/250‑1 (MQ=255)                                                  aaaTGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGCGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGacgt                                                                                                                                                                                             >  1:232841/1‑248 (MQ=255)                                                    aTGTCATGGTCATGGTCGTGTTGAGCGCAGCATAACGCTATCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGGACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGACGGCAACAACCTCTATTGCGAAGTCCCGAGCAACTTCGATATGGCGGCGCAGGTTGGCGAGATCGAAGTAc                                                                                                                                                                                           >  1:140252/1‑251 (MQ=255)                                                    aTGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATATGTACGTTCAGGTTCAGGTTAAACAGCACCAGATTTTCGAGCGTGAAGGCAACAACCTCTATAGCGAAGTCCCGATCAACTTCGCTATGGCGGCCCTGGGTGGCGGAATCGacatac                                                                                                                                                                                           >  1:8338/1‑246 (MQ=255)                                                                   acgggtttaGCGCAGAAAAACGGGGTCCGTGAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGAGCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTGGATGGt                                                                                                                                                                             <  2:258449/242‑1 (MQ=255)                                                                          gAGCGCAGCAAAACGCTCTCCTGTAATATCCCGGCTGGGGTGGCCACTGGTGACCTCATCTGTCTTGCGGGCGAAGGTGTAGCGGGCGAGCTTGGCGCACCTGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACTGCACTCGATTTTAGATCGTGAAGGCAACAGCCTCTATTGCGAAGTCCCGTTCGTCTTCGCTATGGCGGCGCTGGGTGGCGAAATAGAAGTACCGACCCTTGATGGTCGCGtga                                                                                                                                                                      <  2:321790/250‑3 (MQ=255)                                                                                agcaAAACGCTGTCCGTTAAAATCCAGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGAGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGGTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTACCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATAGACATACCAGCCCTTGATGGTGGCGTCAAACTGa                                                                                                                                                                >  1:241550/1‑250 (MQ=255)                                                                                          cgccctgtatgaTCCCGGCAGTGGGGGACACTGGAGACCGCATCCGGCGTGCTGGCGGAGGTGCAGTGGGCGTGCAGTGCGCACCGTCATGCGATCTGTACGTTCATGTTCAGGTTAAACGGCACCCGTTTGTCGATCGTGAAGGCAACATCCTCTAGGGCGAAGTCCCGATCAACTTCGATATGGCGGCGCTGTGTGGCGAAATCGAAGTACTGACTCTTGATGGTCGCGTCAAACTGAAAGTGcccgga                                                                                                                                                     <  2:243670/240‑5 (MQ=255)                                                                                           gTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAAAAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGc                                                                                                                                                     <  1:86521/250‑1 (MQ=255)                                                                                                ttAAAATCCAGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGCGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTAAGGCGGCACTGGGTGGCGGAATCGGAGGAACGACCCTTGCTGGGCGCGCCCAACTGCAAGTGGCTGGCggaaac                                                                                                                                               >  2:167240/1‑246 (MQ=255)                                                                                                            cAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAg                                                                                                                                   <  1:207327/251‑1 (MQ=255)                                                                                                             aGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGCTTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGc                                                                                                                                  >  1:350648/1‑251 (MQ=255)                                                                                                                              aCCGCATCCGGCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGt                                                                                                                 <  1:29701/251‑1 (MQ=255)                                                                                                                               ccGCAGCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGAATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTa                                                                                                                <  1:373753/251‑1 (MQ=255)                                                                                                                                    cccGTCTGGCGGGCGAAGGTGAAGTGGGCGAGCATGGCGACCCGGGAGGCGAGCGGTACGGTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCTACAACCTCGATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGTCCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGTTAAGCTTTTCCGTATGCGCGGTAAAGGc                                                                                                           <  2:336343/250‑1 (MQ=255)                                                                                                                                       gtcgtGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGCGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGAGAAACCCAGAACAGAAAGCTTTTCCGCATGCGCGGGAAAGGCGTc                                                                                                        >  1:182755/5‑251 (MQ=255)                                                                                                                                           ggCGGGCTAGGGTGAAGTGGGCGAGCAGGGCGTACCGGAAGGCGATCTGGTGGTTCAGTTGCATGTTAAACAGCACCCGATTGTCGAGCGTGAAGTCAACTACCTCTATTTCGAAGGCGCGATCAATTTCGTAATGGCGGCGCTGTGTGTCGAAATCGAAGTGCCTACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAATCCCTGACCGGTAATCTTTTCGGTATGCTCCGTAAATGCGTCAAg                                                                                                     <  2:339477/249‑1 (MQ=255)                                                                                                                                                        gAAGCGGGCGAGCAGGACGCATCGGCAGGCGATCTGTACGTTCAGGTCCAGGTTTAACTGCTCCCGATTCTCGAGCGTGACGGCAACAACCTCCATTGCGATGTCCCCATCAACTTCGCTATGGCGGCGCTTGGTGGTGAAATCGAAGTACCTACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTGCCGTATGCGCGGTAAAGGCGTCAAGTTTTTCCGCGGTgg                                                                                       <  2:95927/251‑1 (MQ=255)                                                                                                                                                          aGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGgcg                                                                                     >  1:217816/1‑251 (MQ=255)                                                                                                                                                                     aTGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGAt                                                                          >  1:294729/1‑251 (MQ=255)                                                                                                                                                                                           gTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTgtcg                                                      <  1:167240/248‑1 (MQ=255)                                                                                                                                                                                           gTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGCCCGGTAAACTTTTCCGGATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGtcgtcg                                                   >  2:148881/1‑251 (MQ=255)                                                                                                                                                                                           gTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTgtcgtc                                                    <  1:151267/250‑1 (MQ=255)                                                                                                                                                                                              cGTTCAGGTTCAGGTTAAACAGCCCCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGAGAAGGCCCGATCAACTTCGCTATGGCGGCGATGGGTGGAGAAATCGAAGGACCGACACTTGATGGTAGCGTCAAACTGAAAGAGCCTGGCGAAACCCAGCCCGGTAAGATTTTCCGGATGCGCGGTAAAGGCGTCCAATCTGTCCGCGGGGGAGCACAGGGTGTG‑TTGTTGTGCcgcggttgagggtaa                                                 >  2:221192/1‑239 (MQ=255)                                                                                                                                                                                                       tcaggttAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGACGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGt                                         >  1:130468/1‑249 (MQ=255)                                                                                                                                                                                                         aggttAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGc                                     <  1:141862/251‑1 (MQ=255)                                                                                                                                                                                                           gttAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGGAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTTTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGAGCCGCGTTGTCGTCGAAACACCGGTAGGCCt                                   >  1:329363/1‑251 (MQ=255)                                                                                                                                                                                                               aaCAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAAc                               >  1:326848/1‑251 (MQ=255)                                                                                                                                                                                                                 cAGCACCCGATTGTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACGTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTTCCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAG‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGa                             <  1:41382/251‑1 (MQ=255)                                                                                                                                                                                                                  cgaaaccgatgtccgggCGTGGAGGCTACAACCTCGATTGCGAAGTCCCGAGCAACTACGATATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGGCAAGTCTGTCCGTGGTGGCGCACAGGGTGAT‑TTGATGTGCAGCGTTGTCGTCGAAACTCCGGTAGGCCTGGACGag                            <  2:82556/235‑2 (MQ=255)                                                                                                                                                                                                                                gAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAgct              <  1:162003/251‑1 (MQ=255)                                                                                                                                                                                                                                      gatggCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGAGCTGGGGGGGGGAATCGAAGTACCGTCCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAAACCGGATAGCTATTCCGGATGGGCAGTAAAAGAGGCCAGTCTGGCCGCGGGGGGGCACCGGG‑GAGGTTGCTGGGCCGCGTTGTCGTCTAAACACCGGTAGGCCGGTACGAAAgggaggaaaaagtgcggc         >  2:190544/4‑234 (MQ=255)                                                                                                                                                                                                                                            aacaacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagct  <  1:148881/251‑1 (MQ=255)                                                                                                                                                                                                                                            aacaacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAAGCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagct  <  1:95063/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    CTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCT  >  NC_000913/14713‑15198

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