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	92.5%	135.9 / 0.4	56	S46N (AGC→AAC)	torD	trimethylamine‑N‑oxide reductase‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (24/25);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATC  >  NC_000913/1061691‑1062174                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ggtggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAAcc                                                                                                                                                                                                                                                  <  1:258810/251‑1 (MQ=255)    ggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGAGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAAAAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGCAAAACGAAcc‑gcc                                                                                                                                                                                                                                              >  2:109689/1‑251 (MQ=255)                    aaCGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAAcc‑gcc                                                                                                                                                                                                                                              >  1:496970/1‑235 (MQ=255)                    aaCGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGGGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAAcc‑gcc                                                                                                                                                                                                                                              <  2:496970/235‑1 (MQ=255)                      cGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGa                                                                                                                                                                                                                             >  1:311741/1‑250 (MQ=255)                            ggTGGAGCGGGGGACTGCTTTTAACGGCCCCGTGGAGATGGGGGCGCAGTGAGAATATGTTCCCGCGTCGCAGGGGAAATCAGGCCCACGCTGAAGGCACAACAGATGGCTGGTTTTTACGCCTGGT‑AAGCGCAGTTGTTCTCCGGGG‑AGCGGGACGATGAACAACTGACGCAAATCGCCGGTGCGAAGAGGGCTGAATGGTTTGCGTTGCTAAAAACCGATCC‑GCCGCTCACTGCGGCGGTGAACG‑AGCt                                                                                                                                                                                                                     <  1:465688/250‑1 (MQ=255)                             gTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGCTGAAACACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAGGGTTTTCGTGGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGGACC‑AGCTg                                                                                                                                                                                                                    >  2:463312/1‑251 (MQ=255)                               ggAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGGGGACG‑GGCTggg                                                                                                                                                                                                                  >  2:169892/1‑250 (MQ=255)                                gAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGa                                                                                                                                                                                                                  >  1:444857/1‑250 (MQ=255)                                  gCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGGGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGATCTCCCGTG‑CGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAGGGTTTTCGTTGATGAAAAACGAACC‑GCCGCTCACTGCGGCGGTgtagcggctggaaa                                                                                                                                                                                                                 >  2:428042/1‑237 (MQ=255)                                      gtgggggggttaaccGGCCCCGTGGAGAGGGGGGCGAGGGGCAAATATGTCCCCGGGCCGCAGGTAAAATCAGGCCCCCGCGGACAGCAAAACAGTGTGCTTGTTTTTCCGCCTGGT‑TAGCGAGGTTGTTCCCCCGGG‑AGCTGGACGATAAACACCTGACGAAAATCGCCAGGGCGCAGATGGCTGAATGGTTTTCGTGGCGAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGt                                                                                                                                                                                                           <  2:266823/237‑1 (MQ=255)                                         cgggCGTTAAACGGCCCCGGGGAGAGGGGGGCGCAGCGCGAATATGTCCCCGGGGCGCAGGTAAAATCAGGACCCCGCTGACAGCACAACAGATGGCCTGTTTTTACGCCTGGC‑TAGCGCAGTTTTCCCCCCGGG‑AGCTGGACGATGAACAACTGACGCAAACCGCCAGTGCGCAGAGGGCTGAATGGTTTTCGTGGCAGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTAtt                                                                                                                                                                                                        <  2:356244/249‑1 (MQ=255)                                                  aaCGGCCCCGTGGAGGTGGTGGCGCAGTGCGAATATGTTCCCGCGCCGCAGGTAAAATCATGACCACGCGGACAGCACAACAGATTGCCTGTTTTTACGCTGGGC‑TAGCGCAGTTGTTCCCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGAAGATGGCTGAATGGTTTTCGTTGCTAAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTg                                                                                                                                                                                              <  1:449500/251‑1 (MQ=255)                                                         ccggggggaTGGTGGCGGAGGGCAAATATGTTCCGGGGTCGCGGGTGAAATCATTACCACGCTGACAGCACAACAGAGTGCCTGTGTTTGCGCCTGGC‑TAGGGCGGTTGTTCCCCCGGG‑AGCTGGACGAGGAAAAACTGAGGAAATTCGCTAGTGCGCGGAGGGCTTAGTGGTTTTCGTTGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTAc                                                                                                                                                                                       <  1:268469/244‑1 (MQ=255)                                                         ccGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑GGCTGGAAAACCGTATTGCC‑CCGCTGACAGTAc                                                                                                                                                                                       >  1:212807/1‑251 (MQ=255)                                                             ggAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑CCGCTGACCGTACGTg                                                                                                                                                                                    >  1:273421/1‑250 (MQ=255)                                                               agaTGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCAGGCCCACGCTGCCAGCACAACAGAGTGCCTGGTTTTACGCCGGGC‑TAGCGCATTTGTTCTCCCGTG‑AGCGGAACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGTTTTTCGTTGCGAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTTGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACg                                                                                                                                                                                 <  1:169892/251‑1 (MQ=255)                                                                  tggtggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCAGCCACCACAACAGATTGCCTGTGTTTAAGCCTGGC‑TACCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGCCGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGCCAGTCCGTGACGATg                                                                                                                                                                              >  1:72343/1‑251 (MQ=255)                                                                  tggtggCGCAGTGCGAATATGTTACCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCAGGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGGGAGCG‑AGCTGGAAAACCGTCTTGCC‑ACGCTGCCAGTACgagaacatg                                                                                                                                                                              >  2:530854/1‑243 (MQ=255)                                                                        cgcAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATAGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGGTGCCCGTc                                                                                                                                                                        >  1:109828/1‑251 (MQ=255)                                                                                  atatTTTCCCGCGTCGCAGTGGAAATCATGCCCACTCGTACAGCAAAACAGATGGCCGGTGTTTACGCCGGGC‑TAGGGCAGTTGTTCGCCCGGG‑AGGTGGAAGATGAACACGTGAGGAAAATCTCCGGGGCGCAGATGGCTGAATGGTTTTCGTTGCTAAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGc                                                                                                                                                              <  2:410240/251‑1 (MQ=255)                                                                                  atatGTTCCCGCGTCGCCGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGGGTTTACGCCTGGC‑TAGCGCAGATGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAGAACGAACC‑GCCGCTCACTGCGGCGGTGTACGCAG‑TAGAAAACCGGATTGCCCAC‑CTGACAGTACGTgcaggtggccggcgggaaactgc                                                                                                                                                              >  2:369123/1‑229 (MQ=255)                                                                                    tgtgtccggcgcggagggTGAAATCAGGCCCACGTTGCCAGCAAAAGAGATTGCCGGTGTTTACGCCTGGC‑TGGCGGTGGTGTTCTCCCGTG‑ATCGGGAGGAGGAAAAACTGAAGCAAATCGCCGGTGCGAAGGGGGCGGAAGGGTTTTCGTTGCGAAAAAACGACCC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGcc                                                                                                                                                             <  2:224963/234‑1 (MQ=255)                                                                                    atGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGGGAAAC‑AGCTGGGAAACCGTATTGCC‑ACGCTGACCGTACGTGGCGGTGCCCGTCTGGGACAGGccc                                                                                                                                                            >  2:146890/1‑250 (MQ=255)                                                                                      gTTCCCGCGGCGCAGGTAAAACCATGCCCACGCTGACAGCAAAACAGATTGCTTGTGTTTACGCCTGGC‑TTGGGCAGTTTTTCTCCCTGG‑AGCGGGAAGAGGAAAAACTGACGCAAATCGCCAGTGCGCAGAGGGCTGAATGTTTTGCGTGGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCcgcg                                                                                                                                                          <  1:109689/251‑1 (MQ=255)                                                                                       ttCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATCGCCTGTGTTTCCGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTAACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGCGAACG‑AGCTGGAAAACCGGATTGCC‑ACGCTGACAGAACGTGACGATGCCCGCCTGGAACTGGCCGCgg                                                                                                                                                         >  1:341672/1‑251 (MQ=255)                                                                                           cgcgTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCTTTGGTTTACGCCGGGC‑TAGCGGAGTTGTTCCCCCGTG‑AGCGGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGAGGGCTGAAGGTTTTTCGTTGCAAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACtt                                                                                                                                                     <  1:499973/251‑1 (MQ=255)                                                                                              gTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGAGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTg                                                                                                                                                  >  1:420119/1‑251 (MQ=255)                                                                                                                    cGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑CCGCTGACCGTACGGGGCGATGCCCGGCTTGAACTGGCCGCGGACTTTTGCGGCCTTGTTT‑TGCT‑GACCGac                                                                                                                          >  2:132233/1‑251 (MQ=255)                                                                                                                             cacagtaggttcccgttgtTTACGCCGGGC‑AAGCGATGGTGTTCTCC‑GTGGAGCTGGACGATGAACACGTGACGAAAACCGCCGGCGCGCAGATGGCTGAAGGGTTTGCGTTGTAAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAgcg                                                                                                                 <  2:105366/235‑1 (MQ=255)                                                                                                                                  caGATTGCCTTTGTTTCCGC‑TGGCCTTGCGCAGTTGTCCTCCCGTG‑AGCTGGCCGATGAACAACTGACGCAAACCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCt                                                                                                            <  1:426125/251‑1 (MQ=255)                                                                                                                                      ttgccgggttttACGCCTGGC‑TAGAGCAGTTTTTCTCCCGTG‑AGCTGGACGATGAACACCTGACCCAAACCGCCAGTGCGCAGATGTCTGAATGGTTTTCGTTCCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGcc                                                                                                         <  1:516137/241‑1 (MQ=255)                                                                                                                                           tgtgTTTCCGCCTGGC‑TAGCGCAGTTGTTCCCCCGTG‑AGCTGGACGATGAACAACTGACGAAAATCGCCAGTGCGCAGAGGGCTGAATGGTTTTCGTTGCTGAAAACCGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TAt                                                                                                    <  1:175146/249‑1 (MQ=255)                                                                                                                                                   cGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑CGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGGATTGCC‑ACGCTGACAGTAAGTGACGATGCCCGGCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCCCTGCCG‑CATGCCTCGGCCt                                                                                          >  1:490285/1‑251 (MQ=255)                                                                                                                                                             ggcgcAGTTGTTCTCCCGTG‑GGCGGGACGATGAACAACTGCCGAAAATCGCCGGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAgg                                                                                 <  1:90132/250‑1 (MQ=255)                                                                                                                                                             ggcgcAGTTGTTCCCCCGGG‑GGCTGGACGATGAACAACTGACCAAAAGCGCCAGTCCGCAGTTGGCTGAATGGTTTTGGTTGC‑GAAAAACCAACCCGCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAg                                                                                  <  1:63825/249‑1 (MQ=255)                                                                                                                                                             aGCGCAGTTGTTCTCCCGTG‑AGCTGGGCGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACCAACAgg                                                                                 >  1:69772/1‑251 (MQ=255)                                                                                                                                                                ggtttgtttcccccGTG‑AGCTGGACGATGAACAACTGACGCAAACCGCCTGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACg                                                                              <  1:431253/240‑1 (MQ=255)                                                                                                                                                                gcAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGGACG‑AGCTGGAGAACCGTATTGCC‑ACGACGACAGTACGTTACGATGCCCGTCGGGGACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GGCCGACAAACAAGCGGCGCCGCCGGTA‑GCATCGGCCTACAAACAAGACg                                                                              >  2:487255/1‑251 (MQ=255)                                                                                                                                                                  aGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGGAAACCGTATTGAC‑ACGCTTCCAGTAAGTGACGCTGCCCGTCTGGCACCTGCCGCGGCCTTTTGCGGCCT‑GTTTCTTATTGCCCG‑CAAACAAGCTGCGCTGCCC‑GATGCCGCGGCCTTCAAACaagaag                                                                              >  2:327934/1‑244 (MQ=255)                                                                                                                                                                           cccGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGGACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGGTGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCtgat‑ggacggc                                                                                                                          >  2:70103/1‑191 (MQ=255)                                                                                                                                                                              gTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTTCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTaa                                                                 >  1:580093/1‑250 (MQ=255)                                                                                                                                                                               tG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGCCTGAATGGTTTTCGTTGCTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTTAT‑TACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCCACAAACAGGACGAGCAAGAGATTaaa                                                                >  2:48522/1‑250 (MQ=255)                                                                                                                                                                                   gcgggaagaggaacaccggACGCAAACCGCAAGTGCGCAGAGGGCTGAATGGTTTGCGTTACTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCt                                                            <  1:177565/233‑1 (MQ=255)                                                                                                                                                                                   gCTGGACGATAAACAACTGACGAAAACCGCCAGTCCGCAGATGTCTAAATGTTTTCCGTTGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCt                                                            <  1:149272/251‑1 (MQ=255)                                                                                                                                                                                             gAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTg                                                   >  1:416672/1‑250 (MQ=255)                                                                                                                                                                                              aacaacTGACGAAAACCGCAAGTGCGCAGAGGGCTGAATGGTTTTCGTTGCGAAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAg                                                 <  1:426766/251‑1 (MQ=255)                                                                                                                                                                                                 accggACGAAAATCGCGAGTGCGCAGAGGGTTAAATGGTTTTCGTTGCAAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCa                                              <  1:151061/247‑1 (MQ=255)                                                                                                                                                                                                      aacaaaatggccggcgcgcAGATGGCTGAAGGTTTTTCGTTGCTGAAAAGCAAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGa                                         <  1:408901/235‑1 (MQ=255)                                                                                                                                                                                                              ccccAGTGCGCAAATGGCTAAATGTTTTTCTTTTCTAAAACCCAAACC‑GCCGCTCCCTGCGGCGGTGAACG‑AGCTAGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCa                                <  1:496232/249‑1 (MQ=255)                                                                                                                                                                                                                          gagggCTGAATGTTTTTCGTTGCGAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCAGCGGCAATTTCa                    <  1:275942/248‑1 (MQ=255)                                                                                                                                                                                                                          gATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGGTGGAAAACCGAATTGCC‑ACGCTGACAGTACGGGACGATGCACGTATGGAACAGGCCGCGGACTTTTGGAGCCT‑GTTTCTGAT‑GGAAGAAAAACAAGCGGCGCTGCCG‑GATGGATCGGCCAACAAACAGGACGGGCCAAAGATTAAAAGCTTATTAGTT‑AGGGACGGGATGGAAAACAGCGGCAATTTc                     >  1:560023/1‑250 (MQ=255)                                                                                                                                                                                                                                 gAATGGTTTTCTTTGCGAAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCAGCGGCAATTTCAACGAAcc             <  1:48522/251‑1 (MQ=255)                                                                                                                                                                                                                                   aTGGTTTTCGTTGATGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCAGCGGCAATTTCAACGAACCgg           <  1:146890/251‑1 (MQ=255)                                                                                                                                                                                                                                         ttCGTTGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGat      <  1:368876/250‑1 (MQ=255)                                                                                                                                                                                                                                            gtggctaaaaaACGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGatcatc  <  1:28066/244‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑AGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACG‑AGCTGGAAAACCGTATTGCC‑ACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCT‑GTTTCTGAT‑GACCGACAAACAAGCGGCGCTGCCG‑TATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCA‑GGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATC  >  NC_000913/1061691‑1062174

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