Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,156,549	T→C	21.7%	S180P (TCC→CCC)	argB →	acetylglutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,156,549	0	T	C	21.7%	136.6 / 3.4	69	S180P (TCC→CCC)	argB	acetylglutamate kinase
Reads supporting (aligned to +/- strand):  ref base T (37/13);  new base C (2/15);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 8.59e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGACGGCGACAGCGTCAAAGTGACCCAGCTTGATGAAGAGTTAGG‑TCATGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAG  >  NC_000913/4156323‑4156783                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            cGGTGACGGCGACAGCGTCAAAGTGACCCAGCTTGATGAAGAGTTAGG‑TCATGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCg                                                                                                                                                                                                                                 >  2:851688/1‑239 (MQ=255)      aCGGCGACAGCGTCAAAGTGACCCAGCTTGATGAAGAGTTAGG‑TCATGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGCCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCCCC‑GGCGGCAACGCTGGGCACGGCACTGCATTTGCCCCccgccgctcccgccat                                                                                                                                                                                                                             >  1:504799/1‑228 (MQ=255)          cGACAGCGTCAAAGTGACCCAGCTTGATGAAGAGTTAGG‑TCATGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGAGCTGATTTTGCTCTCCGCCGTAAGCGGCATtctc                                                                                                                                                                                                                         >  2:117708/1‑238 (MQ=255)                 ggcAAAGTGACCCAGCTTGAGGAAGAGTTAGGTTCAT‑TTGGACTGGCGCCGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCa                                                                                                                                                                                                                  <  2:551179/236‑1 (MQ=255)                          aCCCAGCTTGATGAAGAGTTAGG‑TCATGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGACCGGCTTTTGCCCTCCCCCGCCGGCGCCATTCCCGCCGGCCAGAGGCAAca                                                                                                                                                                                                        >  1:570039/1‑238 (MQ=255)                           ccccGCTTGGAGGAGGGTTGGG‑TCTGGTTGGACTGGCGGAGCCGGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATTGGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACg                                                                                                                                                                                                        <  1:434078/238‑1 (MQ=255)                               tttttgtgaAGGTTTGGG‑TCATTTTGGGCTGGCGCAGCCAGGTTCGCCTAAGCTTATAACCTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTg                                                                                                                                                                                                   <  1:884362/233‑1 (MQ=255)                                               gg‑TCATGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGCCTCGGATCTGATTTTGCTCTCCGACGTCAGCGGAATTCCCGACGGCAAAGGGCAACGCATTGCCGAAATGAAcgccgc                                                                                                                                                                                   >  2:732949/1‑239 (MQ=255)                                                     tGTTGGACTTGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCCGGCGGCAACGGCGCT‑GGCGGAAACGCTGAGCGCGGCGCGGATTTTGCTCTCCGAAGGCCGCGGCATTCTCGCAGGCAAAGGGCAACGCATGCCCGCCATGACCGCCGGCAAAg                                                                                                                                                                              >  1:475633/1‑239 (MQ=255)                                                     tGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCGGTATTTGCCTCCCGGCCCTAGCGGGATCCCCGACGGCAAAGGGCACCGCATTCCCGaaagaaaacccccgaaag                                                                                                                                                                              >  1:748280/1‑224 (MQ=255)                                                     tGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGAGCTGATTTTGCTCTCCGCCGCCCGCGGCATTCCCGCCGGCAAAGGGCAAAGCATTGCCGAAAAAAACGGCGCGaaaa                                                                                                                                                                              >  1:366142/1‑238 (MQ=255)                                                     tGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGAGGATCTGATTTTGCTATCCGCCTTACGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGAcacacgcgcag                                                                                                                                                                               >  2:501314/1‑228 (MQ=255)                                                       ttGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCCCCGCGAAAGcc                                                                                                                                                                            >  2:877887/1‑238 (MQ=255)                                                           aCTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCCCGAAAGCAGaac                                                                                                                                                                        >  2:306149/1‑239 (MQ=255)                                                                 ggcgccgtgttggcttagctttATCAACTCCTGGCGGGGGAACGGTTGTCCGGCGGTGGTAAGGTCAAGTGGG‑AGAACAAGCGAAGGGCAACGAAT‑GAACGTCAGGGCCCGCCGGGGGGAAGCGGCGCT‑GGGGGAAACGCTGGGCGCGTATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGCCGGCAAAGGGAAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGa                                                                                                                                                                  <  1:194447/219‑1 (MQ=255)                                                                        ggTTCGCCTAAGCTTTTCAATTCTTTGGTGGAAAACGGTTACTGGCCGGGGGTCGGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCa                                                                                                                                                           <  2:366142/239‑1 (MQ=255)                                                                         gTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAg                                                                                                                                                          >  2:808912/1‑239 (MQ=255)                                                                          ttCGCCTAAGTTTATCACCCTTTGGTTGGAGAGCGGTTATCTGGGGGTGGTCAGCTCCATTGGC‑GTAACAGACGAGGGGCAACTGAT‑TAACGTCAATGCCGAGCGGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCCCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAg                                                                                                                                                          <  1:386710/238‑1 (MQ=255)                                                                            cgctcaaGCTTATTACATCCTTGGCGGAGAGCGTTTATCTGCCGTGGGGCAGCCCCATGGGC‑GTAACAGACGAAGGGAAACTGAG‑GAACTCCAAGGCCGACCAGGCGGCAACGGCGGT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCCCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATGGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAggg                                                                                                                                                        <  1:815255/233‑1 (MQ=255)                                                                              ccTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGCCGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAAAACAACTGATTTAGCAGGGCat                                                                                                                                                     >  2:756678/1‑239 (MQ=255)                                                                               cTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCACGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGaacacatatatagccagggcaaa                                                                                                                                                    >  1:261338/1‑220 (MQ=255)                                                                                       ttccacctcTTGTTGGAGAAGGGTTTGCTTCGGGTGGT‑AGACCCAATGGCGGAAACAGACAAAGGGCAATTGAT‑GACCGTCAATGCCGACCAGGGGGCAACGG‑GCTGGGCGGAACCGCTGGGCGCGGATCTGATTTTGCTCCCCGACGTCAGCGGCATTCTCGACGGCAAAGGGAAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGa                                                                                                                                            <  1:851688/231‑1 (MQ=255)                                                                                       tacaccTCCTTGCGGGAGAACGGTTGTGGGCCGGGTGCCAGTCCCATGGGG‑GAAACAGGCGAGGGGAAAGTGAT‑GAACGTCAATGCCCACCAGGCGGCAACGGGGGT‑GGGGGCAACGCTGGGCGCGGATCGGATTTTGCTCCCCGACGTCAGCGGCATTCTCGACGGAAAAGGGCAACGCATGGCCAAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGa                                                                                                                                            <  1:222890/234‑1 (MQ=255)                                                                                             tCCTTGCTGGAAACCGTTTTTGTGCGGGTGGTCAGCTCCTTTGGC‑CTAACAGACGAAGGGAAACTGAG‑GACCGTAAATGCCGCCCGGGGGGAAACGGGGCT‑GGGGGCAACGGTGGGGGCGGATCTGATTTTGCCCCCCGAGGTCGGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCAt                                                                                                                                      <  1:363178/239‑1 (MQ=255)                                                                                             ctctggCTGGGGAGCGGTTTTCCGCCGGTGGTCCGCCCCAGTGGC‑GTAACAGACGAAGGGCACCGTAT‑GAGCGTCAAGGCCGACCAGGCGGCACCGGCGGT‑GGCGGAAACGCTGGGCGCGGATTTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGAAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCAt                                                                                                                                      <  1:839417/234‑1 (MQ=255)                                                                                              ccTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGCCGTCAGCGGCATTCTCGCCGGCCAAGGGCACCGCATTGCCGAAATTACCGCCGCGAAAGCAAAACAAATGATTGAGCAGGGCATTATTACTGACGGCaag                                                                                                                                     >  1:612243/1‑237 (MQ=255)                                                                                                 tGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGCTTGAGCAGGGCATTATTACTGACGGAATGATa                                                                                                                                  >  2:363178/1‑239 (MQ=255)                                                                                                 tGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACggacggata                                                                                                                                  >  2:160151/1‑232 (MQ=255)                                                                                                 tGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGCCCTGAtc                                                                                                                                  >  2:247056/1‑238 (MQ=255)                                                                                                 tGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACGGATTGAGCAGGGCATTATTACTGACGGCATGATa                                                                                                                                  >  2:542432/1‑239 (MQ=255)                                                                                                  tggggggaACGTTTATCTTCGGGGGGTCAGCTCCATTGGG‑TGAACGGCCGAAGGGCAACTGAG‑GAACGCCAAGGCCCACCGGGCGGCACCGGCGCG‑GGCGGCAACGCGGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATag                                                                                                                                 <  1:503158/233‑1 (MQ=255)                                                                                                      gagaACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCCCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAAACCATTGCCGAAATGACCGTCCCGGAAGCGGAACAACTGAATGGGCAGGGCATTATTACCGAGGGCAAAATagggaa                                                                                                                             >  1:756041/1‑239 (MQ=255)                                                                                                              ttctCTGGCGGTGGGCAGCCTCATTGGC‑GTAACAGACGGAGGGCAACTGAA‑GAACGGCATGTGCGACCAGGGGGCAAGGGCGCT‑GGCGGCAACGCGGGGCGCGGATCTGATTTTGCTCGCCGACGTCGGCGGCATTCTGGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAAcg                                                                                                                     <  1:230959/236‑1 (MQ=255)                                                                                                               tATCTGCCGGTGGGCAGCTCCATGGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAAGGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAAcg                                                                                                                     <  2:710046/238‑1 (MQ=255)                                                                                                                   tGCCGGTGGTTAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAAAAACTGATTGAGCAGGGCATTATTACTGACGGCATGAAagggaaagagaaaaggggg                                                                                                                >  1:235382/1‑232 (MQ=255)                                                                                                                   tGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGGAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCCTTCTCGACGGCAAAGGGCAAGGCATTGCCGAAATGACCGCGGCGAAAGCAAAAACACTGGGTGGGCAGGGACTTTTTACGGAGGGCAAGAAagagaaagggaaaggaggg                                                                                                                >  1:103617/1‑232 (MQ=255)                                                                                                                        gtggtCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGGGAACGCGG‑CGCGGGa                                                                                                          >  2:50633/1‑239 (MQ=255)                                                                                                                        gtggtCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGCCAACGCTGGGCGCGGACCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAAAGCATTGTCGAAATGACCGCCGCGACAGCAGCACAAGCGATTGACGAGGGACTTATTATTAACGGCTGAATTGTGAACGAGAAAGGGG‑CGTTGGa                                                                                                          >  1:733382/1‑239 (MQ=255)                                                                                                                         ggttccgccccATTGGC‑GCAACAGGAGGAGGGACACTGAG‑TAGCCGCGATCCCGGCCAGTGGGCAACGGGGAG‑GGCGGCAACGCTGGGCGGGGATTTGATTTGTCTCTCCGACTACGCGGGAAATCTAGAAGGGAAAGGGCAACGCATCGCAGAAATGGCCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTCCTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGAt                                                                                                         <  1:448312/230‑1 (MQ=255)                                                                                                                           gtCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCCTTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAAGGAACGGGG‑GGCgggcggg                                                                                                       >  1:622654/1‑235 (MQ=255)                                                                                                                                cTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACGTACGGCATGATAAGGAAAGTGAACGCGG‑CGCTGGATGCGGcc                                                                                                   >  2:609232/1‑238 (MQ=255)                                                                                                                                  cttcagtggggtaacgggcaggggCAAGTGGTGGAACGTAAAGTCCGACCGGGCGGCAAGGGGGCG‑GGGCGCAACGCGGGGCGCGGATTTGATTTTGCTCACCGACGCCAGCGGCATTCCCGACGGAAACGGGAAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGccc                                                                                                  <  1:864618/218‑1 (MQ=255)                                                                                                                                   agaggggtcgaaacggaagaGGGCAACGGAT‑ATACGTCAAGGCGGTCCGGGCGGCAGCGGCGGT‑GGGGGCAACCGGGGGCGCGGCTCTGATTTTGCCCCCCGGCGTCGGCGGCATTCGCGACGGCAAAGGGAAACGCATTCCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGccc                                                                                                  <  1:194432/218‑1 (MQ=255)                                                                                                                                     ttGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGGAGAACAACTAGTTTGGCCGGGCATTTTTATTGACGGCTGGATAGGGAAAGTGAAAgcgg‑cgacggtggcgccgcgcccc                                                                                             >  1:534883/1‑221 (MQ=255)                                                                                                                                           gTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGGGGCCCGCACGCggg                                                                                         >  2:685048/1‑238 (MQ=255)                                                                                                                                              aCAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCCTGATAGTGAAAGTGAACCGGG‑CGCGGGATGCGGCCCGCACGCTGGGcgg                                                                                     >  1:143192/1‑237 (MQ=255)                                                                                                                                                   cGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTggg                                                                                        >  2:194432/1‑231 (MQ=255)                                                                                                                                                   cGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCgg                                                                                >  2:266345/1‑239 (MQ=255)                                                                                                                                                                           aTGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGTTCTGATTTTGCCCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACGGCCGCGAAAGCAGAACAACTGATTGGGCAGGGCATTATTATTGACGCACTGATACGAGAAGGGAACCCGGCCGCGGTATG‑GGCCCCCACGCTGGGCCGTCCGGgagaaacccccccccggccgca                                                          >  1:405440/1‑220 (MQ=255)                                                                                                                                                                            tgccccccgggccgcCACGGCGGT‑GGGGGCAACCCTGGGCGGCGTTTTGTTTTTGCCTCCCGACGTCGGCGGCATTCCCGACGGCAAAGGGCAGCGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATg                                                        <  1:266345/230‑1 (MQ=255)                                                                                                                                                                                 aCCAGGGCGAAACGGCGCG‑GGGGTCCACGCTGGGCCGGGATCTGATTTTGCCCTCCGACGTCAGCGGCATTCTCAACGGCAAAGGGCAACGCAGTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGag                                                   <  2:334335/239‑1 (MQ=255)                                                                                                                                                                                 aCCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCGGACTGCGGCCCGCACCCTGGGCCGTCCGGGAGATATCGCCTCCGGGccccatcgggac                                                   >  1:700463/1‑228 (MQ=255)                                                                                                                                                                                            agggggcT‑GGGGGGAAGGGTGGGCGCGGATCTTATTTTGCTCCCCGAAGTCAGCGGCATTCCCGACGGCAAAGGACAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCgg                                         <  2:570039/233‑1 (MQ=255)                                                                                                                                                                                             gggCGCT‑GGGGGAAACGGGGGGCGCGGATTTGATTTTGCCCCCCGAGGTCAGGGGCATTCTCGACGGAAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCa                                       <  2:263302/238‑1 (MQ=255)                                                                                                                                                                                             cGGCGTG‑GGGGGAAACGGTGGGCGCGGTTTTTTTTTTGCCCCCCGACGTCAGCGGAATTCCCGACGGAAAAGGGCACCGCTTTGCCGAAATGCCCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCa                                       <  1:306149/239‑1 (MQ=255)                                                                                                                                                                                             cGGCGCT‑GGGGGCAACGCTGGGCGGGGATCTGATTTTGCTCCCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGCCCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCa                                       <  2:504799/239‑1 (MQ=255)                                                                                                                                                                                             cGGCGCG‑GGGGGCAACGCTGGGCGGGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCa                                       <  2:622654/239‑1 (MQ=255)                                                                                                                                                                                               gcgcT‑GGGGGCAACGCGGGGCGGGGATCTGATTGTGCTCCCCGGCGTCAGCGGCATTACCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCa                                       <  2:157324/237‑1 (MQ=255)                                                                                                                                                                                                   gggggggcAACTCGGGGGGGGGATTTTATTTTCTCCCCCGACGTCAGGGGTATTCTGGACGGCAAAGGGCACCGCATGGCCGAAATGACCGCCGGGAAAGCAGAACAACTGAGTGAGCAGGGCATTATTACTGACGGCATGATAGGGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGtt                                  <  1:37722/234‑1 (MQ=255)                                                                                                                                                                                                    tgggggcacgtcTGGGCGCGGATTTGTTTTTGCTCTCAGAGGTCAGCGGCATTCACGACGGAAAAGGGCAACGCATTGCGGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGttt                                 <  1:232532/228‑1 (MQ=255)                                                                                                                                                                                                       cggcAACGGTGGGCGCGGATCGGATTTTGCTCCCCGACGTCAGGGGAATTCGCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTaa                               <  2:700463/238‑1 (MQ=255)                                                                                                                                                                                                           aaCGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCCCGCTGGGCCGGCCGGTAGATATGGCCTCCTGGGGTCATCCGGGGCCGCTTCCGGCCCTGTTTAACGGTa                          >  1:183342/1‑239 (MQ=255)                                                                                                                                                                                                                tGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACGGTTTAACGGTATGccc                     >  2:284470/1‑238 (MQ=255)                                                                                                                                                                                                                tGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCCGGGCCGTCCGGTAGATAGCGCCTCCTGGCGGCATGCCGGGCAGCTTCCGGCCCTGTTTAACGGTaagacg                     >  1:768399/1‑234 (MQ=255)                                                                                                                                                                                                                  ggCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACTCGG‑CGCTGGATGAGGCCGGCACCCTGGGCCGTCCGGGTGATTTCCCCCCCTGTCGGCCTGCGGACCAGGTTCCggaccgttttaagggaagacgaca                   >  1:517496/1‑217 (MQ=255)                                                                                                                                                                                                                          tttGATTTTGCCCCCCGACGTCAGGGGCATTCACGACGGCAAAGGGCACCGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTAcg            <  2:347319/236‑1 (MQ=255)                                                                                                                                                                                                                               ttttGCGCTCCGACGCCAGCGGCATTCGCGACGGAAAAGGGAAACGCTTTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGAtt      <  2:43358/239‑1 (MQ=255)                                                                                                                                                                                                                                   gCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCGCCTGGCGTCCTGCCGAGCAGCTTCCGGCCCTCTTTAACGGTATCCCGATGGGTACgggggatttgg  >  1:43358/1‑232 (MQ=255)                                                                                                                                                                                                                                   gCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCATCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTa   >  2:112933/1‑238 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            CGGTGACGGCGACAGCGTCAAAGTGACCCAGCTTGATGAAGAGTTAGG‑TCATGTTGGACTGGCGCAGCCAGGTTCGCCTAAGCTTATCAACTCCTTGCTGGAGAACGGTTATCTGCCGGTGGTCAGCTCCATTGGC‑GTAACAGACGAAGGGCAACTGAT‑GAACGTCAATGCCGACCAGGCGGCAACGGCGCT‑GGCGGCAACGCTGGGCGCGGATCTGATTTTGCTCTCCGACGTCAGCGGCATTCTCGACGGCAAAGGGCAACGCATTGCCGAAATGACCGCCGCGAAAGCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGG‑CGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAG  >  NC_000913/4156323‑4156783

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