Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,358,812	A→G	20.6%	S126P (TCC→CCC)	puuP ←	putrescine importer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,358,812	0	A	G	20.6%	157.9 / 4.4	76	S126P (TCC→CCC)	puuP	putrescine importer
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base G (16/2);  total (48/28)
Fisher's exact test for biased strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTCTCGCTGATAAACGGCTGAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAC  >  NC_000913/1358583‑1359049                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              gcgTTCTCGCTGATAAACGGCTGAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTTCCACCCCCACCCCCGGCGGCACTTCCGGGGAGACGGGGGGGagaaagat                                                                                                                                                                                                                                      >  1:919521/1‑239 (MQ=255) gcgTTCTCGCTGATAAACGGCTGAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGGCCTTCCGGGGAGGGGGCGGGGagatagga                                                                                                                                                                                                                                      >  2:1109098/1‑237 (MQ=255)  cgTTCTCGCTGATAAACGGCTGAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACCGCGGCACTTCCGGGGAGAGGGAGGGGAGATAGGtt                                                                                                                                                                                                                                     >  2:189881/1‑239 (MQ=255)  cgTTCTCGCTGATAAACGGCTGAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACCGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGAtt                                                                                                                                                                                                                                     >  1:796184/1‑239 (MQ=255)              tAAACGGCTGAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCcaacaa                                                                                                                                                                                                                          <  1:760023/238‑1 (MQ=255)                       gAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGGTAAAATTGGGAAGAAGGGTAACACCCACACCCCCCGCCGCACTTCCGGGGAGAGGGGGGAGAGAAAGAATTTCGCCAACAACACGGTaaac                                                                                                                                                                                                                >  2:435186/1‑237 (MQ=255)                             aCCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCCCCCACACCCACGGCGGCACTTCCCGGAAGAGGGCGGAGAGATAGATTTTCGCCCACACGAAGTTAATCAACGGc                                                                                                                                                                                                          >  2:1140529/1‑239 (MQ=255)                                  aCGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCaaaaa                                                                                                                                                                                                     >  2:347464/1‑239 (MQ=255)                                            cGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATcc                                                                                                                                                                                           <  1:1032103/239‑1 (MQ=255)                                               cTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCcag                                                                                                                                                                                         <  1:347464/238‑1 (MQ=255)                                                         aTGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGAc                                                                                                                                                                               >  2:588101/1‑238 (MQ=255)                                                               tCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCCAGGCGGCACTTCCGGGGAGAGGGCGGAGAGGTAGATTTTTGCCAACAAAACGTTAATCACGGGCAAAAAGAGATAATGCAGCAGCGATGACCAGCCg                                                                                                                                                                         >  1:602754/1‑238 (MQ=255)                                                               tCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCCACAAGACGATAATAATCGGCAAAAAGAGATAAACCagcagccgatcacagacg                                                                                                                                                                         >  2:757448/1‑226 (MQ=255)                                                               tCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCg                                                                                                                                                                         >  1:1217494/1‑238 (MQ=255)                                                               tCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCg                                                                                                                                                                         >  1:112302/1‑238 (MQ=255)                                                               tCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACCAAGGGTA‑CACCCACCCCCACGGCGGCACTTCCGGGGAGAGGGCGGAGGGATAGAATTTTCCCCACAAGACGTTAATCATCCGCCAAAAGGGATAATCACACAGCGATGACCAGCCgc                                                                                                                                                                        >  2:39405/1‑237 (MQ=255)                                                               tCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGAAGCAGGGTACAACCCACACCCACCGCGGCCCCTCCGGGAAGAGGGCGGAGAGACAGATTTTCCCCCACAAGACGGTAATCATCCGGAAAAAGAGAGAATCCACCAGCGGTGCCCAGCCg                                                                                                                                                                         >  2:1187789/1‑238 (MQ=255)                                                                cccTGAACACCCAGGGAGATAAACACCACCAGTATGGAGATTTGCACCAGTACAAACAGGGTAGTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCTCCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGcc                                                                                                                                                                          <  2:381335/236‑1 (MQ=255)                                                                  cTGAACTACCAGGAAGATAAACACCACCAGGATGGAGAGTTGCAGCAGTACAAACAGGGTATTGAAGTTTGCGACCAGGTTGACGCTCTTCAGATGCGCGGCGGTTAAAATGGCGACAAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCa                                                                                                                                                                     <  1:88102/239‑1 (MQ=255)                                                                         aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCCCTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACCAGACGTTAATCATCGGCCAAAAGAGATAATCCAGCAGCGGTGACCACCCGACCATAAATc                                                                                                                                                               >  2:799738/1‑238 (MQ=255)                                                                         aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGGGGTAGGTTTTCGCCCACAAGACGGTAATCATCGGCAAAAAGAGATAATCCAGCCGCGATGAACAAcctgccctaacac                                                                                                                                                               >  2:751066/1‑227 (MQ=255)                                                                         aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGAGATAGATTTTCGCCCACAAGACGTTAATCATCGGCCAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCCTAAAta                                                                                                                                                               >  2:540554/1‑237 (MQ=255)                                                                         aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCCACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATc                                                                                                                                                               >  2:424977/1‑238 (MQ=255)                                                                         aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCcagcag                                                                                                                                                                                      <  2:369709/215‑1 (MQ=255)                                                                         aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATc                                                                                                                                                               >  1:133592/1‑238 (MQ=255)                                                                         aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGGTAATCATCGGCAAAAAGAGATAATCcagcag                                                                                                                                                                                      >  1:369709/1‑215 (MQ=255)                                                                                        aacaccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAAt                                                                                                                                               <  1:1187789/237‑1 (MQ=255)                                                                                         ccaccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCCCCCACCCCCACGGCGGCACTTCCGGGGCGAGGGGGGGGAGATAGATTTTCGCCCACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCACCCGACCATAAATCCCACGTGCGGGTTAAt                                                                                                                                               >  2:54675/1‑238 (MQ=255)                                                                                            ccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCcag                                                                                                                                                                                         <  1:113605/193‑1 (MQ=255)                                                                                            ccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCCGCCGACCATAAATACGACGTGCGGGTTAATCg                                                                                                                                             >  1:1044830/1‑237 (MQ=255)                                                                                            ccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAAACcag                                                                                                                                                                                         >  2:113605/1‑193 (MQ=255)                                                                                               tgatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAAACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCCACCATAAAACCGACGTGCGGGTTAATCGACttt                                                                                                                                        >  1:971422/1‑239 (MQ=255)                                                                                                gatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                       <  1:840848/239‑1 (MQ=255)                                                                                                gatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTA‑CACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTg                                                                                                                                      <  1:39405/239‑1 (MQ=255)                                                                                                  tgGAGATTTGCCCCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTg                                                                                                                                      <  1:918792/238‑1 (MQ=255)                                                                                                    gagaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGGCGAAGGGTACCCCCCACACCCACGGCGGCAATTCCGGGGAGAGGGGGGGAAGATAGATTTTCGCCAACAAGACGTTAATTATCGGCAAAAAGGGATAATACACCAGCGAAGACCAGCCCACCATAAAACCGACGCGAGGGGTaaagcaatttttcga                                                                                                                                   >  2:1028573/1‑226 (MQ=255)                                                                                                       aTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCCCTTCCGGGGAGAGGGCGGGGAGATAGATTTTCGCCAACAAGGCGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGACGCACATAAATCCGACGTGCGGGGTAATTGACTTTTTcgcgga                                                                                                                                >  2:60273/1‑237 (MQ=255)                                                                                                       aTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAACCCACCATAAATCCGACCTGCGGGTTAAACGACTTTTGCGCGTa                                                                                                                                >  2:268855/1‑239 (MQ=255)                                                                                                        tttGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGt                                                                                                                                 <  1:608480/237‑1 (MQ=255)                                                                                                         ttGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCCACAAGACGTTAATCAACGGGAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGGTAATCGGCTTTTGCGCGGAgg                                                                                                                              >  1:1224965/1‑239 (MQ=255)                                                                                                          tGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCCACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATTAATCCGACGTGCGGGTTAATCGACTTTTTCGCGTAGGt                                                                                                                             >  1:65294/1‑239 (MQ=255)                                                                                                                 gTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGGACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCAACCAGACGTTAATCATCGGCAAAAAGAGAAAATCCAGCAGCGATGACCAACCGACCATAAATCCGACGTGCGGGTTAATCCGCTTTTTGCCGTAGGTAAAGGc                                                                                                                       >  2:1081321/1‑238 (MQ=255)                                                                                                                           ggTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTACCCCCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAACCGACCATAAATCCCACGTGCGGGGTAATCGACTTTTTCGCGTAGGGAAAAGCCCAACCCGcc                                                                                                             >  2:582089/1‑238 (MQ=255)                                                                                                                           ggTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGGTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGcc                                                                                                             >  1:16427/1‑238 (MQ=255)                                                                                                                           ggTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGGTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGAcc                                                                                                             >  1:663483/1‑238 (MQ=255)                                                                                                                            gTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGGGGGGGGGTAGATTTTTGCCAACCAGAAGTTAATCATCGGGAAAAAGAGATAAACCAGCAGCGGAGAACAACCGACCATAAATCCGAAGTGCGGGGTAATAGACTTTTGGGCGGGGGGGTAGGCCCAACCGGccc                                                                                                            >  2:80833/1‑237 (MQ=255)                                                                                                                            gTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCCACAAGACGTTTATCATCGGCAAAAAGAGATAATCCCGCAGCGGTGACCAGCCGACCATAAATCCGACGTGCGGGTTAAATTACTTTTTCGCGTAGGTATAAGCCGAACCCGccg                                                                                                            >  2:1146722/1‑237 (MQ=255)                                                                                                                                   aGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAAc                                                                                                    <  1:435186/239‑1 (MQ=255)                                                                                                                                       aGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAAATGGc                                                                                                >  1:63744/1‑239 (MQ=255)                                                                                                                                                gTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGtt                                                                                       >  1:1219334/1‑239 (MQ=255)                                                                                                                                                 ttGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACTCCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGt                                                                                        <  1:371544/237‑1 (MQ=255)                                                                                                                                                     cGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGATAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCg                                                                                  <  1:1109098/239‑1 (MQ=255)                                                                                                                                                         ctTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGGTAATCGACTTTTTCGCGTAGGGATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTa                                                                                >  1:180140/1‑237 (MQ=255)                                                                                                                                                          ttCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGGACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGAAGACCAGCCGACCATAAATCCGAAGGGCGGGGTAATCGCCTTTTTCGCGTAGGTGTAGGCCGCACCGGCCTCCGGGAACTGGCGAACCAGTTTGCCgga                                                                                >  2:185481/1‑234 (MQ=255)                                                                                                                                                          ttCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCAACCGCAAAAAGAGATAATCCAGCAGCGATGACCAGACGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAAGCCGAACCGGCCTCCGGAAAATGGCGAACCAGTTTGCCGtc                                                                                >  2:636843/1‑235 (MQ=255)                                                                                                                                                            caggctCGCGGGGGTTAAAGTGGCGACGAAGGTTACCAGCCTCACCCTGGGCGGCACCTACGAGACGAGGGCGGAGAGATCGATTTTCGCAACCAAGACGTTTATCTTCGGCAAAAAGGGATAATCCAGCAGCGATGTCCAGCCGACCATAAACCCGATGTGCAGGTTAATCGACATTAGCGCGTAGGTATAGGCCGAGCCGGCCCCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGa                                                                           <  2:718291/234‑1 (MQ=255)                                                                                                                                                            cAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGa                                                                           <  1:423640/239‑1 (MQ=255)                                                                                                                                                                 tCGCGGGGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCg                                                                      <  2:602754/239‑1 (MQ=255)                                                                                                                                                                 tCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCg                                                                      <  1:60273/239‑1 (MQ=255)                                                                                                                                                                    cggcggTTAAAATGGCGACGAAGGGTACCACCCCCACCCACGGCGGCCCTTCCGGGGAGGGGGGGGGGAGGTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGGTAATCGACTTTTGCGCGTAGGGATAGGCCCAACCGGCCTCCGGAAACCGGCGAACCCGTTTGCCGTTGCCGATAGCggga                                                                   >  2:581396/1‑237 (MQ=255)                                                                                                                                                                     ggcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGCTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTa                                                                   <  2:493277/238‑1 (MQ=255)                                                                                                                                                                        ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACa                                                               <  1:799738/239‑1 (MQ=255)                                                                                                                                                                        ggTAAAAACTGCGACGAAGGTTACCCCCCACCCCCGCGGCGGCACTTCCGGGAAGACGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTGGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACa                                                               <  2:112446/239‑1 (MQ=255)                                                                                                                                                                                         aGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACTGGCGAACCAGTTGGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgc                                               <  1:636843/238‑1 (MQ=255)                                                                                                                                                                                                     acacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCg                                                                      >  1:563652/1‑203 (MQ=255)                                                                                                                                                                                                     acacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCg                                                                      <  2:563652/203‑1 (MQ=255)                                                                                                                                                                                                     acacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGAc                                  >  1:680375/1‑239 (MQ=255)                                                                                                                                                                                                            cggcggCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGa                            >  1:83936/1‑238 (MQ=255)                                                                                                                                                                                                             ggcggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAAc                          >  2:649530/1‑239 (MQ=255)                                                                                                                                                                                                                 gctcTTCCGGGAAGAGGGCGGCGAGATAGATTTTCGCCAACAAGACGTTCATCATCGGCTAAAAGAGATAATCCAGCAGCGATGACAAGCCGACCCTAACTCCGACGTGCGGGTTAATCGACTTTTGCGAGTAGGTATAGGCCGAACCATCGTCCGGTAACTGGCGAACCAGGTTGCCGTAGCTGATAGCGGTAAACCGCACACCCCCCAGACACAGCAACTACAACGCCGGAAAGTg                       >  2:959077/4‑238 (MQ=255)                                                                                                                                                                                                                   ccTTCCGGGAAGAGTGCGAGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGAAAAAAGAGATAACCCAGCAGCGATGACCACCCGACCATAATTCCGACGGGCGGGTTAACCGGCTTTTGCGCGTAGTTATAGGCCGACCCGTCCGCCGGAAACTGGCGAACCAGTTTGCCGTAGCGGATAGCGGTAAACAGCCCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGAcc                    <  2:919521/238‑1 (MQ=255)                                                                                                                                                                                                                   aCTTCCGGGAAAGGGGCGAGGAGAAAAATTTCCGCAAAAAAGACGTTAACCATCGGAAAAAAGAGATAATCAAGCGGCGATGACCAGCGGACAATAAATCCGACGTGCGGGTTAACCGCCTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACGGGAGAACAAGTTTGCCGTAGCTGATAGGGGAAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGAcc                    <  2:1224965/239‑1 (MQ=255)                                                                                                                                                                                                                           gAAGAGGGGGGAAGAAAAAATTTCCGCAAAAAAGACGTTAATAATCGGAAAAAAGAGAAAACCCAGAACGGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTa             <  2:1219334/238‑1 (MQ=255)                                                                                                                                                                                                                                gggCGGAAAGATCGATTTTCGCCAACAAGACGTTAACCATCGGCAAAAAGAGAAAATCCAGCAGCGATGACCAGCCGACCATAAATCCGCCGTGCGGGTTAACCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAACCGGGCGAACCAGTTTGCGGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGcc        <  2:83936/238‑1 (MQ=255)                                                                                                                                                                                                                                 ggCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCgg      >  1:907254/1‑239 (MQ=255)                                                                                                                                                                                                                                 ggCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCgg      >  1:807457/1‑239 (MQ=255)                                                                                                                                                                                                                                      agagaTAGATTTTCGCCAAAAAGACGTTAATCATCGGAAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc  <  1:424977/238‑1 (MQ=255)                                                                                                                                                                                                                                      agagaTAGATTTTCGACAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAACCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGAAAACTGGCGAACCAGTTGGCCGTAGCTGATGGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc  <  1:54675/238‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              GCGTTCTCGCTGATAAACGGCTGAAGTGACCAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAC  >  NC_000913/1358583‑1359049

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