Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 94,078 G→A 14.5% G305S (GGC→AGC)  murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091394,0780GA14.5% 1808.1 / 236.2 686G305S (GGC→AGC) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (309/276);  new base A (59/40);  total (369/317)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 696 total aligned reads displayed.

CGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGT  >  NC_000913/93930‑94226
                                                                                                                                                    |                                                                                                                                                    
cgcgGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATgg                                                                                                                                                      <  3:1360532/149‑1 (MQ=255)
  cgGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGgc                                                                                                                                                    <  2:3139340/149‑1 (MQ=255)
  cgGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGgc                                                                                                                                                    <  4:532183/149‑1 (MQ=255)
    gTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGgcgc                                                                                                                                                  >  1:3032321/1‑149 (MQ=255)
       gCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCttt                                                                                                                                               >  1:1582574/1‑149 (MQ=255)
        cGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCtttt                                                                                                                                              <  1:4611882/149‑1 (MQ=255)
           tATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAAc                                                                                                                                           <  3:253037/149‑1 (MQ=255)
           tATCAATGGAAGATCATATTAATCCGAACTGTAGCGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAAc                                                                                                                                           <  4:4717286/149‑1 (MQ=255)
               aaTGGAAGATCATGTTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGATTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACTTCa                                                                                                                                       <  4:4523358/149‑1 (MQ=255)
               aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCa                                                                                                                                       <  4:4885764/149‑1 (MQ=255)
               aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCa                                                                                                                                       >  1:2491413/1‑149 (MQ=255)
                aTGGAAGATCATATTAATCCGAACCGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCtt                                                                                                                                                <  3:4490383/139‑1 (MQ=255)
                  ggAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCa                                                                                                                                    >  2:1133257/1‑149 (MQ=255)
                   gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCaa                                                                                                                                   >  2:2901177/1‑149 (MQ=255)
                   gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCaa                                                                                                                                   <  3:3104165/149‑1 (MQ=255)
                   gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCaa                                                                                                                                   <  2:6886595/149‑1 (MQ=255)
                         cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACctgct                                                                                                                             <  4:1019819/149‑1 (MQ=255)
                             ttAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTc                                                                                                                         >  2:4425512/1‑149 (MQ=255)
                             ttAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACCATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTc                                                                                                                         <  3:3340356/149‑1 (MQ=255)
                                aTCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCTCGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTcgcg                                                                                                                      >  4:3861402/1‑149 (MQ=255)
                                 tCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCgcgc                                                                                                                     <  3:788315/149‑1 (MQ=255)
                                 tCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCgcgc                                                                                                                     <  3:4011440/149‑1 (MQ=255)
                                 tCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCAGCTGCTCgcgc                                                                                                                     <  3:6158383/149‑1 (MQ=255)
                                   cGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTg                                                                                                                   >  3:6319390/1‑149 (MQ=255)
                                     aaCTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGc                                                                                                                 <  4:6917753/149‑1 (MQ=255)
                                       cTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGa                                                                                                               <  1:3791600/149‑1 (MQ=255)
                                          tCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGacac                                                                                                            >  4:1600367/1‑149 (MQ=255)
                                          tCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGACATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGacac                                                                                                            >  4:443449/1‑149 (MQ=255)
                                             cGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGt                                                                                                         >  4:1977578/1‑149 (MQ=255)
                                             cGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGt                                                                                                         <  4:2105473/149‑1 (MQ=255)
                                             cGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGt                                                                                                         >  1:4066057/1‑149 (MQ=255)
                                              ggACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGtt                                                                                                        <  3:2324368/149‑1 (MQ=255)
                                               gACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTg                                                                                                       <  3:3718847/149‑1 (MQ=255)
                                               gACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTg                                                                                                       <  2:1999035/149‑1 (MQ=255)
                                               gACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTg                                                                                                       >  3:198048/1‑149 (MQ=255)
                                                       ttGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCgg                                                                                               >  2:1639508/1‑149 (MQ=255)
                                                        tGAAAGCGACCGAAGTGAACTATCACGACCGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGc                                                                                              >  1:5779892/1‑149 (MQ=255)
                                                         gAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTACTGCTCGCGCTGGCGACACTGTTGGCACTCGGCt                                                                                             >  4:6579217/1‑149 (MQ=255)
                                                         gAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCt                                                                                             <  2:7273631/149‑1 (MQ=255)
                                                          aaaGCGACTGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTa                                                                                            >  2:4610350/1‑149 (MQ=255)
                                                           aaGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTAt                                                                                           >  3:5742169/1‑149 (MQ=255)
                                                            aGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATc                                                                                          <  3:5289673/149‑1 (MQ=255)
                                                            aGCGACCGAAGTGAACTATCACGACACGGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTGGCGCTGGCGAAACTGTTGGCAATCGGCTATc                                                                                          >  2:966279/1‑149 (MQ=255)
                                                             gCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATcc                                                                                         >  4:2992365/1‑149 (MQ=255)
                                                                aCCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACt                                                                                      >  2:2237288/1‑149 (MQ=255)
                                                                   gAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGc                                                                                   <  2:6210348/149‑1 (MQ=255)
                                                                    aaGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCt                                                                                  <  1:631546/149‑1 (MQ=255)
                                                                    aaGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCt                                                                                  <  2:4240754/149‑1 (MQ=255)
                                                                    aaGTGAACTATCAAGACAGCGGTGCGACGATTCGCTGTAGCTCAAGTTGGGGCGATGGCGAAACTGAAAGCCATCGGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCt                                                                                  >  2:6547155/1‑149 (MQ=255)
                                                                      gTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGa                                                                                >  1:2860197/1‑149 (MQ=255)
                                                                       tGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGAt                                                                               >  1:6678141/1‑149 (MQ=255)
                                                                       tGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGAt                                                                               <  2:4006682/149‑1 (MQ=255)
                                                                           cTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATctgc                                                                           >  1:2080506/1‑149 (MQ=255)
                                                                           cTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCCGCTGCTCGCGCTGGc                                                                                                                 <  1:5084794/111‑1 (MQ=255)
                                                                           cTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCCGCTGCTCGCGCTGGc                                                                                                                 >  2:5084794/1‑111 (MQ=255)
                                                                            tATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATctgct                                                                          >  3:4828185/1‑149 (MQ=255)
                                                                             aTCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATctgctg                                                                         >  3:1688929/1‑149 (MQ=255)
                                                                              tCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGa                                                                        >  1:984170/1‑149 (MQ=255)
                                                                                          tGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCCCTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCcgcgc                                                            <  4:4128556/149‑1 (MQ=255)
                                                                                            cgacgaTTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGtt                                                                                                        >  2:2476581/1‑103 (MQ=255)
                                                                                            cgacgaTTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGt                                                          <  1:2357160/149‑1 (MQ=255)
                                                                                             gacgaTTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTc                                                         >  3:689060/1‑149 (MQ=255)
                                                                                             gacgaTTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTc                                                         >  3:4177373/1‑149 (MQ=255)
                                                                                             gacgaTTCGCTTTAACTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTc                                                         >  3:549207/1‑149 (MQ=255)
                                                                                              acgaTTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCt                                                        >  4:2545162/1‑149 (MQ=255)
                                                                                                  ttCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCACTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCaa                                                    >  3:6586109/1‑149 (MQ=255)
                                                                                                       tttAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCGACCGGt                                               <  3:4575912/149‑1 (MQ=255)
                                                                                                       tttAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGt                                               >  3:2786051/1‑149 (MQ=255)
                                                                                                       tttAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGt                                               >  1:1007141/1‑149 (MQ=255)
                                                                                                         tAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGttt                                             <  1:328022/149‑1 (MQ=255)
                                                                                                         tAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGttt                                             <  4:3941723/149‑1 (MQ=255)
                                                                                                          aGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTg                                            >  1:4485535/1‑149 (MQ=255)
                                                                                                           gCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGc                                                                                   <  1:5626643/109‑1 (MQ=255)
                                                                                                           gCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGc                                           <  2:1339835/149‑1 (MQ=255)
                                                                                                           gCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGc                                           >  2:4846314/1‑149 (MQ=255)
                                                                                                           gCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAATCCGCCGCGCGTCTGCAACCGGTTTGc                                           <  2:6677440/149‑1 (MQ=255)
                                                                                                             tCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTc                                                         <  2:3365652/133‑1 (MQ=255)
                                                                                                               aaGTTGGGGTGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGAc                                       <  2:5250967/149‑1 (MQ=255)
                                                                                                               aaGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGAc                                       <  1:1873198/149‑1 (MQ=255)
                                                                                                                   tGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTAt                                   <  4:1228528/149‑1 (MQ=255)
                                                                                                                    ggggCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATg                                  <  4:5928930/149‑1 (MQ=255)
                                                                                                                       gCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGaa                               <  4:5742169/149‑1 (MQ=255)
                                                                                                                         gATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAgt                             <  1:1331621/149‑1 (MQ=255)
                                                                                                                           tGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGt                                                                                                         >  4:3527155/1‑71 (MQ=255)
                                                                                                                              cGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCa                        <  4:6222413/149‑1 (MQ=255)
                                                                                                                              cGAAATTGAAAGCCATCTGATAGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCa                        <  4:6722250/149‑1 (MQ=255)
                                                                                                                               gAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCt                                                        >  1:7132836/1‑116 (MQ=255)
                                                                                                                               gAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCAc                       >  3:1175166/1‑149 (MQ=255)
                                                                                                                               gAAATTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCAc                       >  4:5031274/1‑149 (MQ=255)
                                                                                                                               gAAATTGAAAGCCATCAGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCAc                       >  2:1838179/1‑149 (MQ=255)
                                                                                                                                aaaTTGAAAGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACt                      >  1:5335250/1‑149 (MQ=255)
                                                                                                                                 aaTTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTg                     >  1:4864307/1‑149 (MQ=255)
                                                                                                                                  aTTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTgc                    >  1:1155216/1‑149 (MQ=255)
                                                                                                                                  aTTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTgc                    <  4:3242749/149‑1 (MQ=255)
                                                                                                                                     gAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGcc                 >  1:3569470/1‑149 (MQ=255)
                                                                                                                                        aGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAgg              >  2:5488066/1‑149 (MQ=255)
                                                                                                                                        aGCCATCTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAgg              >  2:3786798/1‑149 (MQ=255)
                                                                                                                                            aTCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCaaa          >  3:5060073/1‑149 (MQ=255)
                                                                                                                                              cTGATGAGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAAcc        <  4:5575515/149‑1 (MQ=255)
                                                                                                                                                 aTGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGAc     >  1:4336701/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGT  >  NC_000913/93930‑94226

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