Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,612 A→T 17.4% intergenic (+86/+454) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,6120AT17.4% 138.3 / 2.8 57intergenic (+86/+454)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (14/32);  new base T (9/2);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 4.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGGTGTGGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAA  >  NC_000913/1286379‑1286840
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
cGATGGTGTGGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATTGGCACGGGTGGGGGGGGGTTGTCGGGTGCGGTTCCCTTTTCC‑gg                                                                                                                                                                                                                                   >  2:644300/1‑238 (MQ=255)
    ggttttgtgTTTTTCTTCCGTCTGCCCCTGGTGCTGGGGATGACGTTATTCCTGCGGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTa                                                                                                                                                                                                                               <  2:554757/231‑1 (MQ=255)
        tgGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGGGCGGGGGGAGTGTTCGGTTGCGCTTCGCTTATCC‑GACCTCCAgg                                                                                                                                                                                                                           >  1:19203/1‑239 (MQ=255)
         ggCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAggg                                                                                                                                                                                                                          <  1:644300/239‑1 (MQ=255)
          gCGTTTATCTTCCGTCTGCCCCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAggg                                                                                                                                                                                                                          <  1:729527/238‑1 (MQ=255)
               tATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCTGGGTGAGTTGTCGGGTGCGCCTCCCTTTTCC‑GGCCTTCAGGgggggag                                                                                                                                                                                                                    >  1:252787/1‑238 (MQ=255)
               tATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGAC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTTTCGGAAGCGCTTCGCTTATAC‑GACCCAAAGGggcggga                                                                                                                                                                                                                    >  1:170335/1‑237 (MQ=255)
                aTCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGTTGCGCTTCCCTTTTCC‑GACCTACAGGGGAGGatt                                                                                                                                                                                                                   >  1:554757/1‑238 (MQ=255)
                        tCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTTGAGGATCTGACACGTACGTACCAGCTGGTGCGCGCCCGTCACTAAGCGCATTTTAGTTCACATAGACCCGGCTTCGGCGGGGGGTTTTTCATGGTC‑CGGTCCGGGGGGAGGGGTGGGATGCGGTTAGCTTTACC‑GACCGACGAGGGGGGATAGTgaaagg                                                                                                                                                                                                           >  2:605451/1‑234 (MQ=255)
                          tGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGGGCTTCGCTTTTCC‑GACCTCAAGGGGGCGGTTTTTTAGGccc                                                                                                                                                                                                         >  1:846128/1‑239 (MQ=255)
                             cccgtgTGCTGGGGATTGCGTTACTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCgg                                                                                                                                                                                                       <  1:488988/233‑1 (MQ=255)
                                   tgttGGGGGTGGCGTTTTTCCTGCGTTTCCCGTTTTCGCGTCTGATACACATCTGGGGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGGGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACg                                                                                                                                                                                                 <  2:756076/235‑1 (MQ=255)
                                             gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACAGTGCGGGGTGAGTTGTCGGATGAGCTTCGCTTTTCC‑GACCTACAGGGGGGGTGATTGTGGGGCCGGGAAACGGGGTG‑cccaag                                                                                                                                                                                     >  2:441219/1‑235 (MQ=255)
                                                gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGgga                                                                                                                                                                                                                        <  1:777293/202‑1 (MQ=255)
                                                          gtgTTCCCGTTCTCGGGTTTTATACACATCGGGAGGGTACCGGTGGAGTATCGGACACGTAAGTACCAGCGGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCt                                                                                                                                                                       <  2:511813/238‑1 (MQ=255)
                                                          ctgTTCCCGTTTCCGCGTTTGGTTCCCATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCt                                                                                                                                                                       <  2:222010/239‑1 (MQ=255)
                                                          ctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATTGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCCCTTTTCC‑GACCTACAGGGGGGGATATTGTAGGCCCGGTAAACCGGGGG‑CCA‑TCGGGCAGAATTAGc                                                                                                                                                                        >  1:23639/1‑238 (MQ=255)
                                                           tgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTCGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCt                                                                                                                                                                       <  1:826164/238‑1 (MQ=255)
                                                            gTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCg                                                                                                                                                                     <  1:778441/239‑1 (MQ=255)
                                                                       cgcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACCGGGGAGGATATTGTAGGCCCGGGAAACGTGGGG‑CCT‑TCGGGCCAAATTTGCTCGGGGGAACAggg                                                                                                                                                          >  1:696775/1‑238 (MQ=255)
                                                                            agagacagcaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                                                                                                                                                     <  1:52031/231‑1 (MQ=255)
                                                                                     ctttGGGGCGTTCCGGTGGAGTATCTGAACCGTAAGTACCAGCTGGTGGGCGCTCGTCACTAATGGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg                                                                                                                                            <  2:751763/236‑1 (MQ=255)
                                                                                     aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAAc                                                                                                                                                                                                  <  1:207874/187‑1 (MQ=255)
                                                                                     aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCt                                                                                                                                                                       <  1:195136/212‑1 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGGGTTGTCGGATTCGCTTCGCTTATCC‑GCCCCCACGGGGGGGGTATTGTAGGCCCGGGGAAAGTGGGGCCCA‑TCGGACAAAATT‑GCCCGAGGTAAACGGTTTTAATTTTTTAAGgaagggaagg                                                                                                                                  >  2:232938/1‑239 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTCTTTTTATCGGCACGGTGCGGGGTGCGTTGTCGGGGGTGCTTCGCTTTTCC‑GACCTACCGGGGGGGGTTTTGTAGGCGCGCGCAACGCGGTG‑CCA‑TCGGGCACAATTAGCCCGGGGGAACAGGTTTTTTTTATTTGAAGAAGGGagag                                                                                                                                  >  2:777848/1‑236 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCC‑GGCCTACAGGGGAGGATATTGTAGGCCCCGTAAACGTGGGG‑CCA‑ACGGGCAAAATTAGATCGGGGTAACAGGGTTTGATTATTTAAGaaatggtcct                                                                                                                                  >  2:52031/1‑232 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAAATTTAGTTCACCTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGCGAGTTTTCGGCTGGGCCTCGCCTATCCGGACCAACA‑GGGAGGATATTTTAGGCCCGGTCAACGTGGTG‑GCACTCGGCGAAATTTACCGCGG‑GTAAACGGTTTTGTTTTTTTAAGGAAGCggggg                                                                                                                                  >  2:731739/1‑235 (MQ=255)
                                                                                               tACCGGGGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtg                                                                                                                                   <  1:884447/238‑1 (MQ=255)
                                                                                                aCCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCt                                                                                                                                                                       <  1:191758/201‑1 (MQ=255)
                                                                                                  cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggtgg                                                                                                                               <  1:232938/239‑1 (MQ=255)
                                                                                                   gggtgggTATTTGACACGTAAGTTCCAGCTGGTGCGCGCTAGTCACTAAGGGAATTTTAGTTGGCATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggtggt                                                                                                                              <  2:23639/233‑1 (MQ=255)
                                                                                                         gtttctgcaccGTTAGTACCCGCTGGTGCGCGCTCGGTATTAAGGGAATTTTAGTTCACATAGGCCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTggggg                                                                                                                         <  2:252787/228‑1 (MQ=255)
                                                                                                         gTATCTGACAGGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTggggg                                                                                                                         <  1:631088/238‑1 (MQ=255)
                                                                                                             cTGACACGTAAGTCCCAGCTGGGGCGCTGGCGTCCCTAAGTGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGa                                                                                                                    <  1:6885/239‑1 (MQ=255)
                                                                                                                    gTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGATTATTTAAGGAAGCGCTGGGGGTGGGGGAAGGGTTACTCa                                                                                                             >  1:576086/1‑239 (MQ=255)
                                                                                                                                   tGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCCCGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTTTCC‑GCCCTCCAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCACAATTAGACCGGGGTAACCGGTTTTGATTATTTAAAGAAGGGGTGGTGGTGGGGGGAGGAGTTATCAGCGCcgcgcggtcg                                                                                               >  2:209370/1‑234 (MQ=255)
                                                                                                                                         gtCGTCCCTAAGCGAATTTTAGTTCCCATAGGCCCTGATTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                                        <  1:605451/238‑1 (MQ=255)
                                                                                                                                                      gAATTTTAGTTCACATGGACCCTTCTTCGGC‑GGGTTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                                                             <  1:506148/237‑1 (MQ=255)
                                                                                                                                                               gtCACATAGGCCCTGCTTCGGC‑GGGGTTTTTTTAGGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                  <  2:696775/238‑1 (MQ=255)
                                                                                                                                                                            tGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGgc                                                      >  1:512163/1‑238 (MQ=255)
                                                                                                                                                                             gCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                                                                                                                                                     >  1:521045/1‑142 (MQ=255)
                                                                                                                                                                             gCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                                                                                                                                                     <  2:521045/142‑1 (MQ=255)
                                                                                                                                                                                  ggc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAg                                               <  1:209370/239‑1 (MQ=255)
                                                                                                                                                                                         gTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCCACGCCCTCTCCCTGGCCCCCAAGTCGaa                                          >  1:6513/1‑238 (MQ=255)
                                                                                                                                                                                          tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAc                                         <  1:557829/238‑1 (MQ=255)
                                                                                                                                                                                                taggggCACGGTGGGGGGGGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTc                                  <  1:441219/236‑1 (MQ=255)
                                                                                                                                                                                                                    aGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGa                                                                                                                                     <  1:858761/120‑1 (MQ=255)
                                                                                                                                                                                                                    aGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGa                                                                                                                                     >  2:858761/1‑120 (MQ=255)
                                                                                                                                                                                                                       tggcGGATGCGCTTGTCTTTTGC‑GACCTACAGGGGAGGATTTTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCGCGGAGTAACTGGTTTGGATTATTGAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc            <  2:691956/235‑1 (MQ=21)
                                                                                                                                                                                                                       tGTCGGGTTGGCTTCGCTTTTCC‑GACCGACAGGGGAGGGTTTTGTTGGCCCGGTAAACGGGGTG‑GCA‑GCGGGAAAAATTGGCTCGGAGTAACGGGTTTTGATTTTTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc            <  2:585404/238‑1 (MQ=18)
                                                                                                                                                                                                                           ggtctcgcttcgtttATCC‑GGCGTAGAGGGGAGGATTTGGTAGGCCTGGTTAACGTTGTG‑CCA‑TCGGGCAATATTAGCTCGGAGTAACAGGTTTGTATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTggg       <  2:250809/234‑1 (MQ=21)
                                                                                                                                                                                                                              tGCGCTTGGCTTATCC‑GACCTAAGGGGGGTGATTTTGGAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCa     <  1:908161/238‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCGGGTAAAGGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCa     <  1:410381/238‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGGACCTTGGTCGAAACTTCTCATCCTTCCCCGCTTGGGCAg    >  1:585404/1‑239 (MQ=255)
                                                                                                                                                                                                                              tGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAg    >  1:133795/1‑239 (MQ=255)
                                                                                                                                                                                                                                cgcttcgcttATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTCCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCTCAAGTCGAACCTTGGCCGAACCTTTCCTTCCTTCCCCCCTTTGGCGGag  >  1:202456/1‑238 (MQ=25)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
CGATGGTGTGGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCC‑GACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTG‑CCA‑TCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAA  >  NC_000913/1286379‑1286840

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