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_0009131,061,9350GA91.4% 84.0 / ‑6.1 37S46N (AGC→AAC) torDtrimethylamine‑N‑oxide reductase‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (19/13);  minor base C (0/3);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 8.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTG‑ACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGG  >  NC_000913/1061692‑1062177
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
gtggAAATTGAGAAGTACAACGGAACAGGGGAGCAGGTG‑ACGGCGTTTAACTGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTATCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAAccg                                                                                                                                                                                                                                             <  2:25385/251‑1 (MQ=255)
        tGAGAAGTACAACGGAACAGTGGAGCAGGTG‑ACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTACCGCGTCGCAGGTGCAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCCAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGATGAAAAACGAACCGCCGCTccc                                                                                                                                                                                                                                     >  1:247214/1‑249 (MQ=255)
            ccgtcacctgctccAGTGGAGCAGGTG‑ACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACAGCCGCTCACTgcg                                                                                                                                                                                                                                 >  2:303979/14‑251 (MQ=255)
                  aaCGGACCGGTGGAGC‑GGTGGACGGGGTTTAACGTCCCCGTGGGGATGGTGGCGCAGGGCAAATATTTCCCCGGGCCGCAGGGGAAATCATGACCACGCTAAAAGCACAAAAGATTGCCTGTTTTGACGCGTGGTTGGCGCAGTTGTTCTCCCGTGAGCTGGACGATAAACACCTA‑ACGCAAACC‑GCC‑GTGGCGCATATGGCTGTATGTTGTTCTTTGATGAAAAGCAAACCGCCGCTCACTGCTGCGGTg                                                                                                                                                                                                                           <  2:181773/251‑1 (MQ=255)
                          ggTGGAGCGGGTG‑ACGGCGTTGACCGGCCCGGTGGAGAGGGGGGCGCAGTGCAAATAGGTCCCCGCGTCGCGGTTAAAACCATGCCCACGTTGACGCCACAACAGTTTGCCTGTGTTTACCCCTGGCTAGCGCATTTTTTCCCCCTGGAGCTGGACGATGAAAACCTG‑ACGCAAATC‑GCCAGGG‑CGCATATGGCTGAATGTTTTTCGTTGTAAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCt                                                                                                                                                                                                                   <  2:224107/250‑1 (MQ=255)
                                  cggTG‑ACGGCGTTTAACGGCCCCGTGGGGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGCTTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGCACGGGGTGGAAAAcc                                                                                                                                                                                                           >  2:159508/2‑251 (MQ=255)
                                  aGGTG‑ACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGTAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGGGCTGGAAAAcc                                                                                                                                                                                                           >  2:129432/1‑251 (MQ=255)
                                          gggttttAACGGCCCCGGGGAGAGGGTGGCGCAGTTCGAATATTTTCCCGCGTCGCAGGGAAAATCATGACCACTTTGACAGCCCAACAGATTGCCTGGGTTTACGCCTGGTTAGCGCAGTTGTTCCCCCGTGAGCTGGAGGATAAACAACGG‑AC‑CAAATCCGCCAGTG‑CGCAGATGTCTGATGGGTTTTCTTTGCTGAAAAACAAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGc                                                                                                                                                                                                    <  1:142647/247‑1 (MQ=255)
                                                      ccccGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGCTTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGGGAACGAACTGGAAAACCGTATTGCCCCGCttcaagg                                                                                                                                                                                        >  2:297151/1‑245 (MQ=255)
                                                      ccccGGGGAGAGGGTGGCGCAGGGCGAATATGTTCCCGCGTCGCAGGAGAAATCATGCCCCCGCTGACAGCACAACAGATTGCTTGTTTTTACGCCGGCCTAGCGCAGTTGTTCTCCCGGGAGCTGGACGATGAACAACTG‑ACGAAAATC‑GCCAGTG‑CGAAGATGGCGTAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGt                                                                                                                                                                                        <  1:317969/251‑1 (MQ=255)
                                                          gTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAAACGTATTGCCACGCTGACAGTACgg                                                                                                                                                                                    >  1:332173/1‑250 (MQ=255)
                                                               gtgggtggCGCAGTGGGAATATGTCGCGGCGTCGCAGGTGAAATGATGACCACGCTGACAGCACACCAGATGGCTGGTGTTTACGCCTGGCGAGCGCAGTTGTTCTCCGGGGGGCGGGACGATGAACAACTG‑CCGAAAACC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTGCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGa                                                                                                                                                                               <  2:100603/248‑1 (MQ=255)
                                                                aTGGTGGGGCAGTGCAAATATGTCCCCGCGCCGCAGGTGAAATCATGACCACGCTGACAGCACAACATATTGCCTGTGTTTACGCCGGTCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAAAACCTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAACCGACCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGAt                                                                                                                                                                              <  1:103499/251‑1 (MQ=255)
                                                                    tggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATAccc                                                                                                                                                                          >  1:220006/1‑251 (MQ=255)
                                                                     ggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCg                                                                                                                                                                         >  1:353792/1‑251 (MQ=255)
                                                                          agggCGAATATGTCCCCGCGGCGCAGGTGAAATCAGGCCCACGCTGACAGCACAACAGATTGCTTTTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGAGGAACAACTG‑ACGAAAATC‑GCCAGGG‑CGCAGATGGCTGAATGTTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTgg                                                                                                                                                                    <  1:391248/248‑1 (MQ=255)
                                                                                    tGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGACCTGGCcgc                                                                                                                                                          >  1:10753/1‑251 (MQ=255)
                                                                                                   ggTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCt                                                                                                                                           >  1:367305/1‑251 (MQ=255)
                                                                                                    gtaaaagcaccaCCACTCTGACAGCACAAGGTATTGCCTTTGATTACGCGTGGCACGCGATGTTGTTCTCCGGTCGGCTGGATGTTGAACAACTG‑ACGAAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTGCGTTGATAAAAAACGAGCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTg                                                                                                                                          <  1:79411/240‑1 (MQ=255)
                                                                                                           cATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGCAAAACGAACCGCCGCTCACTGCGGCGGTGAACGCGCTGGAAAACCGTATTGCCACGATGACAGGACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCtg                                                                                                                                    >  2:359013/1‑250 (MQ=255)
                                                                                                            atgcccACGCTGACAGCACAACAAATTGCCTGGGTTTACGCCTGGCTAGCGCATTTGTTCTCCCTGGAGCTGGACGATGAACAACTG‑ACGAAAATC‑GCCAGTG‑CGCAGCGGGCTGAATGGTTTTCGTTGCTGAAAATCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTGGTGCGGCCTGTTTCtgat                                                                                                                                  <  1:129432/247‑1 (MQ=255)
                                                                                                                       gACAGCACAACAGTTTCCTTGTTTGTCCGCCGGGCTAGCGCAGTTTTTCTCCCGGGAGCTGAACGATGA‑CAAATGGACGCAAATC‑GCCAGTG‑CGCAGAGGCCTGAATGTTTTTCGTTGCTAAAAAACAAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGacaaa                                                                                                                        <  1:14315/250‑1 (MQ=255)
                                                                                                                           gCACAACAGATTGCCTGTTTTTACGCCGGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCAGAAAACCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAgc                                                                                                                   <  2:367305/251‑1 (MQ=255)
                                                                                                                                    aTTGCCTGTGTTTACCCCTGGCTAGCGCAGTTGTTCCCCCGTGAGCTGGACGATTAACAACTG‑ACGCAAATC‑GCCAGGG‑CGCAGATGGCTGAATGGTTTTCGTTGCTAAAAACCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGcc                                                                                                          <  1:323182/251‑1 (MQ=255)
                                                                                                                                           ggtttcaCGCCGGGCTAGCGCAGTTGTCCTCCCGTGAGCTGGACGATAAACAACTG‑ACGCAAACC‑GCCAGTG‑CGCAGAGGGTTGAAGGGTTTTCGTTGCGAAAAAACGACCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGc                                                                                                    <  2:322727/244‑1 (MQ=255)
                                                                                                                                                      tGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGGGACGATGCCCGGCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGGGCCGCCGTATGCATCGGCCTACa                                                                                         >  1:128010/1‑250 (MQ=255)
                                                                                                                                                                 gttgttCTCCCGGGAGCGGGACGATAAAAAACTG‑CCGCAAATC‑GCCAGTG‑CGCAGATGGTTAAAGGGTTTTCGTTTCGAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGc                                                                             <  1:297151/251‑1 (MQ=255)
                                                                                                                                                                        tcCCGTGAGCTGGACGATGTACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGACCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGCCTGGAACTGGCCGCGGACCTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTCTGCCTCGGCCTACAAACAGGACGAGccat                                                                          >  1:342192/1‑244 (MQ=255)
                                                                                                                                                                                   ggccGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGCGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAACGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGt                                                           >  1:175531/4‑251 (MQ=255)
                                                                                                                                                                                   ggACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGAGGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTg                                                            <  1:28617/250‑1 (MQ=255)
                                                                                                                                                                                                                  gTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCCGCGGc                              >  1:309043/1‑251 (MQ=255)
                                                                                                                                                                                                                  gTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGGGAACGAGCTGGAAAACCGTATTGCCACGCTGACCGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGGATGCATCGGCCTACAAACAGGACGAGCAAGAGAATAAACGCTTGTTAGTTGGGGAAGGGATGGAAAcccacggg                              >  2:46604/1‑245 (MQ=255)
                                                                                                                                                                                                                             ggCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACg                    <  1:124364/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggagtgtttttCGTTGCAAAAAAACGACCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAAc                 <  2:277848/244‑1 (MQ=255)
                                                                                                                                                                                                                                            gTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATTCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGatcatt     >  1:11276/1‑250 (MQ=255)
                                                                                                                                                                                                                                             ttGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCCGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAAATCATCt    >  1:191987/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAAACGGATTGCCACCCTGACAGCACGTGACGATGCCCGTCCGGAACTGGCCGACGACTTTTGCGCCCTATTTCTGCTGACCGGCAAACCAGCGGCGCTGCCGGATGCATCGGCCTACAAACAGGACGAGCAAGAGATTATACGATTGTTATTTAAGGCCCGGATGGAAACCAACGGCAATCTCAAAGCACCGGCAGATCAACTgg  >  1:37374/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTG‑ACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTG‑ACGCAAATC‑GCCAGTG‑CGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGG  >  NC_000913/1061692‑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 < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑