Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,581,316	G→A	G60E (GGA→GAA)	yhhY →	N‑acetyltransferase YhhY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,581,316	0	G	A	97.2%	217.3 / NA	72	G60E (GGA→GAA)	yhhY	N‑acetyltransferase YhhY
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (35/35);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGG  >  NC_000913/3581071‑3581562                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       tGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGAcg                                                                                                                                                                                                                                                   >  1:439064/1‑251 (MQ=255)      ttGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGAcgtcgt                                                                                                                                                                                                                                               <  2:152384/250‑1 (MQ=255)        gACTCGCATTTGAAAGGCAATGATGGTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGGTAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTgg                                                                                                                                                                                                                                             <  2:462245/250‑1 (MQ=255)            cGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCCGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGAGCCTAATCATTCCGATCATATGTGTCAGGAGCG‑AATCGCCGTTCGGCCCGGCATCAAGCAACTAGTAGCATGTATTGATGAAGACGTAGGGGGCCac                                                                                                                                                                                                                                        >  2:390190/1‑250 (MQ=255)                   gAAAGTCAATGATGTTGAAAGGGACATTTACCCCACAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCCGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGAG‑ACTCGACGGTCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGCCGTCGTGGGCCATCTCAcc                                                                                                                                                                                                                                  >  2:74733/1‑250 (MQ=255)                   gAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCAcc                                                                                                                                                                                                                                  >  1:80483/1‑250 (MQ=255)                    aaaGTCCATGAAGTTTAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGACAAGCAGAAACAAGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGAGGATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCTTAAGTGGCAGGAACG‑CCTCGCCGATAGTCCCGGCACAAAGCAACTCGTCTCGTGTATTGATGAAGACGTCGGGGGcagacccacaat                                                                                                                                                                                                                                >  2:492521/1‑240 (MQ=255)                     aaGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACACAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCAtt                                                                                                                                                                                                                               >  1:154201/1‑251 (MQ=255)                       gTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGCTAGTAATCCGCCACGCCGAAACACGGGATAACGAGGCCATCAGGCAGATTCACGCCCAGCCGACGGTGAATTGCAACACACTAC‑AGGTGCATCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCCCCGATCGTCACGGCCTCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGGGGGCCCTCTCACCATTg                                                                                                                                                                                                                              >  2:364766/1‑250 (MQ=255)                         cAATGATGTTGAAAGGTACATTTACCCCCAAGCGGACAAAGGAATGAATGCGATAGTAATACGCCAAGCAGAAACACGGGAATACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTCTTGCAACACACTCC‑AGGTGACTCATCCTCCCGATCATATGTGGCAGGCGCG‑TATCGCCGCTCGTCCCGGCACCATGCAACCCGTCGCCTGTATTTCTGGAGCCGACCTGGGCCATCACCCCATTTACg                                                                                                                                                                                                                           >  2:295325/1‑251 (MQ=255)                            tgatgTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGc                                                                                                                                                                                                                        >  1:66987/1‑251 (MQ=255)                             gatgTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCCCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTTACGTGc                                                                                                                                                                                                                        >  2:78197/1‑250 (MQ=255)                                    aaaGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGACCAGCCTGAGGTGTATTGCAACACACGACCAGG‑GCCTCATCCTTCCGATCATATGTGGCAGGAGCGCAC‑CGCAGGTCGTCCCGGCATCAATCAACACGTCGCCTTGATTGATGAAGGCGGCCTGGGCCAGCTCAcccttaccgggcaaacgcg                                                                                                                                                                                                                 >  2:314080/1‑233 (MQ=255)                                            aattACCGCAAAGAGGAAAAAGTAATGAGTGATATAGTAATACGCCCAGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCACCACACTAC‑AGGTGCCGCATCCGTCCGATCATATGTGGCATGAGCG‑TCTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGGCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAcgccg                                                                                                                                                                                                        <  2:328947/249‑1 (MQ=255)                                                      aaGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACCCTAA‑AGGTGCCTCATCCTTCCGCTCATATGTGGCAGGAGCG‑CCTCGCCGCTCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACCTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTccagtt                                                                                                                                                                                              >  2:202214/1‑246 (MQ=255)                                                                aGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTg                                                                                                                                                                                    <  1:171566/251‑1 (MQ=255)                                                                 ggAATGAGTGAGATAGTAATACGCAACGCAGAAACACGGGATTACACGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACCCTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGACCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGACGATTTTgg                                                                                                                                                                                   >  2:419119/1‑251 (MQ=255)                                                                      gagtgagATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCt                                                                                                                                                                              >  1:382385/1‑251 (MQ=255)                                                                           agaTAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCCTCTCACCATTTGCGGGCAACAGCCCCCCCGCCGCAGTCATGTTGCCGGTTTTGGTATCtggttc                                                                                                                                                                         >  2:422944/1‑247 (MQ=255)                                                                             aTAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCg                                                                                                                                                                        >  1:346274/1‑250 (MQ=255)                                                                                atAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGActc                                                                                                                                                                    <  1:74733/250‑1 (MQ=255)                                                                                 tAATACGCCATCCCGAAACACGGGATTACGAGGCCATCAGGCACATTCACGCCCAGCCGGAGGTGCATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGCTCATATGTGGCAGGCGCG‑ACTCGCCGATCGCCCCGACATAAAGCAACTCGTCGCCTGCTTTGATGAAGCCGTCGTAGGACCTCTCACAATTGACGACCGACAGCGCCCACGACGCCGGCATGTTGCCGCATTTGTAATCTGTGTcacttcc                                                                                                                                                                   >  2:243842/1‑245 (MQ=255)                                                                                    tACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGt                                                                                                                                                                >  1:152355/1‑251 (MQ=255)                                                                                     ccGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACTCCCAGCCGGAGTTGTATTGCAACACACTAC‑AGGTTCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGACCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGtt                                                                                                                                                               <  2:309755/250‑1 (MQ=255)                                                                                         cACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGga                                                                                                                                                            >  1:150195/1‑250 (MQ=255)                                                                                                       ggATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAc                                                                                                                                                                                                            >  1:119868/1‑188 (MQ=255)                                                                                                       ggATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCg                                                                                                                                             <  1:534417/251‑1 (MQ=255)                                                                                                       ggATGACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAc                                                                                                                                                                                                            <  2:119868/188‑1 (MQ=255)                                                                                                          ttccGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGGGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑CCTCGCCGATCGTCCCGGCCTCAAGCCACTCGTCGCCTGTATTGATGACTACGTCGGGGGCCATCTCCCCATTGACGTGCAACCACGCCCCCGCCGCAGTCATGTTGCCGATTTTGGAATCTGTGTGTACTCTCGTTGGAGGAACCGCGGCGTCGccc                                                                                                                                          >  2:178734/4‑250 (MQ=255)                                                                                                           tACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGTTCATATGGGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCa                                                                                                                                          <  1:492521/250‑1 (MQ=255)                                                                                                             cGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAgcg                                                                                                                                       <  1:525530/251‑1 (MQ=255)                                                                                                              gAGGCCATCAGGCAGATTCACACCCAGACGGAGGAGTATTGCAACACACTAC‑AGGGGCATCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATAGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAATACGGCATGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTTCCGATTTTGGTTACTGTGCAGACTCTAGTTGGAAGAACAGCGGCGGCGCCAgcgc                                                                                                                                      >  2:460751/1‑251 (MQ=255)                                                                                                                  ccATCAGGCAGTTTCACGCCCAGGCGGAGGTGTATTGCAACACACTAC‑AGGGGCCTCATCCTGCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGCCGTGCAACAGCGCCCACGCCGCAGTCATGGTGCCGATTTTGGTATCTGTTTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTg                                                                                                                                  <  2:328374/251‑1 (MQ=255)                                                                                                                      cAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGTAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATg                                                                                                                               <  1:58359/250‑1 (MQ=255)                                                                                                                      cAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCTTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGc                                                                                                                              <  1:55248/251‑1 (MQ=255)                                                                                                                         gCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgag                                                                                                                           <  1:78197/251‑1 (MQ=255)                                                                                                                         gCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCg                                                                                                                             >  1:85744/1‑249 (MQ=255)                                                                                                                           aGATTCACGCCCAGCCGGGGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCAGCTCACCATTGACGTGCAACTGCGCCCACGCCGCATTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgagag                                                                                                                         <  1:335609/251‑1 (MQ=255)                                                                                                                                cACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatgat                                                                                                                    <  1:202214/251‑1 (MQ=255)                                                                                                                                 ccGCCAAGCCGGCGGTGTTGTGCAACACCCTAA‑AGGTCCCACATCCGTCCGTACATTGGTGGCAGGTGCG‑ACTCGCCGAAAGTCCCGGCATCAAGCAACTGGATGCTTGTATTGATGAAGACTTCGTTGGCCATCTCGCCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGAGTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGAtt                                                                                                                   <  2:131067/250‑1 (MQ=255)                                                                                                                                      cAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAt                                                                                                              >  1:501491/1‑251 (MQ=255)                                                                                                                                               gtgtCTTGCAACACACTAC‑AGGTGCCCCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATGGATGAAGACGTCGGGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACaa                                                                                                     <  2:152355/251‑1 (MQ=255)                                                                                                                                               gtgtATGGCAACACCTTAC‑AGGTGCCTCATCCTTCCTAGCATATGTGGCAGGAGCG‑ACTCGCCTATCGTCCCGGCATCAAGCAACTCGTCTCCTGTATTGATTAAGACGTCGTGGGCCATCGCACCATTGACGTGCAACAGCGCCCGCGCCGCAGTCATGTTGCTGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGAGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACaa                                                                                                     <  2:421722/251‑1 (MQ=255)                                                                                                                                                 gtATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACt                                                                                                   <  1:404471/251‑1 (MQ=255)                                                                                                                                                   aTTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTg                                                                                                  <  1:422944/250‑1 (MQ=255)                                                                                                                                                      gCAACCCACTAC‑AGGTGCCTCAACCTCTCGATCATATGTGGCAGGAGCG‑ACTCGCCGAGCGTCCCGGCATCAAGCAACTCGGCGCCTGTATTGATTAAGACGTCGTGGGCCATCTCTCCATTGACGTGCACCAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTg                                                                                              <  2:501491/251‑1 (MQ=255)                                                                                                                                                      gCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCGTCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTg                                                                                              >  1:145660/1‑251 (MQ=255)                                                                                                                                                       cAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGc                                                                                             <  1:178653/251‑1 (MQ=255)                                                                                                                                                          cacacTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGACGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCAGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAACTGTGCGACAACTGGTTGCggg                                                                                          >  1:528345/1‑251 (MQ=255)                                                                                                                                                                aC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCAACTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTATCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGAt                                                                                     <  1:314080/250‑1 (MQ=255)                                                                                                                                                                    ggtgcctaatgtttCCGGTCATTTGTGACAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGAAACTACTGGCCTGTATTGACGAAAACGAACTGGGCCATCTCACCCTTGACGTGCAACAGAGCCCCCGCCGCAGTCAAGTTGCCGACTTTGGTATATGTGTTGAACATCGTTGGAAGAACAGCGGAGTCGCCAGCTCCATGATGCGGGAAGTGTTTGAAATGTGCGATCACTGGTTGCGGGtgagtcgc                                                                                  >  2:360961/13‑243 (MQ=255)                                                                                                                                                                    ggTGCCTCATCCTTCCGATCAGAGGGGGCAGGAGCT‑TCTCGCCGATCGTCCCGGCATCAAGCAACGCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCAGTGACGTGCAAAAGCGCCCACGCCGCAGTCATGTTGCCGATTGTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCa                                                                                 <  1:32328/251‑1 (MQ=255)                                                                                                                                                                     gTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCAt                                                                                <  1:422852/251‑1 (MQ=255)                                                                                                                                                                           catccgcccGATAATATGTGGCAGGAGAG‑CCTCGCCGATCGGCCCGTCATCAAGCAACGCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTTACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGaa                                                                            <  2:439064/242‑1 (MQ=255)                                                                                                                                                                                  ccGATCATAGGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCgtg                                                                   <  1:530920/251‑1 (MQ=255)                                                                                                                                                                                   cgatcagttgtgGCAGTAGCG‑ACGCGCCGATCGTCCCGGCATCAATCAACTCGTCGCCTGTATTGAGGAAGACGTCGTGGGCCATCTCCCCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCgtgt                                                                  <  2:198980/243‑1 (MQ=255)                                                                                                                                                                                    gATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGtt                                                                 >  1:45885/1‑251 (MQ=255)                                                                                                                                                                                     aTCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGttt                                                                <  1:364766/251‑1 (MQ=255)                                                                                                                                                                                       cATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGCACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGCTTAAAATGTGCGCCAACTGGTTGCGGGTCGATc                                                                                    >  2:153428/1‑229 (MQ=255)                                                                                                                                                                                       cATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATc                                                                                    <  1:153428/229‑1 (MQ=255)                                                                                                                                                                                       cATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCATGTTTg                                                               <  1:417154/250‑1 (MQ=255)                                                                                                                                                                                        ttatGTGGCAGGCGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCTTCTCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCTACTCTCGTTGGCAGAACCGCGTCGTCGCCAGCGCCCTGATGCGTGAGATGATTGTAATGTGCGTCACCTGGTTGCGGGTAGATCTCCTTGTACTTACCGTGtttttc                                                             >  2:273358/2‑248 (MQ=255)                                                                                                                                                                                            gtgGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGAt                                                          <  1:292935/250‑1 (MQ=255)                                                                                                                                                                                                         aCTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGc                                               <  1:36917/249‑1 (MQ=255)                                                                                                                                                                                                         aCTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTAATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTAGCCAGCGACCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCCTTGGACTAAACGTGTTTGTCGATACCGCCCCGGc                                               >  2:36917/1‑249 (MQ=255)                                                                                                                                                                                                                 aTCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTc                                     >  1:128016/1‑251 (MQ=255)                                                                                                                                                                                                                     tCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGGGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCGTTGAACTAACCGTGTTTGGCGATAACGCCCCGGAAATTAAGGTCtat                                  >  2:451761/1‑250 (MQ=255)                                                                                                                                                                                                                           cATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTGGCGGGTAGATCGCATTGAACTAACCGTGTTAGTCGATAACGCCCCGGCAATTAAGGTCTATaaaaaa                            >  2:170102/1‑250 (MQ=255)                                                                                                                                                                                                                                        gtcgCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAAt              <  1:486382/251‑1 (MQ=255)                                                                                                                                                                                                                                             cTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGATGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGAATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATCAAAAATACGGATTTGAAATTGa           >  1:530740/1‑249 (MQ=255)                                                                                                                                                                                                                                                tATTGATGAAGACGTCGTGGGCCATCTCACAATTGACGTGCAACAGCGCCCCCGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGTAATTGAAgga       >  1:395936/1‑249 (MQ=255)                                                                                                                                                                                                                                                    gATGAAGACGTAGTGGACCCTCTCACCATTGACGTGCAAAAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGCAATGTGCGACAACTGGTTGCGGGCAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCACGGCAATTACGGTCTATAAAAAATACGGATTTGAAATTGAAGGGACTgg  >  2:291874/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       TGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTAC‑AGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCG‑ACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGG  >  NC_000913/3581071‑3581562

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