Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,842,367	A→C	30.8%	V12G (GTG→GGG)	yicN ←	DUF1198 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,842,367	0	A	C	30.8%	61.4 / 6.7	51	V12G (GTG→GGG)	yicN	DUF1198 family protein
Reads supporting (aligned to +/- strand):  ref base A (18/15);  new base C (2/16);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 2.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGT  >  NC_000913/3842142‑3842585                                                                                                                                                                                                                                    |                                                                                                                                                                                                                            cTACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGAACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAACCGCTTTTCTGGCCCCGGATGTCGGCACCCCAAAACCAACCACAAacacc                                                                                                                                                                                                                           >  1:19653/1‑227 (MQ=255)  tACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTTATCGACCATCGAATCCCC‑GGGTACAACATCGATGTTCGACCGATCGCTGCGACGGCGAATCGCTTTTCTGGCCCCGGATGTCCGCACCCACAAACCACCCACAAacaaaa                                                                                                                                                                                                                          >  1:130150/1‑225 (MQ=255)       tCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGCTGTTCAGCCTATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGTAGTCACCACCCGAAAAACAACcacacacaaacccccc                                                                                                                                                                                                                     >  1:404223/1‑220 (MQ=255)       tCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGAACGACATCGATGTTCAGCCGATCGCTGAGCGGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAACCCTACCACAAACACTACCTcc                                                                                                                                                                                                                     >  2:645375/1‑228 (MQ=255)       tCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAAAGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGAACGACATCGATGTTCAGCCGATCGCTGCGACGGCGAATCGCTTTTCTGGCCCCGGACGTCAGCACCCGAAAAACAACCACAAACACCACCGcc                                                                                                                                                                                                                     >  1:402797/1‑228 (MQ=255)              tgacaatTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTg                                                                                                                                                                                                               >  2:494833/1‑227 (MQ=255)              tgacaatTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGAGCTCAGCACCCGAAAACCAACCACAACCACCACcccccacgcgc                                                                                                                                                                                                              >  1:420226/1‑221 (MQ=255)               tccaatTTGCCAAATG‑GCAACCCCCTGCGGGGGTTTTTGCGGGTGGGCTCCCACCGGACGATGAAAATAACGTAAAAATCCGTCACCAGCTGATTTCCCCATTGGAGCGACCCTCGACCCCAC‑GGGTACGCCATCGATGTTCGGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCAAAACCCAACCACAACCCCTACCGCCAGCGTGGc                                                                                                                                                                                                             <  1:284827/226‑1 (MQ=255)                   ttttaaaggacaggcgaccacCTGCGCGGCTTTTGGCGGGTGCGACTCACCCGAACGATGCAAATAAGCTAAAAATTCGTCACCAGCTGACTTTCCCTTTTGATCGCCCACCGACCCCCC‑GGGTACGACACCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAg                                                                                                                                                                                                          <  1:12106/210‑1 (MQ=255)                         cAGATGAGCAA‑CACCTGCGTGGGGTTTTGCAGGGGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATc                                                                                                                                                                                                   <  2:391647/228‑1 (MQ=255)                            aTGAGCAA‑CACCTGCGCGGCTTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAg                                                                                                                                                                                                <  2:386398/228‑1 (MQ=255)                                 cAA‑CCCCTGCGCGGCTTTTTGCAGGTGCAACTCATCCGGCCGATGCAAAAAACGAAAAAATTCGTCCCCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCGGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCAt                                                                                                                                                                                           <  2:456010/228‑1 (MQ=255)                                  aa‑CACCTGACGGGCTTTTTTCAGGTGCGACTCATCCGGACGATGCAAAAACCGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGACCGACCATCGACTCCAC‑GGGTACGCCATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATg                                                                                                                                                                                          <  2:219789/228‑1 (MQ=255)                                       ccggggcggcgTTTTGCAGGTTCAGCTCACCCGTGCGATGCAAATAACGTAAAAATCTGTCCCCAGCCGACTTTCCCATTTGTTCTCCCATCGACTCCAC‑GGTGACGACATCGATGTTCAGCCGATCGCGGAGACGCCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGtt                                                                                                                                                                                      <  1:30279/223‑1 (MQ=255)                                           cgcgGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCCACCGAAAACCAACCATAAACACAACCGCCAGCGGGCGGGGCATTAACCAGATCATGCATTCTcc                                                                                                                                                                                  >  1:521399/1‑228 (MQ=255)                                           cgcgGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGa                                                                                                                                                                                                           >  2:12106/1‑203 (MQ=255)                                            gcggcgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTCGCCCCGGATGTCAGCACCCGACAACCAACCACAAACACTCCCGCCAGCGTGGCGAGCATTTTCCAGCTCATGACtttccct                                                                                                                                                                                 >  2:163099/1‑223 (MQ=255)                                                cgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGCATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCCACCACCAAACCTACCGCCAGCGGCGCGAGCATTATCACGATCATGGGTTCCCCTGCgg                                                                                                                                                                             >  1:345647/1‑227 (MQ=255)                                                cgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGt                                                                                                                                                                             >  2:177478/1‑228 (MQ=255)                                                   tttGCGGGTGCGCCTCTCCCGGATGTTGCAAATAAAGTAAAAATTCGTCACCCGCTGACTTTCCCATTTGACCGACCATCGACTCCGT‑GGGTACGAACTCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAACCCAACCACAACCCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAAt                                                                                                                                                                          <  1:177478/228‑1 (MQ=255)                                                     ggaagggtcgccTCAACCGGAACAGGCAAATAACGTAAAAATTGGGCACCAGTTGGCTTTCCCATTTGACCGCCCATCGCTCCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCa                                                                                                                                                                        <  2:586425/217‑1 (MQ=255)                                                        aGGTGCGACTCACCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTCCCCATTTGATCGACCATCGACGCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGaa                                                                                                                                                                     <  2:463382/228‑1 (MQ=255)                                                                    tCCGGACGATGCCAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGCCTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGAACCCGAAAAACAACCAAAAACCATACCGCCAGCGTGGCGAGCATTATCACGATCATGAATTCACCTCCGTAATCACAACtactagcatc                                                                                                                                                          >  1:564362/1‑218 (MQ=255)                                                                             ggaaaaTAACGTAAAAATTCGTCACCAGCTGTCTTTCCAATTTGAGCGACCAGCGACTCCAG‑GGGTACGCAATCGATGTTCGCCAGATCGCGGAGCCGGGAAATCGCTTTCCGGGCCCCGGAGGTCAGCACCCAAAAACCAACCCAAACCCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCt                                                                                                                                                <  1:645375/225‑1 (MQ=255)                                                                                aaaCAACGTAAAAATTCGTCACCAGCTGCTTTTCCCATTTGAGCGACCATCGACTCCCG‑GGTTACCACATCGATGTGCAGCCGATCGCTGAGACGGCGAATCGCTGTTCTGGCCCGGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGc                                                                                                                                             <  2:586913/228‑1 (MQ=255)                                                                                       gTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCAAAACCATAAGTTCTCCTGCGGAATCAGAACTCAAGATATTGTGAAACCCCCGCCaagtac                                                                                                                                      >  1:589166/1‑223 (MQ=255)                                                                                        tAAAAAGTCGTCCCCAGCTGATTTTCCCATTTGATCGACCATCGACCCCAC‑GGGTACGACATCGATGTTCAGCCGACCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCACCCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTAc                                                                                                                                      <  2:19653/227‑1 (MQ=255)                                                                                              ttggtcccccccttattttCCCTTTTGATCGACCATCG‑CCCCACGGGGTCCGCAATCGAGTGTCCGCGGATCCCTGAGACGGGCAATCGCTTTTCTGGCCCCGGATGTCAGCCCCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATcca                                                                                                                               <  1:494833/212‑1 (MQ=255)                                                                                                cGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGCCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTATTccac                                                                                                                              >  2:492123/1‑227 (MQ=255)                                                                                                   tccccgCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGa                                                                                                                          <  2:125038/223‑1 (MQ=255)                                                                                                            gcTTTCCCATTTGATCGCCCTTCGACCCCAC‑GGGTACGACCTCGATGGTCCGCCGATCGTTGAGTCGGCGAATCGGTTTTCTGGCCCCGGGTGCCAGCACCCGAAAACCACCCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTc                                                                                                                 <  1:7343/227‑1 (MQ=255)                                                                                                                tCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCa                                                                                                                >  2:355055/1‑225 (MQ=255)                                                                                                                      ttGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACa                                                                                                       >  2:11110/1‑228 (MQ=255)                                                                                                                          tCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACagag                                                                                                    >  2:643612/1‑227 (MQ=255)                                                                                                                                agcggcTCCAC‑GGGTACGACATGGATGTCCAGCCGATCAGTGAAACGGCGAATCGCTTTTCTGGCCCCGGATGTACGCACCCAAAAACCAACCACAATCACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGac                                                                                             <  2:556242/223‑1 (MQ=255)                                                                                                                                   caCTCCAG‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCATAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccac                                                                                          <  2:521399/227‑1 (MQ=255)                                                                                                                                    aCTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCCGAATAATCTTCAGTGCAAACAGACCTGCGaccacc                                                                                         >  2:50464/1‑228 (MQ=255)                                                                                                                                             ggtaacgactCGCTGTCCCGCCGATCGCTGAGGCGGGGAATCGTTTTGCTGGCCCGGGATGTCAGCCCCCGAAAACCAGCCACAAACCCTCCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCCCCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacac‑a                                                                                <  1:435033/219‑1 (MQ=255)                                                                                                                                             ggatcGACATCGATGTTCAGCCGTTCGCTGAACCGGCTAATCGCTTTTCTGGCCCCGGATGGCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacac‑a                                                                                <  2:589166/224‑1 (MQ=255)                                                                                                                                                aCGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCCCcacaacaa                                                                               >  1:528922/1‑223 (MQ=255)                                                                                                                                                                    aTCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTACTCAGAACTCATGATATCTGGAAACCCTCGCCAGTTCATAATCCACCAGAATAATCTTAATGGCAAACCGAGCTGCGACCCCCACC‑CACAACAGGTTCaaaaccggccagaccc                                                         >  1:530062/1‑213 (MQ=255)                                                                                                                                                                     tcccTGAGACGGCATATCGGTTTTCCTGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGcc                                                         <  2:404223/224‑1 (MQ=255)                                                                                                                                                                     ccGCTGAGCCGGCAAATCGCTTTTCTCGCCCCGGATGTCCGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGcc                                                         <  2:402797/226‑1 (MQ=255)                                                                                                                                                                           agccGGCGAATCGCTTTTCGGGCCCCGGATTTCAGCACCCGAAACCCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                  <  2:653044/225‑1 (MQ=255)                                                                                                                                                                                   acccggttttCGGGCCCCGGTTGGCAGCCCCCAAAACCCAAACAAAAACCCCACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAACCCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCa                                          <  1:50464/222‑1 (MQ=255)                                                                                                                                                                                         ttttCTGCCACCGGATGTGACCCCCCGAAAACCAACCCCAAACCATACCGCAAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCatgatg                                     <  2:326518/227‑1 (MQ=255)                                                                                                                                                                                         ttttCGGGCCCCGGATGTCAGCACCCGAAAACAAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAGCCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCatgatg                                     <  2:432490/227‑1 (MQ=255)                                                                                                                                                                                             cgggCCCCGGGAGGTAGCCCACGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGAATTAACCAGATCATGAGTTCTCCTGCATAATCAGAACTCATGATAGCTGGAAACCCGCGCCAGTTTCTAAGCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCTCCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTa                                <  1:643198/226‑1 (MQ=255)                                                                                                                                                                                             cTGGCCCCGGATGTCAGCACCCGAAAACCAACCCCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTa                                <  2:420226/228‑1 (MQ=255)                                                                                                                                                                                              aggCCCGGGATGTCAGCACCCGAAAACCAACCCCAAACCCCACCCCCAGCGTGGCGAGCATTAACCAAATCACTAGTTCTCCTGCGTTATCAGAACCCATTATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCCTCAGTGCAAACAGAACTGGGACCACCACCACAC‑ACAGGTTCAAATCGCGCCAGCGCCCGGTGCATGCTTTCATGATGCAGTa                                <  2:530062/226‑1 (MQ=255)                                                                                                                                                                                                         ttcaaccaccGAAAACCACCCAAAAACCCTCCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTTCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCTCAC‑TCAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGcc                    <  2:102445/220‑1 (MQ=255)                                                                                                                                                                                                               cccaaaaaaCCACCCCAAAACCCTCCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc              <  1:11110/223‑1 (MQ=255)                                                                                                                                                                                                                 ccGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCaa            >  2:29893/1‑228 (MQ=255)                                                                                                                                                                                                                           accaCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATCATTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGt  >  2:356091/1‑228 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                            CTACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCAC‑GGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACAC‑ACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGT  >  NC_000913/3842142‑3842585

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