Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,040,129 A→C 19.1% T218P (ACA→CCA)  yedY → membrane‑anchored, periplasmic TMAO, DMSO reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,1290AC19.1% 111.8 / 2.4 56T218P (ACA→CCA) yedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base A (14/29);  new base C (6/6);  total (21/35)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAA  >  NC_000913/2039894‑2040362
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
gCGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTcaca                                                                                                                                                                                                                                             <  1:93252/238‑1 (MQ=255)
 cGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCGTTATCGGCGGTGGGCGGAAAAATCCTTTTGTCGAAGGGGTGTGTTTCGGCGGAGGAATAGATCCGCTCacacc                                                                                                                                                                                                                                           >  1:514380/1‑238 (MQ=255)
    atgtaacGTTTGCGCCGCGCGGCAGAGTGGGCGGTTGTGGTGCCGTGGATTGTTCATGCGGTGCACAAATTGCGGGCGCTTGCCGAACC‑TACGAGCAATGCGAAGTATGGCGCTTCCGAAACAA‑GTTAGGCACCGGAA‑CA‑GATGCCGGGCCAGCAGGACCGCTTTATGGGCGGTGGGCTGAAATATCCTTATGTCGAGGGATTGCGTCTGGACGAAGCAATGCAGCCGATCACACtgca                                                                                                                                                                                                                                        <  2:831364/233‑3 (MQ=255)
    aTTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgat                                                                                                                                                                                                                                        <  1:144214/239‑1 (MQ=255)
         tCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑Cg                                                                                                                                                                                                                                 <  1:163060/239‑1 (MQ=255)
         tCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCGGAAATATCCTTATTTCGGAGGATTGTGTTTCCGCGGAGGAATACAACCCCTCACACttat‑gaaag                                                                                                                                                                                                                                  >  1:63743/1‑236 (MQ=255)
             aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAGATACTTATGGCGAAGGGAGTCGTCTCGACGGAACCATGAAACCGCCAACACagat‑gagccttgg                                                                                                                                                                                                                              >  1:485206/1‑232 (MQ=255)
             aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAg                                                                                                                                                                                                                              <  1:322227/238‑1 (MQ=255)
             aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGGGGAATATCCTTTTGTCGAAGGAGTGCGTCTCCACGCCGGCATGCATCCCCTCCCACCGAT‑GTG‑CGTAgg                                                                                                                                                                                                                             >  2:554658/1‑239 (MQ=255)
                   tGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTTTCTGAGGGGTGGGTCTCCGCGCAGCAAAGCATCCCCTCTCACTGAT‑TAC‑Cgggtgggtttt                                                                                                                                                                                                                       >  2:695712/1‑229 (MQ=255)
                   tGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAAAATACTTATGTCGAAGGATTTCGTCTCGCCGAAGCAATGCATCCGCCCCCACTGTT‑GAC‑CCTAGGGG‑GTTa                                                                                                                                                                                                                      >  1:885367/1‑239 (MQ=255)
                     cGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCCTTTTTCGCAGGATTTCGTCTCGACGAAACAAAGCATCCCCTCACACTGAT‑GGC‑CCTAGGAG‑GTTAtt                                                                                                                                                                                                                    >  1:889083/1‑238 (MQ=255)
                        ggAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGc                                                                                                                                                                                                                  <  1:961168/238‑1 (MQ=255)
                         gAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCa                                                                                                                                                                                                                 <  2:427109/238‑1 (MQ=255)
                          aaGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCaa                                                                                                                                                                                                                <  1:1153249/238‑1 (MQ=255)
                                  tCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGCCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgat‑gtgagcgg                                                                                                                                                                                                                               <  1:368676/217‑8 (MQ=255)
                                    gATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAATATCCTTATGTCGGAGGGTTGCGTCTCGACGAAGCAATGGATCCGCCCCCACTTCT‑TAA‑CGTAGGGG‑TTTATGGCAAGGCGGTACCg                                                                                                                                                                                                      >  1:1304048/1‑238 (MQ=255)
                                            tgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAACATACTTATTTCGAAGGGTTTGGTCGCCACGAAGCAAGGCAACCCCCCACACTGCT‑GAG‑CCTAGGGG‑GTTTTGGGAAAGCGTTATAGCCACcaaaa                                                                                                                                                                                             >  2:826643/1‑238 (MQ=255)
                                                    gttGGTTTCCACCTGCACAAATTGCGGGCTC‑TGCCCAACCTTACAAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTGTGCACCGGAA‑CA‑GAGCCCAGGCCAGAGGGACCGCTTTAGCGGCGGTGGGCTGAAATATCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GCC‑CGTAGGTG‑TTTATGGAAAGGTGTTACCGCCACAAAATGgcgcgc                                                                                                                                                                                      <  2:1351486/237‑1 (MQ=255)
                                                          tttCCGCGGCACAAATTGCTGGCGCTTGCCGACCC‑TACCAGAAAGGCGAAGTTGGTCGCTTCCGAAACAA‑TTTATGCCCCGGAA‑CA‑GATGCCAGGCCGGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATGGCGTCTCGACGAAGCAATGCACCCGCTCACACTGAG‑GAC‑CGTGGGTG‑TTTAGGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCg                                                                                                                                                                               <  2:711761/239‑1 (MQ=255)
                                                                       aaTTGCTGGCGCTTGCCGAACC‑TACCGGAAATGCGAATTATGTCGTTTTCGAAACAT‑TTTATGC‑CCGGAACCA‑GATGCCAGGCCAGCAGGCCCGCTTTACCGGGGGTGGGTTAAATTATCTTTATGTCGAGGGATTGCGTCTCGACGAAGCAAGCCATCCGCTCCCATGGAT‑GCC‑CGTAGTTT‑TTTATGGCAAGCCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCg                                                                                                                                                                 <  2:1275795/239‑1 (MQ=255)
                                                                                cgctgcccGAACC‑TCCCAGAAATGCGAGGTATGGCGTTTTCAAAAAAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGCCCGTTTTATCGGGGGTGGGCTGAAATACCCTTATGTCGAAGGATGGCGTCTCGACGAAGCAATGCACCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAata                                                                                                                                                         <  2:1304048/232‑1 (MQ=255)
                                                                                cgcTGGCCAACCC‑TACCAGAAAGGCAAAGTATGTCG‑TTTCGAAAAAATTTTAGGCCCCGGAA‑CA‑GATGCCGGGCCAGAAGGCCCGCGTTATCGGGGGGGGGCTAAAATATCCTTATGTCAAAGGTTGGCGTCTCGACGAAGAAATGCATCCGCTCACCCTGAT‑GAC‑CGTAGGTG‑TTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAata                                                                                                                                                         <  2:514380/238‑1 (MQ=255)
                                                                                               aCCAGTAATGGTAAGTAGGTCGCTTTAGAAAAAA‑TTTAGGCACCGGAA‑AG‑GAGGCCAGGCGGGAGGTACCTCTTTATGGGTGGGGGGATGAAATATCCTTAGGTGGAAGGATGGCTCCTCGACGAAGAAAGGCATCGGCTCACATTGAT‑GCC‑GGTAGGTT‑TTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGAt                                                                                                                                          <  2:370213/239‑1 (MQ=255)
                                                                                                      cgcggAAGTAGTTCGTTTTCAAAAAAA‑TTTATGCACCGGAA‑CG‑TGTGCCGGGCAGGGAGGCCCGTTTGTCCGGCGGGGGTCTGAAAAATCTTTATGTCGAAGGTTTGCGTCTCGACGAAGAATTGCTTCCGCTCCCGCTGAT‑GCC‑CGTGGGTG‑TTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCg                                                                                                                                   <  2:145301/235‑1 (MQ=255)
                                                                                                      agtggAAGTATGTCCCTTTCAAAACAT‑TTTATGCACCGAAA‑CA‑GATCCCAGGCGAGCAGGCCCTCTTTCGCGGCGGTGGGCTGAAATATCTTTATGTCGAAGGTTTGCGTCTCGACGAAGAAATGCATCCGCTCACACTGAG‑GAC‑TGTAGGTT‑TTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCg                                                                                                                                   <  2:615172/235‑1 (MQ=255)
                                                                                                      aTGCGAAGTATGCCGTTTTGAAAACAA‑TTTATGCACGGGAA‑CA‑GATGCCAGGCCAGCAGGACCGTTTTAGCGGGGGGGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCGCGACGAAGAAATGCATCCGCTCACACTGAT‑GAC‑CGTGGGTG‑TTTAGGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCg                                                                                                                                   <  2:1291526/239‑1 (MQ=255)
                                                                                                       tGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATCCTTTTTTTGGAGGGATGGGTTTCGGCGGAGGAATGCATCCCCTCCCCCCGGT‑GAG‑CGGAGGGG‑GTTTTGGGAAAGGGGTAACGCCCCCAAAAGGCGCGGCGGGGGGGCTGCTTTTTGCGGGGGAAAATTGCTTTTAAAGGGGTaaaaaga                                                                                                                                   >  2:394296/1‑236 (MQ=255)
                                                                                                       tGCGAAGGATGTCGCTTTCAAAACAT‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCGGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTGGGTG‑TTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGa                                                                                                                                  <  1:826643/239‑1 (MQ=255)
                                                                                                           aaGTATTGCGCTGTCTAAACAA‑TTTATCCCCCGAAA‑CA‑GATGCCGGCCCGGCAGGACCGTTTTATCGGTGGTGGGCTAAAATACCCTTATGGCGAAGGATTGCGTCTCGACGAGCCACTGCACCCGCTCACACGGAG‑GAC‑CGTAGGTG‑TTTATGCAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:1269567/239‑1 (MQ=255)
                                                                                                            ggTATGTCGTGTTCGAAACAA‑TTTTTGCACCGGAA‑CA‑GATGCCAGGCCGGCGGGCCCGCTTTAGCGGCGGTGGGCTAAAATATCTTTATGTCGAGGGATTGCGTCTCGACAAAACAATGCACCCGCTCACACGGAT‑GAC‑TGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:917108/237‑1 (MQ=255)
                                                                                                            attATGCCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATCTCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  1:1305774/236‑1 (MQ=255)
                                                                                                            aggaTGTGGCTTGCAAAACAA‑TGTATGCGCCGAAC‑CA‑GATCCCGGGCCACGGGGACCGCTTTCCCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGTTTGCGTCTCGGCGAAGCAATGCACCCGCTCACACTGAT‑GAC‑GGTGGGTT‑TTTAGGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:349639/235‑1 (MQ=255)
                                                                                                            aGTATGTCGCTTGCGAAACAA‑TTTATGCCCCGAAA‑AA‑GATGCCAGGCCAGGAGGCCCTCTTTATCGGCGGTGGGCTGAAATTTCTTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  1:502751/238‑1 (MQ=255)
                                                                                                            aGTATGGCGCTTTCAAACCAA‑TTTATGCGCCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTGTGGGCGGTGGGCTGAAATATCCTTATGGCGAAGGATGGCGTCTCGACGAAGAAATGCATCCGCGCCCACTGAG‑GAC‑CTAAGGTT‑TTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  1:281454/238‑1 (MQ=255)
                                                                                                            aGTATGCCGTGTTGGAAACAA‑TTTATGCCCCGAAA‑CA‑GATCCCAGGCCAGCAGGCCCGCTTTATCGGCGGTGGGCTAAAATATCCTTATTTCGAAGGATTGCTTCTCGACGAAGCAATGCACCCGCTCACACTGAT‑GAC‑CGGAGGTT‑TTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:9506/238‑1 (MQ=255)
                                                                                                                tGTCGCTTTCAAAAAAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCGGAAATATCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCCCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTa                                                                                                                         <  1:806743/239‑1 (MQ=255)
                                                                                                                      tttCGAAACAA‑TTTATGCCCCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                                                                                                                   <  1:868977/239‑1 (MQ=255)
                                                                                                                      ttagaaaaCAA‑TTTATGCACC‑GAA‑CAGGATGCCAGCCCACCAGAACCGCTTTATGGGGGGTGGGCTGAAATATCTTTATGTCGAAGGATGGGGTCTCGAGGAAGCAATGCACCCGCTCACACTGAT‑GAC‑CGTGGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                                                                                                                   <  2:29910/234‑1 (MQ=255)
                                                                                                                        tCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGCGGGAATAACCCTATGTCGAAGGACTTCGTCCCGCCGGAGGAAAGCCACCGCTCACACTGCT‑GACAGGTGGGGG‑TTTATGACAAGCGGTTACCG‑CACAAAAAAGCGCGGCGCGGCGCGCG‑GTTGTGGCGGTGGAATATTGGTTTCAAGGGGTTAAAACCAAAGTCTGCAtttaattg                                                                                                                 >  2:1127620/1‑233 (MQ=255)
                                                                                                                        acGAAACAA‑TTTGTGGCCCGGAT‑GA‑GATGCCTGGCCAGGCGGACGTGTTTGTAGGGGGTTGGCAGAAATATCTTTATGAGTAGGTATCGCGTCTCGCGAAAGCAAGCCCTCCCCGCCCAGTTAT‑CCC‑CATAGGCG‑TGTATGCGAAGGCTTAACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTg                                                                                                                 <  2:53289/238‑1 (MQ=255)
                                                                                                                                          ccGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAATCCTTTTTTCGGAGGATTTCGTCTCCACGGAGCAAAGCATCCCCCCCCACTGTTGGAC‑CGAAGGGT‑TTTTTGGAAAGG‑GGTACCCCCCCAAAAAAGGGCGCCGGGGGGGCTG‑GGTGTGGCGGGGGGAAATGGCTTTAAAGGGGATAAATAGATAAACAGGTTTAAAGTGACCcgcgcgggttccc                                                                                                 >  2:197477/1‑229 (MQ=255)
                                                                                                                                                     atgcagggccggCAGGACCGTTTTACCGGCGGTGGGCGAAAATATCCTTATGTCGAAGGATGGGTTCTCGACGAAGCAATGCACCCGCTCCCACGGTT‑GCC‑CGTAGGTT‑TTTATGGCAGGGCGTTCCCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAaccacc                                                                                         <  2:219254/231‑1 (MQ=255)
                                                                                                                                                                                   ggTGGGCTGGAAAATCCTTATGTCGAAGGGTTGCGTGTCGACGAAGCAATGCAACCGCTCCCCCTCCT‑GAC‑CGTAGGGG‑TTTTTGGCCAAGCGGTACCGCCCCAAAAAAGCGCGCCGGGGCGGCTG‑TTTGTGCCGCGGGAATATGGCTTTATAGGGGTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGGCCCCCAACCCCCCGGGATCTCGCAGCGGCTGACGCAAACGGgtt                                                         >  2:595253/1‑239 (MQ=255)
                                                                                                                                                                                            aaaTATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGGGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCa                                                 >  2:397743/1‑238 (MQ=255)
                                                                                                                                                                                                ctCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAGTGCATCCGCTCCCACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                            <  1:554658/238‑1 (MQ=255)
                                                                                                                                                                                                 tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGAG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                           >  2:1025057/1‑239 (MQ=255)
                                                                                                                                                                                                   ctatggTCGAAGGATTGGGTCCCGACGAAGCAATCCACCCCCCCACACTGAT‑GAC‑GGTAGTTT‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                         <  1:394296/234‑1 (MQ=255)
                                                                                                                                                                                                      aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCCCCACAAAATGGCGCGCCGGTGCGCCTG‑ATTGTGCCGTGGACATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCACGCGTCCGCCAACAACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                      >  2:324902/1‑239 (MQ=255)
                                                                                                                                                                                                            aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                >  2:844960/1‑239 (MQ=255)
                                                                                                                                                                                                            aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                >  1:1245871/1‑239 (MQ=255)
                                                                                                                                                                                                                  tGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCa                           >  1:292242/1‑238 (MQ=255)
                                                                                                                                                                                                                  tGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAAACACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACCTTAATCCGTATGTTGATCa                           >  2:368232/1‑238 (MQ=255)
                                                                                                                                                                                                                        gcgacgaGACAATGCACCCGCTCACCCTGAG‑GCC‑GGTGGGTT‑TTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCgcgc                    <  2:889083/238‑1 (MQ=255)
                                                                                                                                                                                                                         cgacgaAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCcgc                      >  2:568578/1‑236 (MQ=255)
                                                                                                                                                                                                                                  aaTGCATCCGCTCACCGTGAG‑GAC‑CGTAGGTT‑TTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcaca            <  1:814387/237‑1 (MQ=255)
                                                                                                                                                                                                                                         ccGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGa   >  2:585763/1‑239 (MQ=255)
                                                                                                                                                                                                                                          cGCTCACATGGAG‑GAC‑GGT‑GGTGTTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGaa  <  2:485206/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
GCGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACC‑TACCAGCAATGCGAAGTATGTCGCTTTCGAAACAA‑TTTATGCACCGGAA‑CA‑GATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAT‑GAC‑CGTAGGTG‑TTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTG‑ATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAA  >  NC_000913/2039894‑2040362

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