Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,708,814	A→G	9.8%	A632A (GCA→GCG)	rsxC →	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,708,814	0	A	G	9.8%	144.4 / 7.6	61	A632A (GCA→GCG)	rsxC	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base G (4/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

GCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGC  >  NC_000913/1708527‑1709093                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                          gCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCa                                                                                                                                                                                                                                                                                           <  1:104638/288‑1 (MQ=255)      aCAGGCAGAACGGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                                                                                                                                                                                                                                                                                     <  1:304279/289‑1 (MQ=255)          gCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                                                                                                                                                                                                                                                                 >  1:494210/1‑289 (MQ=255)          gCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGCCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCAGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTTCCCAAGCACGCAAGCCgg                                                                                                                                                                                                                                                                                 >  1:183572/1‑289 (MQ=255)                  gcagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                                                                                                                                                                                                                                                                          >  1:590129/1‑288 (MQ=255)                  gcagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGACCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCaa                                                                                                                                                                                                                                                                         <  1:396772/289‑1 (MQ=255)                   cagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCaa                                                                                                                                                                                                                                                                         <  1:447567/288‑1 (MQ=255)                    agcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACa                                                                                                                                                                                                                                                                       >  1:437863/1‑289 (MQ=255)                        aaCTAACGACGAAGCAACCGTTGATGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt                                                                                                                                                                                                                                                                   <  1:471193/289‑1 (MQ=255)                                   cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACCGGCTAATGCCGCGcc                                                                                                                                                                                                                                                        >  1:267839/1‑289 (MQ=255)                                     gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCCGGAACAGCAACAGGCTAATGCCGAGCCAg                                                                                                                                                                                                                                                      >  1:158182/1‑289 (MQ=255)                                       aaCCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCa                                                                                                                                                                                                                                                       <  1:565977/286‑1 (MQ=255)                                                    ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                                                                                                                                                                                                                        >  1:81536/1‑288 (MQ=255)                                                       cGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCCATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCACAGCCGCCGTCGAAGCCGCTATTGCCCGGGCACAAGCACGCAAGCTGGAACAGCAACAGGTTACTGCCGGACCAGAACCACAAGTcacatcg                                                                                                                                                                                                                                     >  2:461207/1‑282 (MQ=17)                                                                gCCGTGCAAGCAGATATCGCCCGCGCCCAAGCGGGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAGCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccg                                                                                                                                                                                                                           <  1:279127/289‑1 (MQ=255)                                                                 ccGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCCAGCAAGCTGGAACAGCAACAGGCTAATGACTGACCACAAAAACAAGTCGAAcctgccaaagact                                                                                                                                                                                                                           >  2:303069/1‑277 (MQ=14)                                                                  cGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCCAAGCACGCAAACTGGAACAGCAACAGGCTAATGCCGAGCCAGAAAAACAGGTCGATCCGCGCAAAgccgcc                                                                                                                                                                                                                         >  1:301921/1‑289 (MQ=255)                                                                                       cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACCGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCCCGCAAGCCGGAACAGCAACAGGCTAATGCCGAACCAGAACAACAGGTAGATCCGCGCAACGCCGCCGCAGAAGCCGCCATTGCCCGa                                                                                                                                                                                                    >  1:223558/1‑289 (MQ=14)                                                                                       cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACCGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCACAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                                                                                                                                                    >  1:243900/1‑289 (MQ=255)                                                                                       cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                                                                                                                                                    >  1:586923/1‑289 (MQ=255)                                                                                            aaGCGCGCAAGCTGGAGAAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                                                                                                                                                               <  1:483503/289‑1 (MQ=255)                                                                                                                 aaCAGGGTAATGCGGAACAAGAAAAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCa                                                                                                                                                                         <  1:461207/289‑1 (MQ=255)                                                                                                                          aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCTGGCCAGCAAACCCAGTCGCACCGCGCAAAGCCGCCGCCCGACCCGCTATTGGCAGTACCAAGGCGCGGCAAAC‑GAAACCGCCACCGGATaa                                                                                                                                                                >  2:627607/1‑289 (MQ=14)                                                                                                                          aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCAAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCACAAGCCA‑GCAAACGGGCACAGCAACCGGCtca                                                                                                                                                                >  1:670582/1‑287 (MQ=255)                                                                                                                           tGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGCACACGTCGATCCGCGCCAAACCGCCGCCGAAGCCGCTATTGCCTGTGCCAAAGCGCGCACACTGGAACAACAACAGGCTAATGCCGAGCCCGACCAACAGGCCGATCCCCGCAAAGCCGCGCCCGCAGCCGCTATTGCCCGAGCCAACGCGC‑GCAAACTGGacctccacaccgcc                                                                                                                                                                  >  2:451884/1‑273 (MQ=14)                                                                                                                              ggAACCAGAACAACAGGTCGATCCGCGCAAAACCGCCGTCGAAGCAGCTAGTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCg                                                                                                                                                            <  1:397360/289‑1 (MQ=255)                                                                                                                              gcgaCCAGAACAGCAGGTCGAGCCGCGCAGAGCCGGCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGc                                                                                                                                                             <  1:451884/285‑1 (MQ=255)                                                                                                                                aaCCAGACAAACAGATAGACCCGCGCAAAGCCGCCGTCGGAGCCGCTATGACCCGTGCCACAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAGAGGTCGATCCGCGTAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGa                                                                                                                                                          <  1:678275/289‑1 (MQ=17)                                                                                                                                         aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCaa                                                                                                                                                                        <  1:67135/266‑1 (MQ=255)                                                                                                                                            gaccagaacgacggaaaGCCGCCGTTGAGGCCGCGATGGCACATGTGTAACAGCGCAAGCGGGAACAGCGACAGGCTAATGCGGAACCAGAAGACCGGGCCGATCCACGTAAAGCCGCCGTTGAAGCGGCTATGGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                                                                 <  2:590129/272‑1 (MQ=14)                                                                                                                                                                 ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCCAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCAAATGCCGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCAATTACCCGGAGCAAAGCAC‑GCAAACGGGCACAGCAACCGGCCAACGCAGAACCAGAAGAACAGGTTGACCCGCGCAAAggcg                                                                                                                          >  2:268810/1‑285 (MQ=14)                                                                                                                                                                         gcagCTATTGCCGGTGCCAAAGCGCGCAAGCAGGAACAGCAACAGGATAATACGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACGGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAAAGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGa                                                                                                                   <  1:57917/284‑1 (MQ=255)                                                                                                                                                                                    cacggCCAAACAGCCCAAGCCGGAACAACAACAGGCAAATGCGGAACCAGAAGAACAGGTAGATCCACGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCa                                                                                                      <  1:583353/285‑1 (MQ=17)                                                                                                                                                                                                  caaaGCAGGAACAGCAGCAGGCTACTGCGGAACCAGAAGACCCGGCCGATCCGCGCAAAGCCGCCGTCGAAGCTGCTATTGCCCGTGCAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACg                                                                                       <  1:31281/287‑1 (MQ=255)                                                                                                                                                                                                              aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGAAAGCTGGAACa                                                                           >  1:136581/1‑289 (MQ=255)                                                                                                                                                                                                                  acgggCTAATGGGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCaa                                                                       <  1:303069/286‑1 (MQ=255)                                                                                                                                                                                                                  aCAGGCTAATGCGGAACCAGAAGAACAGGACGATCCGCGCCAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCaa                                                                       <  1:627607/289‑1 (MQ=255)                                                                                                                                                                                                                          aTGCGGAACCGGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTCGCAGCCGCTATTGCCCGAGCCAAAGCCC‑GCAAACGGGAACAGCAACCGGCGTATGCGGGACCCGAAAAACAGGGTGATCCGCGACAAGCTGCCGGCGAAGGGGGCTTTGCCCGGCCCCAAG‑CAACAAAACTGGAACAGCAAACGGCtga                                                               >  2:70488/1‑287 (MQ=14)                                                                                                                                                                                                                           tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAAt                                                              >  1:68060/1‑289 (MQ=255)                                                                                                                                                                                                                                aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCTAAGCCGCTATTGCCCGCGCCAAAGCGC‑GCAAACGGGAACAGCAACAGGCTAATGAGGAGCCAGAAGAGAAGGTTGATCCGCGCACAGCTGCCCTCGAAGCGGCTATTGCACGCGCCAAAG‑CACCCATGCGCGCACAGCAAAAGGATacagagca                                                         >  2:188754/1‑282 (MQ=17)                                                                                                                                                                                                                                aaCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCCATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCAC‑GCAAACGGGAACAGCAACCGGCTAATGCCGATCCAGAAGAACAGGTTGATCCGCGCAAAGCCGCCGTCGAAGCCGCTCTTGCCCGCGCCCAAG‑AACCGAAGCTGGCAcaacaaacgggtaaagccgt                                                         >  2:357571/1‑275 (MQ=11)                                                                                                                                                                                                                                aaCCAATAGACCAAGTCGATCCGCGCAAACTCGGCGTCAAAGCGGCTATTGCCCGTGCCAAAGAACCCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGt                                                         <  1:117618/289‑1 (MQ=17)                                                                                                                                                                                                                                 aCCAGAAAAACAGGTCGAGGCGCGCAAAGCAGCCGTCGAAGCGGCTAATGCCCGTGCCAAAGCACGCGAGCAGGAACAGCAACAGGCTAATGCTGAGCCAGAACAACAGGTCGAGCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGt                                                         <  2:68060/288‑1 (MQ=17)                                                                                                                                                                                                                                  ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAc                                                       <  1:268810/289‑1 (MQ=255)                                                                                                                                                                                                                                   cagcagAACAGGTCAATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCAAAAGCAGGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAcc                                                      <  1:358657/289‑1 (MQ=255)                                                                                                                                                                                                                                             ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                            >  1:290964/1‑289 (MQ=255)                                                                                                                                                                                                                                              gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAgg                                           >  1:393654/1‑289 (MQ=255)                                                                                                                                                                                                                                                cGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAAAACAGGtt                                         >  1:426002/1‑289 (MQ=255)                                                                                                                                                                                                                                                   tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTCATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCCCGCCCCAAAGCCACG‑AAGCTGGAACAGCAAGAGGATAAATAGGTACCCGAAGAACAGGTTgc                                       >  2:597420/1‑287 (MQ=17)                                                                                                                                                                                                                                                   tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCC‑GCAAACGGGAACAGCAACAGGCTGATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCCGCCGTCGCAGCGGCTATTGCACGCGCTAAGG‑CACACAAACCGGAACAGCAAAAGGCTAATACGGTGCCagaagacgcagtcgat                                      >  2:489915/1‑279 (MQ=14)                                                                                                                                                                                                                                                   acccggcAAAGCGGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAGTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGAt                                      <  1:357571/284‑1 (MQ=255)                                                                                                                                                                                                                                                          aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCCGTCGCAGCGGCTATTGCACGCGCCAAAG‑CACGAAAGACAGAACAGCAACAGGCTAATGCGGTCCCTTAAGCACAGCttacgctgccca                               >  2:197997/1‑279 (MQ=17)                                                                                                                                                                                                                                                                                ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGATCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTCCCCGAGCCAAAGCGC‑GCAAGCGTGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAt         <  1:159834/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                ttGCCCGGGCCAAAGCACGCAAGCGGGAACAGAAAAAGGCTAATGCCGAGCCAGAAGAACAGGATGAGCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGACCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAt         <  1:364510/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                   cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCT‑ATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAAAG‑CACGCAAACTGGAACAACAACAGGCTAATCCGGTACCAGAAGAACAGGATGAACCGCGCAAAGCGGCAGTTTACTCGGCTATTGc      >  2:624132/1‑288 (MQ=18)                                                                                                                                                                                                                                                                                       tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                                                                                                                                                  >  1:60908/1‑83 (MQ=25)                                                                                                                                                                                                                                                                                       tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                                                                                                                                                  <  2:60908/83‑1 (MQ=25)                                                                                                                                                                                                                                                                                       tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCACACGGCTAATACGGAGCCAGAAGAACAGGCTGAACCGCGCAAAGCAGCCGTCGAAGCGGCTATTGACCGGTCCAACG‑CACGCAAGACGGACCAGCAACAGGCTCACGAGGTCCCAACAGACCACGTTGAACGCAGCAACGGTGCAGGTGACGCGGccgatcaccc   >  2:451154/1‑279 (MQ=14)                                                                                                                                                                                                                                                                                       tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgc  >  1:407239/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                       tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCCATTGCCcgc  >  1:581582/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                       agCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGC‑GCAAACTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCCAAAG‑CCAGAAAGCTGGAACAGCAACAGGCTAAAGCCGAGGCAGAAAAACAGGTTCATCCGCACAAAGCGGCCGTCGCagcagaagttgcccaa  >  1:389375/2‑276 (MQ=11)                                                                                                                                                                                                                                                                                        gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCTGGCAAAGCGGCAGGTGCCGCGGCTACTGCAcgc  >  1:410226/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                          GCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAG‑CACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGC  >  NC_000913/1708527‑1709093

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