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.7%	92.0 / NA	41	S46N (AGC→AAC)	torD	trimethylamine‑N‑oxide reductase‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/22);  minor base C (0/3);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAG‑ATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGG  >  NC_000913/1061691‑1062177                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    ggtggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAG‑ATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAAcc                                                                                                                                                                                                                                               <  1:352669/251‑1 (MQ=255)    ggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAG‑ATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGGTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAAccgcc                                                                                                                                                                                                                                            >  1:100418/1‑251 (MQ=255)           gagaAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAG‑ATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCTCCTGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACGGCCGCCCAcg                                                                                                                                                                                                                                     >  2:28270/1‑250 (MQ=255)             gaAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAG‑ATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑GGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGCTGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCCCTgc                                                                                                                                                                                                                                   >  2:42890/1‑251 (MQ=255)                  acatcGGAACAGAGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAG‑AGGGTGGCGCAGTGCGAATATGTCCCCGCGTCGCAGTTGAAATCATGCCCACGCTGACAGCACAACATATTGCCTG‑TTTTTC‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTG‑GCATGGACGATGACCATCTGCGGCAAATCGCCAGTGCCCAGAGGGTTGATTGGTTTT‑CGTTGCTTAAAACCGAACCGCCGCTCACTgcggcg                                                                                                                                                                                                                               <  1:342227/246‑1 (MQ=255)                     agggaaaggtgggggaGGTGACGGCGTTTACCGGCCCGGTGGAG‑ATGGTGGCGAAGTGGGAATTTGTCCCCGCGTCGAAGGGGAAATCATGCCCACGCTGACAGCAAAACAGATTGCCGG‑TTTTTA‑CGCCGGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGAAAACCGCCGTTGCGCAGATGGCTGAATGGTTTT‑CTTGGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGa                                                                                                                                                                                                                           <  1:346101/236‑1 (MQ=255)                                gagctgtggACGGCGTTTA‑CGGCCCCGTGGGGGATGGTGGCGCAGTGCGAAGAGTTTCCCGCGTCGCGGGTGAAATCTGGCCCACGCTAACAGCACAACGGTTTGCCTG‑TGTTTT‑CGCCTGGC‑TAGCGCAGTGGTTCTCCCGTGAGC‑TGGACGATGAACATCTGACGCAAAGCCCCAGTGCCCAAATGGGTGAATGTTTTG‑CGTTGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGaa                                                                                                                                                                                                                <  1:362084/243‑1 (MQ=255)                                                ttAACGGCCCCGTGGAG‑ATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGCACGAGCTGGAAAACCGTA‑TTGCCACGc                                                                                                                                                                                               >  1:230361/1‑251 (MQ=255)                                                    cGGCCCCGTGGAG‑ATGGTGGCGCAGTACGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGATGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGCGAACGTGCTGGAAAACCGTA‑TTGCCACGCTGAc                                                                                                                                                                                           >  1:433425/1‑251 (MQ=255)                                                            gggAG‑ATGGTGGCGCAGTGCAAATATGTCCCCGCGTCGCAGGTAAAATCATGACCCCGCTGACAGCACAACAGTTTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCATTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTg                                                                                                                                                                                   <  1:42890/250‑1 (MQ=255)                                                            gggAG‑AGGGTGGCGCAGTGGGAATATGTTCCGGCGCCGCAGGGAAAACCATGCCCACGCGGACGGCACAACAGATT‑CCTGGTGTTTC‑CGCCTGGC‑TAGCGCAGTTGTTCCCCCGTGAGT‑GGGACGATGACCAACTGACGAAAATCGCCGGTGCGCAGATGGTTGAATGTTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTg                                                                                                                                                                                   <  1:186927/250‑1 (MQ=255)                                                               gggagggtggggcattgtgattatGTCCCCGCGTCGCAGGGAAAATCAGGCCCACTCGGCCAGCCAAACAGAT‑GTCTG‑GGTTTAACCCTTGGT‑TAGCGCAGTTGTCCCCCCGGGGTC‑TGGAAGATAAACAACGTACCAAAATCGCCAGTGCGCAGATGGCTGAATGTTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGAc                                                                                                                                                                                 <  1:308870/230‑1 (MQ=255)                                                                         cgcAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGGACGGGCTGGGAAACCGTA‑TTGGCCCGCTTGCAGTACGGGACGATGACCggc                                                                                                                                                                       >  2:416574/1‑249 (MQ=255)                                                                            ggggcgaaatggtccccGCGCCGCAGGTGAAATCATGCCCACGCTGACAGCAAAACAGATTGCCTG‑TGTTTA‑CGCCGGGC‑TAGCGCATTTGTCCCCCCGTGAGC‑GGGCCGATGAACAACTGACGAAAATCGCCAGGGCGCAGATGGCTGAATGTTTTT‑CGTTGAGGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGGGACGATGCCCGTCTgg                                                                                                                                                                    <  1:386881/237‑1 (MQ=255)                                                                              ggCGAATATGTTCCCGCGTCGCAGGTAAAACCATGCCCCCCCTGACAGCACAACAGATTCCCTG‑TGTTAC‑CGCCTGCC‑TGGCGCAGTTGTTCCCCCGGGAGC‑GGGACGAGGAACAACGGACGAAAACCGCCGGGGCGAAGAGGGCTGAATGGTTTT‑CGTTGCAGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGaa                                                                                                                                                                  <  1:38754/250‑1 (MQ=255)                                                                                       gtccccgcttcgcagtataaaTCATGCCCCCGTTGACAGCACAACAGATTCCCTG‑TTTTTA‑CGCTTGCC‑TAGCGCAGTTGTCCCCCCGTGAGC‑TGGCCAATAAACAACTTACGCAAATCCCCAGTGCCCAGTTGGCTGATTGTTTTT‑CGTTGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCcgc                                                                                                                                                          <  1:38696/232‑1 (MQ=255)                                                                                          cccGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGAGGCCCGTCTGGAACTGGCCGCGGAc                                                                                                                                                      >  1:86108/1‑251 (MQ=255)                                                                                                       agAAACCATGACCTCGCTGACGCCACAACAGATTGCTTG‑TGTTTA‑CGC‑TGGCCTAGGGAGGTTGTTCCCCCGTGAGC‑TGGACTATGAACAACGGCCGCAAATCGCTAGTGCACAGATGGCGGAATGTTTTT‑GGTTGCTAAAAAAGGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTg                                                                                                                                          <  1:279378/249‑1 (MQ=255)                                                                                                        gAAATCATGACCACGCTGACAGCACAACAGCTTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGGCTTTTGCGGCCTGtt                                                                                                                                        >  1:375999/1‑251 (MQ=255)                                                                                                         aaaTCATGACCACGTGGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCTGTTGTTCTCCCGTGAGC‑TGGACGAAGAACAACGGACGAAAACCGCCAGTGCGCAGAGGGCTGAAGGGTTTT‑CGTTGCTGAAAACCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGAGGCCTGttt                                                                                                                                       <  1:236233/251‑1 (MQ=255)                                                                                                                gACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCCCCCGTGAGC‑GGGACGATGAACAACTGACGAAAATCGCCAGTACGCAGATGGCTGAATGTTTTT‑CGTTGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCtgatga                                                                                                                                <  1:336340/251‑1 (MQ=255)                                                                                                                   cACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGGCGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCtgatgaacg                                                                                                                             >  1:328049/1‑248 (MQ=255)                                                                                                                       cTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGGAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATTCCCGCCTGGAACTGGCCGCGGACTTTTGTGGGCTGTTTCtgcttgccgcca                                                                                                                          >  2:205225/1‑240 (MQ=255)                                                                                                                                    ggatggcttg‑tgTTTC‑CGCCGGGC‑TAGCGCAGTGGTCCCCCGGGGGGC‑GGGACGATGAAAAACTGAGGCAAACCGCAAGTGCGCAGATGGGTGAATGTTTTT‑CGTGGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTg                                                                                                            <  2:235791/243‑1 (MQ=255)                                                                                                                                       ttGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGa                                                                                                                                                                                                                 >  1:327503/1‑148 (MQ=255)                                                                                                                                                   tA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGc                                                                                              >  1:194167/1‑251 (MQ=255)                                                                                                                                                          tGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACa                                                                                         >  1:25006/1‑250 (MQ=255)                                                                                                                                                             gcgaGCGCAGTTGTGCCCCGGTGGGC‑GGGACGATGAAAAACTGACGCAAATCGCCAGTGCGCAGAGGGTTGAATGGTTTT‑CGTGGCGGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAAc                                                                                      <  1:248500/248‑1 (MQ=255)                                                                                                                                                                   acAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTCAACGAGCGGGAAAACCGTA‑TTGTCACCCTGCCAGTACGTGAAGATTCCCGTCTGGAACTTGCAGCGGCCCTTTTCGGCATGGGTCTGATACCCGCCAAACACGCGGCGCCTACGGATGCATCGGCCTACAAACaatagc                                                                                >  2:257174/2‑246 (MQ=255)                                                                                                                                                                        gtttcccccGGGAGC‑TGGACGATGAACAACGGACGCAAACCCCCAGTGCGCAGAGGGTGGATTGGTTTT‑CTTTGCAAAAAAACAAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCaa                                                                           <  1:416574/245‑1 (MQ=255)                                                                                                                                                                          gttctcccgggagcggaaagagaaacaacTGACGCAAAGCGCCAGGGCGCAGAGGGCGGAAGGGTTTG‑CGTGGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAag                                                                          <  1:83952/228‑1 (MQ=255)                                                                                                                                                                               ccGTGAGC‑TGGACGATAAACAACTGACGAAAATCGCCAGGGCGCAGTTGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTa                                                                    <  1:130949/251‑1 (MQ=255)                                                                                                                                                                                      gctgtaggataaaaactggACGAAAATCGCAAGTGCGCAGAGGGCTGAATGTTTTC‑CGTTGCGAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCt                                                              <  1:298521/233‑1 (MQ=255)                                                                                                                                                                                               tGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCCTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGGACTGGCCGCGGGCTTTTGCGGCCTGTTTCTGATGAACGGCAAACAAGCGGCGCTGCCGTATGCATCGGCCTCCAAAACGGACGAGCAAAAGGTTGAAAGCTTTTTAGTTg                                                     >  2:451697/1‑251 (MQ=255)                                                                                                                                                                                                     gctggcgtaattccgtaGTGCGCAGTTGG‑TGAATGGTTTTTAGTGGGTAAAAACCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGACCTGGCCGCGGACTTTTGCGGCCTGTTGCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAg                                               <  2:166433/235‑1 (MQ=255)                                                                                                                                                                                                                   ccAGTGCGCAGATGGCTGAATGGTTTG‑CTTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGc                                 <  1:184156/251‑1 (MQ=255)                                                                                                                                                                                                                       ggcgcAGATGGGTGAGTGGTTCC‑CGTTGAGGAAAAACGACCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCa                             <  2:120606/250‑1 (MQ=255)                                                                                                                                                                                                                             gATGGCTAAAGGGTTTT‑CGTGGCGGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGt                                                           <  1:336696/215‑1 (MQ=255)                                                                                                                                                                                                                                   ggattGGTTTT‑CGTTGCAGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAAc                 <  1:375925/247‑1 (MQ=255)                                                                                                                                                                                                                                                 ttGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTAGTTACCA‑GCTGACAGTACGTGCAGATGTCCGTCTGGAACGGGCCGACGACTTTTGCGGCCAGGTTCTCATGCACGACGAACAGACGGCGCTGCCGGCAGCAGCTGCCTACCGACAAGGCGAGCAAACCATTGAACGCGTGGTAGTGGAGGCAGGACTTGAAACCAGCGGCAAATGTACCGAACCGGCAGaccctcc    >  1:296543/1‑245 (MQ=255)                                                                                                                                                                                                                                                   gtgaaaaaACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTgg  <  2:86108/247‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAG‑ATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTG‑TGTTTA‑CGCCTGGC‑TAGCGCAGTTGTTCTCCCGTGAGC‑TGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTT‑CGTTGCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTA‑TTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGG  >  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: ‑