Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,061,935	G→A	S46N (AGC→AAC)	torD →	trimethylamine‑N‑oxide reductase‑specific chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,061,935	0	G	A	88.0%	117.5 / ‑5.4	51	S46N (AGC→AAC)	torD	trimethylamine‑N‑oxide reductase‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (25/19);  minor base C (1/5);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 9.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGATGG‑AAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGG  >  NC_000913/1061691‑1062177                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         gggggAAATTGAGAAGGAAAACGGAACAGTGGAGCCGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAAGCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCGATCGCAGTTGTTCTCCCGTGAG‑CTGTACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGaa                                                                                                                                                                                                                                                      <  2:578512/246‑1 (MQ=255)             gaAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTgc                                                                                                                                                                                                                                       >  1:270067/1‑251 (MQ=255)                         aaCAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGCACAACTGCCGCAAATCGCCCGTG‑CGCC‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCTCCCACTGCGGCGGttgacg                                                                                                                                                                                                                            >  2:285137/1‑245 (MQ=255)                                gAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑a                                                                                                                                                                                                                    >  1:168739/1‑250 (MQ=255)                                   cAGGTGAGGGCGTTTAACGGCCCGGGGGAGAGGGGGGCGCAGTGCGAATATGTCCCCGC‑GTCGCGGGTGAAATCATGCCCACGTTGACAGCACAACAGATTGCTTGT‑TTTTACGCCGGGCTAGCCCAGTTGTTCTCCCGTGAG‑CGGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAAc                                                                                                                                                                                                                <  1:65473/251‑1 (MQ=255)                                    aGGTGACGGAGTTTAACGGCCCCGAGGAGATGGTGGCGCAGAGCGAATAAGTTCCCGC‑CTCGCAGGGGAAATCATGACCACGCTGACAGCACAACCGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGCG‑CTGGCCGATGAACAACTGCCGCAACTCGCCCGTG‑CGAC‑GATGGCTGAG‑CGGTTTTC‑G‑TTGCTGAAAACCGAACC‑GCCGCTCCCTCcggcggcgtaacagcagggaaacc                                                                                                                                                                                                                >  2:680089/1‑233 (MQ=255)                                                   ccGGCCCCGGGGAGATGGTGGCGCAGTGCGAATATGTTCCCCC‑GTCGCAGGTAAAATCATGCCCACGCGGACAGCACAACAGATTGCCTGT‑TTTTACGCCTGGCAAGCGCAGTTGTTCCCCCGTGAG‑CGGGACGATGAACAACTGACGAAAAGCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑T‑TTGCTGAAAACCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTg                                                                                                                                                                                                 <  1:682804/249‑1 (MQ=255)                                                   aCGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGCCCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGCACGAGGTGG‑AAAAACGAATTGCCACgacgg                                                                                                                                                                                                >  2:682804/1‑247 (MQ=255)                                                      gCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGGACGAGCTGG‑AAAACCgattttccacgcgtac                                                                                                                                                                                               >  2:196826/1‑234 (MQ=255)                                                         ccGTGGAGATGGGGGCGCGGTGCGTATATGTCCCCGC‑GTCGCAGGTAAAATCATGCCCACGCTGACAGCACAACAGATTGC‑TGTTGTTTACGCCTGGCTGGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGTTTTGC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCGGG‑AAAACCGTATTGCCACGCTGACAGTAc                                                                                                                                                                                          <  1:251315/251‑1 (MQ=255)                                                          cGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAAAGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACg                                                                                                                                                                                         >  1:416957/1‑251 (MQ=255)                                                          cGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGATGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTTGACGAGCTGGGAAAAcctctttcccccctgccagggcg                                                                                                                                                                                          >  2:31821/1‑230 (MQ=255)                                                               agagggtggCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCGGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GAGGGTGGAA‑TGGTTTTC‑G‑TTGTTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACg                                                                                                                                                                                    <  1:273318/247‑1 (MQ=255)                                                                   ggaggCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGAAAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTGTCCCGTGAG‑CTGGACGATGAACAACTGACGCCAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAGCGAACC‑GCCGCTCGCTGCGGCGGTGAACGGGCTTG‑AAAACCGTATTGCCACGCTGACAGTACGGGACGATGc                                                                                                                                                                                >  1:48916/4‑251 (MQ=255)                                                                       gctcAGTGCGAATTTGT‑CCCGCCGCCGCAGGGAAAATAAGGACCCCGCTGACAACACAACAGATTGCCGGT‑GTTTACGCCGGTCTAGGGCAGTTGTTCTCCCGGGAG‑CTGGAGGATAAACAACTGACGCAAAGCGCCAGGG‑CGCA‑GATGGCTGAA‑GGTTTTTC‑G‑TTGCTAAAAAACGACCC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGt                                                                                                                                                                            <  1:277716/248‑1 (MQ=255)                                                                               caaaTATGTCCCCGC‑GTCGAAGGTAAAATCAGGCCCCCGCGGACAGCACAACATATTGCCTGT‑GTTTCCGCCGGGCTAGCGCAGTTGTCCCCCGGGGAG‑CGGGACGATGAACACCGGACGAAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGTTTTTC‑G‑TTGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACt                                                                                                                                                                    <  1:285137/249‑1 (MQ=255)                                                                                gAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGGGACGATGCCCGTCTGGCACTg                                                                                                                                                                   >  1:336092/1‑251 (MQ=255)                                                                                  atatGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGAGGGGAACGAGCTGG‑GAAACCGTTATGCCCAGCTGCCAGTACGTGGAGATACCCGTCTGGaccaggc                                                                                                                                                                 >  2:255397/1‑245 (MQ=255)                                                                                   tatGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAAAGAGATGG‑GAAACAGTATTGCCACGCTGACAGAAAGTGACGATTGCCGTCTGGAACTGGcc                                                                                                                                                                >  1:666022/1‑251 (MQ=255)                                                                                          ccGC‑GCCGCAGGTAAAACCATGACCACGCGAACAGCAAACCAGATTGCCTGT‑TTTTCCGCCTGGCAAGCGCAGTTTTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGGG‑CGCA‑GATGGCTGAA‑‑GGTTTTCCG‑TTGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACt                                                                                                                                                         <  1:278346/251‑1 (MQ=255)                                                                                               gTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTg                                                                                                                                                                                                 >  1:320482/1‑207 (MQ=255)                                                                                                      gtaaaaaattgacccggCTGCCAGCCCAAAAGATTGCTTGT‑GTTTACGCCTGGCAAGCGAGGTTTTTCTCCCGTGGG‑CTGGCCGAGAAACAACTGAAGAAAACCGCCAGGG‑CGCA‑GATGGTTGAA‑GGGTTTTC‑G‑TGGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTg                                                                                                                                              <  1:139875/235‑1 (MQ=255)                                                                                                        gAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGGGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGGTGCCCGTCTGGGACTGGCCGCGGACTTTTGCGGCCTGtt                                                                                                                                            >  1:129118/1‑251 (MQ=255)                                                                                                          aatctgaaccgcgcggACAGCACACCAGATTGCTTGT‑GTTTACGCGGGGCGAGCGCAGTT‑TTCCCCCCTGGGTCGGGACGATAAACAATTGACGCAAATGGCC‑GTGGCGCA‑GATGGTTGAA‑TGGTTTGC‑G‑TTGTTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTc                                                                                                                                          <  1:433305/236‑1 (MQ=255)                                                                                                           aTCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCAAGCTGACAGTACGTGACGAAGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCt                                                                                                                                         >  1:455764/1‑251 (MQ=255)                                                                                                                        tGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCACGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAAAGAGCTGG‑AAAACCGTATTGCCACGCTGCCAGTAAGTGACGATGCCCGTCTGGGAATGGCCGCGGACTTTTGCGGCCTGTTTCTTAT‑TACCGacaa                                                                                                                            >  1:286227/1‑250 (MQ=255)                                                                                                                                   caGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCt                                                                                                               >  1:147405/1‑251 (MQ=255)                                                                                                                                    cggtgggctggtttttccGCCTGGCTAGCGCAGTTGTTCTCCCGGGAG‑CGGGACGATAAACACCGGACGCAAACCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAACCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCt                                                                                                               <  1:255397/234‑1 (MQ=255)                                                                                                                                    aagtgtgcttgtttttccGCCTGGCTGGGGGAGTTTTTCCCCCGTGAG‑CGGGACGATGAACACCTGACGCAAATCGCCAGTG‑CGCA‑GA‑GGGTGAAATGGTTTTC‑G‑TTGGTGAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCt                                                                                                               <  1:350264/234‑1 (MQ=255)                                                                                                                                          ccTGT‑GTTTCCGCCGGGCTAGCGCAGTTGTTCTCCGGGGAG‑CTGGACGATGAAAAACGGACGCAAATCGCCGGTG‑CGAA‑GGTGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTAt                                                                                                        <  1:388614/251‑1 (MQ=255)                                                                                                                                                              aGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGAA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑CAAACCGGATTGCCACGCTGACCGTACGTGACGATGCCCGTCTAGAACTGGTCGCGGAACTTTACGGCCTGTTTCTTATTGAC‑GACAAAAAA‑ACGGGGCTGCCGTAATCATCgggcccaaaaaa                                                                                       >  2:536319/1‑239 (MQ=255)                                                                                                                                                                cgaaGTTGTTCTCCCGGGAG‑CTGGAGAATAAACAACTGACGAAAATCGCCGGGG‑CGCA‑GATGGCTGAT‑TGGTTTTC‑G‑TTGCGGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGAc                                                                                   <  1:421529/248‑1 (MQ=255)                                                                                                                                                                       gttCCCCCGTGAG‑CTGGACGATGAACAACTGACGCAAGGCCCCAGTG‑CGCA‑GATGGTGGAA‑TG‑TTTTC‑GGTTGCTAAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCaa                                                                             <  1:147413/250‑1 (MQ=255)                                                                                                                                                                         tctcCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTTCGGCCTGTTTCTGTT‑GACCGACAAACAA‑g                                                                                                                      >  1:417039/1‑207 (MQ=255)                                                                                                                                                                             cctgtgagtgggactagaaacaacTCACGCAAATCCCCAG‑G‑CCCAAGAGGGCTAAA‑TGTTTTTC‑G‑TTGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAgaga                                                                         <  1:105726/231‑1 (MQ=255)                                                                                                                                                                                gtgagtggtaggataaacatcttaCGCAAAGCGCCTGGG‑CGCA‑GAGGGCTGAA‑TGTTTTCC‑G‑TTGCGGAAAACCAAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTaaa                                                                    <  1:529329/228‑1 (MQ=255)                                                                                                                                                                                          cGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAAAGCTTGtt                                                            >  1:234311/1‑249 (MQ=255)                                                                                                                                                                                              gAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCt                                                                                                               >  1:678623/1‑194 (MQ=255)                                                                                                                                                                                                aaatcggccgaaactCGCCAGGG‑CGCA‑GATGGCTGAA‑TGGTTTTT‑G‑TTGCTAAAAAACAACCC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGTCTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑Agg                                                   <  1:440660/237‑1 (MQ=255)                                                                                                                                                                                                    cTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGCACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGGATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGccgg                                               >  1:597443/1‑248 (MQ=255)                                                                                                                                                                                                        gccgaaaatcgccggcgtgaaaagggCTGGA‑TGGTTTGC‑G‑TGGC‑GAAAAACGGACCCGCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGa                                             <  2:23626/227‑1 (MQ=255)                                                                                                                                                                                                          aaaaccGCCGGTG‑CGCA‑GATGTCTGAA‑TGTTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGATGG‑a                                        <  1:15080/246‑1 (MQ=255)                                                                                                                                                                                                          aaaaTCGCCAGTG‑CGAA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTAAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGATGG‑a                                        <  1:196826/250‑1 (MQ=255)                                                                                                                                                                                                                             gATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGGATTGCCCCGCTGACAGTACGGGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGAAAAACAA‑GCGGCGCTGCCGTAAGCATCGGCCTACAAACCAGCCGAGCAAGCGGTTAAAACCTTTTTAGTTG‑GGGCAGGGGTGGGAAACC‑GCGGGAAtttt                        >  2:330196/1‑249 (MQ=255)                                                                                                                                                                                                                                 gCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTCCGGGACGCTGCCCGGCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGCT‑GACCGGCAAACAAAGCG‑CGCTGCCGTGAGCATCGGCCTACAAACAGGGAGAGACAGGGATTAAAAGATTGTTTGTTGGAGG‑AGGGGTGG‑GAAACAGCGGGAATTTccaagg                   >  2:394294/1‑246 (MQ=255)                                                                                                                                                                                                                                  cTGAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCCCGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAAAAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGCGATTAAACGCTTGTTAGTTG‑TGGCAGGGATGG‑AAACCAGCGGCAATTTCAACGa                   >  1:474206/1‑250 (MQ=255)                                                                                                                                                                                                                                   gggtgatggttttccgggggaaaaaaCCGACCC‑GCCGCCCACTGCGGCGGGGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGGACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGATGG‑AAACCAGCGGCAATTTCAACg                    <  2:168739/229‑1 (MQ=255)                                                                                                                                                                                                                                    gAA‑TGGTTTTC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAAAGAGCTGG‑AAAACCGGATTGCCACGCTGACAGTACGTGACGATGACCGTCTGGAACTGGACGCGGACTTTTGCGGACTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCCTCGGCCTACAAAAAGGACGAGCAAGAGGTTAAAAGCTTGTTAGATG‑AGGCAGGGATGG‑AAACCAGCGGCAATTTCAACGAAcc                >  1:103300/1‑251 (MQ=255)                                                                                                                                                                                                                                         ggTTTTC‑G‑TGGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCGGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACGGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGATGG‑AAACCAGCGGCAATTTCAACGAACCGGCa            <  1:31821/251‑1 (MQ=255)                                                                                                                                                                                                                                             ttC‑G‑TTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGAAGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCCTCGGCCTAAAAAAAGGACGGGCAAGAGATTAAACGCTTGTTAGGTG‑TGGCGGGGGTGG‑GAACACGCGGGAATTTACACGCACCGGCacggc        >  2:107542/1‑247 (MQ=255)                                                                                                                                                                                                                                                     gctaaaaaCCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGATGG‑AAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTgg  <  1:536319/247‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGC‑GTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGT‑GTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAG‑CTGGACGATGAACAACTGACGCAAATCGCCAGTG‑CGCA‑GATGGCTGAA‑TGGTTTTC‑G‑TTGCTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGG‑AAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGAT‑GACCGACAAACAA‑GCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTG‑AGGCAGGGATGG‑AAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGG  >  NC_000913/1061691‑1062177

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