Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,879,458	A→C	16.7%	V43G (GTG→GGG)	yoaG ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,879,458	0	A	C	16.7%	171.8 / 4.0	76	V43G (GTG→GGG)	yoaG	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (43/18);  new base C (3/12);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 6.88e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCCAAAAATAATCCATGACAGAATGCCCATT‑GTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAAT  >  NC_000913/1879223‑1879696                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                aaGCCCAAAAATAATCCATGACAGAATGCCCATT‑GTTTCCTCACTTATCTGTGTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCa                                                                                                                                                                                                                                              <  1:62132/239‑1 (MQ=255)   tcctaaaaaTATACCATGTAAGAATGCCAATT‑GTTTCCGCAATTATTTGTTTTGCATTAGCGGGTTAGGCGCTGATAAAAAGCATAGCACAAAATCGGGAGGGGAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGa                                                                                                                                                                                                                                            <  2:383382/235‑1 (MQ=255)        aaaaTAATCCATGACAGAATGCCCATTGGTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGGATTAACCAGATCtt                                                                                                                                                                                                                                        >  1:58747/1‑239 (MQ=255)          aataaaCCATGACCTAATGCCCCTT‑TGTTCCTCACTTATCTGTTTTTCATTAGCGGGTTAGGCGCTGATAAAAAGCATAGCTAAACATCGGGAGGGCAAGATTTGTTACGGGCAT‑CACGGAGGTTTTTTTTGCGATGGCTCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTa                                                                                                                                                                                                                                     <  2:843187/239‑1 (MQ=255)              atCCATGACAGAATGCCCATT‑GTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGCACTAAGAACGTACGGAATTCACCAGATCTTTTTTc                                                                                                                                                                                                                                   >  1:1088761/1‑237 (MQ=255)                   tGACAGAATGCCCATT‑GTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCTGGGGGGCAAGACTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATAAACGATCGGTGATAATTACCAACCACTAACATCC‑AGTTCGTTTTCCGGGGCAAAAAAAAGAACGGAATTACACACAAATTTTATCACTTcg                                                                                                                                                                                                                             >  2:319788/1‑237 (MQ=255)                    gACAGAATGCCCATT‑GTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGCTCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGAACGGTATTCCCCAAATCTTTTAtacattcag                                                                                                                                                                                                                            >  1:1080588/1‑230 (MQ=255)                    gACAGAATGCCCATT‑GTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAAAACGTACGGAATTCACCAAACCTTTTATCACttaac                                                                                                                                                                                                                            >  2:1337732/1‑235 (MQ=255)                                aTT‑GTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATAAACGATCAGTGATATTTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGAATTCACCAGTGCTTTTGTAACTTC‑ACCCGCCATTTCTg                                                                                                                                                                                                               >  1:403486/1‑238 (MQ=255)                                 tt‑GTTCCCTCCCTTATCTGTTTTGCATTCGGGGGTTTGTTGTTGAAAAAAAGCATAGCCCAACATCGGGAGGGCAAGATTTGTTGCGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCCCTATCAACGATCAGTGATAATTCCCAACCAAAAACATCA‑TGTTTGTTTTCCGTGTAATAAGAACGTACGGTATTCACCAGATCTTTTATCATTTC‑AGCCGCCACTTCTGGc                                                                                                                                                                                                             <  1:1100991/239‑1 (MQ=255)                                      ttCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCAcc                                                                                                                                                                                                          >  1:268609/1‑238 (MQ=255)                                       tCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAAAGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCa                                                                                                                                                                                                         >  1:472460/1‑238 (MQ=255)                                        ccTCCCTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAAATTTTTTATCACTTC‑AGCCGCCACTTCTGGCAcaaca                                                                                                                                                                                                       >  2:451580/1‑235 (MQ=255)                                                tttttttttGTGTTAGGGGGTTAGTCGATAAAAAAACGCATACCACAACATCGGGAGGGCAAGTTTTGGGACGAGCAT‑CACGGAGGTTTTTTTTCCGATGGCACAGAATTGCCCCCAACACCAACCAGTAAAAATTCCAAACCACAACCTTCA‑TGTTCTTTTTCCGTGTCATAAGAACGTCCGGTATTCACCAGATTTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCAt                                                                                                                                                                                               <  2:852696/234‑1 (MQ=255)                                                aacTGTTTTGCATTAGGGGGTTAGTCGCTGAAAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCCCCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCAt                                                                                                                                                                                               <  1:1339102/237‑1 (MQ=255)                                                    gTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCCCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGCCATACGAACGGACGGTATTCACCAACTCTTTTATCCCTTC‑CGCCGCCAATTCTGGCACCCGCCAAACCACCGGc                                                                                                                                                                                           >  2:230128/1‑239 (MQ=255)                                                      tttGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGCCATAAGACCGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGCATTCGGCg                                                                                                                                                                                          >  1:825211/1‑238 (MQ=255)                                                        tGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGGAGAGCAT‑CACGGAGGTTTTTTTTGCGAGGGCGCAGAAATTGAGCCATCAACGATCAGTGATATTTACCAAACACAAAACACT‑AGTTCGTTTTCCGGGTACTAAGAACGGACGGAATTGCCCAGATCTTTTATATCTTA‑AGCAGCAACTTCTGGCACCAGCAAAGTcctgggagtac                                                                                                                                                                                       >  2:1278985/1‑229 (MQ=255)                                                        tGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGTTATTCACCAGATCTTTTATCACTTC‑CGCCGCCACTTCTTGCCCCACCAAATTCATCGGCgcc                                                                                                                                                                                        >  1:1005252/1‑236 (MQ=255)                                                        ggaaTTAGCGGGATAGTCGCTGAGAAAAAGCATAGCACAACATCGGGTGGGCAAGATTTGTGTCGAGCAT‑CACGGAGGTTTTTTTTTCGATGGCGCAGAAATTGCGCCATCAACGCTCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGGTTCATAAGAACGTACCGCATTCACCAGATCTTTTATCCCATC‑AACCGCCACTTATGCGAACAGCAAAGACAtacgggcc                                                                                                                                                                                        >  2:1278959/4‑231 (MQ=255)                                                        ggCATTAGGGGTTTAGTGGCTGATAAAAAGCAAAGCACAACATCGGGAGGGAAAGATTTGTGCCGAGCAT‑CACGGAGTTTTTTTTTGCGATGGCGCAGAAATTGGTCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGtct                                                                                                                                                                                       <  1:1148202/238‑1 (MQ=255)                                                          cATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCCGCCAAGTCATCGGCGTctct                                                                                                                                                                                     >  1:1320173/1‑239 (MQ=255)                                                               gCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGGTAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAAATCTTTTATTACTTC‑AGCCGCCACTTCTGGCACCAGACAAGCAAACGGCGTCTCTTTTTc                                                                                                                                                                                >  2:1224193/1‑239 (MQ=255)                                                               gCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTc                                                                                                                                                                                >  1:1268647/1‑239 (MQ=255)                                                                             gATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGCCCCTGTTTCATAATCGCCAGaa                                                                                                                                                                   >  1:275204/1‑238 (MQ=255)                                                                                        tAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCCCCAGCAAAGTCATCGGCGTCTCTGTTTCATAAACGACAGAAACGCCCTTGcgg                                                                                                                                                       >  2:1264359/1‑237 (MQ=255)                                                                                         aGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAATCCATCGGCGTCTCTGTTtaataaccacccaaaaacccatggcggt                                                                                                                                                      >  2:502478/1‑217 (MQ=255)                                                                                           gacaacaTGGGGGGGGAAAGTTTTGGGAAGAGCAC‑CGGGGAGTTTTTTTTGGCGATGGCGCAAAAATTGGACCATCAACGATAAGGGATAATTCCCAACCACAAACATCA‑TGTTCTTTGTCCGTGTAAAAAGAACTGACGGTATTCAACAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTa                                                                                                                                                    <  2:357642/238‑1 (MQ=255)                                                                                                     ggAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCCCCAGCAAAGTCATCGGCGTCCCTGTTTCATAATCGACAGAAACGCCATTCCTGTTATTGGGGACgg                                                                                                                                          >  2:58709/1‑239 (MQ=255)                                                                                                     ggAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGAATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGAACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCCGTTCTTGTTTACg                                                                                                                                           >  1:1103628/1‑238 (MQ=255)                                                                                                      gggggAAAGTTTTGTGACGGGCAT‑CAGGGGGGTTTTTTTTGCGATGGCGCAAAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTCCGGTGTCATAAGACCGTACGGTATTCACCAGATCTTATATCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACgg                                                                                                                                          <  1:177064/236‑1 (MQ=255)                                                                                                       gggtcAAGATTGGTAGCGAGCAT‑CACGGAGGTTTGTTGTGCGATGGCGCAGAAATTGCCCCAACAAGGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCATTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGt                                                                                                                                         <  1:1358/234‑1 (MQ=255)                                                                                                       aGGGCAAGTTTTGTGCCGA‑CATACACGGAGGTTTTTTTTGCGATGGCGCAAAAATTGCCCCACAAACGATCAGTGATAATTACCAACCAAAAACATCA‑TGTTCGTTTCCCGTGTCATAAGAACGTACGGTATTCCCAAGATTTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGt                                                                                                                                         <  1:1264359/238‑1 (MQ=255)                                                                                                        gggCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTggggcaggtac                                                                                                                                       >  2:629070/1‑230 (MQ=255)                                                                                                        gggCAAAATTTGGGACGACCAT‑CCCGGAGGTTTTTTTTGCGATGGCGAAAAAATTCCGCCACCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTc                                                                                                                                        <  1:16592/238‑1 (MQ=255)                                                                                                         ggCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACCGAAACGCCATTGCTGTTATTGGTGACGGTCAc                                                                                                                                      >  1:802846/1‑239 (MQ=255)                                                                                                          gCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACg                                                                                                                                     >  1:188702/1‑239 (MQ=255)                                                                                                                  tttgggacagaatggcggggtttttttttGGGAGGGACGAGAATGTGGCCCATCAATGGTAAGGTATCATTAACAACCAAAAACTACG‑TGGTCGTTTTCCGGGTAAAAAGACCGTGCGGTATTCCCCAGATCTTTTATCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACg                                                                                                                              <  2:1080588/218‑1 (MQ=255)                                                                                                                   gtgACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGtt                                                                                                                            >  1:781022/1‑239 (MQ=255)                                                                                                                      aCGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCCAAGTCATCGGCGTCTCTGGTTCATAATCGCCAGAAACGCCCTTGCTGTTATTGGTGACGGTACCGGGATAAGTTGCt                                                                                                                         >  1:1273364/1‑239 (MQ=255)                                                                                                                                 gggggtttttttttGCGAGGGCAAAAAAATTGCCCCAACAACGAAAATTAAAAATCCCCACCCA‑AAACCTCACTGTTTGTTTTCCGTGTCAAAAGACCTCACGTTTTTCCCAAGATCTTTTATCACTGC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGa                                                                                                               <  2:825211/233‑1 (MQ=255)                                                                                                                                         ttttttGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATa                                                                                                                                                                             >  1:349021/1‑169 (MQ=255)                                                                                                                                         ttttttGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGGATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCCAAGTCATCGGAGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGGGGAGGCAACGGGATACGTTGCTTGGCCCagggttacttt                                                                                                        >  2:806708/1‑228 (MQ=255)                                                                                                                                            tttGCGATGGCGACAAAATTGCCCCATCAACGATCATGAATAATTCCAAACCAAAAACCTCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCCCCAGATCTTTTATCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCg                                                                                                    <  1:109373/239‑1 (MQ=255)                                                                                                                                            ttggCGATGGCGCAAAAATTGCGCCATAAACAATCGTTAATACTAACCAACCACAAACATTA‑TGTTCTTTTCCCGTGTAAAAAGAACGTCCGGTATTCACAAGATCTTTTATCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCg                                                                                                    <  1:601235/236‑1 (MQ=255)                                                                                                                                                cGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCAACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACCGAAACGCCATGGCTGTTATTGGTGACGGTCACGGTATAAGTTGCTTTGCCCATTATTAATTTCCCGTTat                                                                                                >  2:220064/1‑239 (MQ=255)                                                                                                                                                             attgcccccacaaagtacAGGGAAAATTCCCAACCAAAAACATCA‑TGTTCTTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATTTTTTATCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCg                                                                                   <  1:58709/222‑1 (MQ=255)                                                                                                                                                                gcgcCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCgttg                                                                                >  1:1305524/1‑239 (MQ=255)                                                                                                                                                                gcgcCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCgttg                                                                                >  1:451902/1‑239 (MQ=255)                                                                                                                                                                     ctcacattttagtgtaaaTTCCCACCCAAAAACCTCT‑TGTTCGTTTCCCGTGTCATAAGAACGTACGTTATTCACCAGATCTTTTACCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTgcg                                                                           <  2:231333/223‑1 (MQ=255)                                                                                                                                                                     aTCAATGATCATTAATACTTACCAACCCAAAACTTCT‑TTTTCTTTTTCCGTCTCATAAGAACGTACGGTATTCACCAGTTCTTTTATCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTgcg                                                                           <  2:451902/239‑1 (MQ=255)                                                                                                                                                                              cAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTTTTGCGCACCTTCCa                                                                  >  1:842722/1‑239 (MQ=255)                                                                                                                                                                               aGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATTGGCGTCTCTGTTTCATAATCGACAGAAACACCATTGCTGTTATTGGTGCCGGTCACGGTATACGTTGCTTTGCCCCTGATTCATTTCCCGTAATGAAGGACTTTCCGTTTTTGGGCACCTTCCaa                                                                 >  2:687678/1‑238 (MQ=255)                                                                                                                                                                                gtattAATTACCAACCCAAACCTTTT‑TTTTCTTTTTCCGTGTCATAAAAACTTCCGGTTTTCACCAGATTTTTTACCCCTTC‑AGCCGCCCCTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATc                                                                <  2:1268647/235‑1 (MQ=255)                                                                                                                                                                                  taTAATTCCCAACCAAAACCTTCA‑TGTTCGTTTCCCGTGTTAAAAAACCGTCCGGTATTCCCCAGTTTTTTTACACCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAg                                                              <  1:1337732/238‑1 (MQ=255)                                                                                                                                                                                           caaccaCAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAgg                                                     >  1:780148/1‑239 (MQ=255)                                                                                                                                                                                           caaccaCAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAgg                                                     >  1:362772/1‑239 (MQ=255)                                                                                                                                                                                             accaCAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTTCGCACCTTCCATCAGGACTTAAGGa                                                    >  2:61212/1‑238 (MQ=255)                                                                                                                                                                                               cacaAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGcc                                                 >  2:337773/1‑239 (MQ=255)                                                                                                                                                                                               caaaaaCCTGA‑TGTTCGTTTCCCGTGTCATAAGACCGTCGGGTTTTCCCCAGATCTTTTAACACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGcc                                                 <  1:181023/236‑1 (MQ=255)                                                                                                                                                                                                    acaatctggttctgttctggtctcatcaGAACTTACTGTTTTCCCCAGATCTTTTATCTCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACga                                             <  1:1278959/217‑1 (MQ=255)                                                                                                                                                                                                     actctcgtttCCTTTTCCGTGTAAAAAGGCCGTACGGTATGCCCCAGTTTTTTTACCCCTTC‑AGCCGCCACTTAGGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACgaa                                            <  1:687678/231‑1 (MQ=255)                                                                                                                                                                                                         a‑tGTTCGTTCCCCGTGTAAAAAGAGCGTACGTTTTTCACCAGATTTTTTATCCGTTC‑GGCCGCCACTTCCGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGt                                       <  1:230128/239‑1 (MQ=255)                                                                                                                                                                                                             ttCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCCGGAGCCACGAAGTAGCCaa                                    >  1:560716/1‑239 (MQ=255)                                                                                                                                                                                                                   ttCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGa                              >  1:896406/1‑239 (MQ=255)                                                                                                                                                                                                                       gtgtCATACGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCCCCTTCCATCCGGACTTCCGGAGCACCGAAACAATCAACGTTGAAATa                          >  2:120782/1‑239 (MQ=255)                                                                                                                                                                                                                        tttCATAAGAACGTACGGTTTTCCCCATTTTTTTTACCCATTG‑AGCCGCCACTTCTGGCCCCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATaa                         <  2:780148/237‑1 (MQ=255)                                                                                                                                                                                                                          ttataagcagattcGGTATCCCCCAGTTTTTTTAT‑ACTTCGAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGt                       <  2:268609/226‑1 (MQ=255)                                                                                                                                                                                                                               aGACCGTTCGGTTTTCACAAGTTTTTTTACACCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCgt                  <  1:1224193/239‑1 (MQ=255)                                                                                                                                                                                                                                     tcgggTATTCCCCAGATTTTTTACCCCTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTCGAAATAAGTATCgtcagtca            <  1:502478/236‑1 (MQ=255)                                                                                                                                                                                                                                     tACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATg                                  >  1:1329561/1‑217 (MQ=255)                                                                                                                                                                                                                                     tACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGAATCgtcagcaa            >  1:294487/1‑236 (MQ=255)                                                                                                                                                                                                                                     tACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGACGTCTCTGCTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATGCATTTCCCGTTATGAATGACTTTCCGTTGTTCTGCACCTTCGACACGGCCTTCAGGATGCCCCAAGAATTTCATGTTGAAACaaattgcgctag               >  2:1191412/1‑226 (MQ=255)                                                                                                                                                                                                                                      aCGGTTTTCACCAGATTTTTTAGCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCAt           <  1:451580/239‑1 (MQ=255)                                                                                                                                                                                                                                       cGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGAGAGTCATg          >  2:976124/1‑239 (MQ=255)                                                                                                                                                                                                                                            ttCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCCCGAAGAAGTCCATGTTGAAATAAGATTCGTCAGTCATGGCTTc     >  1:559165/1‑239 (MQ=255)                                                                                                                                                                                                                                               aCCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAGCAAGTCAATGTTGACATAAGTATGGTCAGTCATGGCTTAAAt  >  1:708962/1‑239 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                AAGCCCAAAAATAATCCATGACAGAATGCCCATT‑GTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCAT‑CACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCA‑TGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTC‑AGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAAT  >  NC_000913/1879223‑1879696

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