Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,919 A→T 5.3% E48V (GAA→GTA)  mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947635,9190AT 5.3% 230.8 / 2.8 75E48V (GAA→GTA) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base T (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGCCCCCAGTGGCTGGCCGA  >  NC_002947/635773‑636066
                                                                                                                                                   |                                                                                                                                                   
gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGaa                                                                                                                                                    <  1:140715/148‑1 (MQ=255)
   cATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGa                                                                                                                                                <  6:173984/149‑1 (MQ=255)
   cATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGa                                                                                                                                                <  5:100457/149‑1 (MQ=255)
   cATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGATGGCGGTGCTGTCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGa                                                                                                                                                <  8:323274/149‑1 (MQ=255)
    aTGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGAc                                                                                                                                               >  6:1895/1‑149 (MQ=255)
      gCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGATATCGAACTGAcc                                                                                                                                              >  3:227540/1‑148 (MQ=255)
           ttCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAgg                                                                                                                                         <  4:231802/148‑1 (MQ=255)
                          tttGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAg                                                                                                                         >  5:245512/1‑149 (MQ=255)
                          tttGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAg                                                                                                                         >  5:317357/1‑149 (MQ=255)
                              cGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTgg                                                                                                                     <  1:55972/149‑1 (MQ=255)
                               ggCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGa                                                                                                                    <  6:127791/149‑1 (MQ=255)
                               ggCACTTTCATGGGTTCGCTGGCGGGGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGa                                                                                                                    <  8:201288/149‑1 (MQ=255)
                                   cTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCg                                                                                                                >  2:360351/1‑149 (MQ=255)
                                    tttCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCgg                                                                                                               >  2:170893/1‑149 (MQ=255)
                                    tttCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCgg                                                                                                               >  2:320669/1‑149 (MQ=255)
                                        aTGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTAt                                                                                                                                     <  8:257702/123‑1 (MQ=255)
                                        aTGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTAt                                                                                                                                     >  7:257702/1‑123 (MQ=255)
                                         tGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                                                                                                            >  2:120635/1‑147 (MQ=255)
                                         tGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCAGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCa                                                                                                          <  6:245512/149‑1 (MQ=255)
                                          gggTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                                                                                                            >  7:77198/1‑146 (MQ=255)
                                             ttCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGCACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGAc                                                                                                       >  1:88606/1‑149 (MQ=255)
                                             ttCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGAcc                                                                                                      >  7:310904/1‑149 (MQ=255)
                                                 cTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCCGCTGGATCCGGCGCCAGACCtggt                                                                                                  >  8:322329/1‑149 (MQ=255)
                                                  tGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                                                                                                 >  3:115152/1‑149 (MQ=255)
                                                    gCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtgg                                                                                                <  7:106930/148‑1 (MQ=255)
                                                      ggTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                                                                                                  <  6:177265/128‑1 (MQ=255)
                                                      ggTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                                                                                                  >  5:177265/1‑128 (MQ=255)
                                                          cTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg                                                                                         <  6:317357/149‑1 (MQ=255)
                                                          cTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCg                                                                                          <  3:291570/148‑1 (MQ=255)
                                                                  ggAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCa                                                                                 >  2:3232/1‑149 (MQ=255)
                                                                  ggAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCa                                                                                 <  5:104317/149‑1 (MQ=255)
                                                                     aCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGt                                                                              <  5:1895/149‑1 (MQ=255)
                                                                            cACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCgg                                                                       >  2:292650/1‑149 (MQ=255)
                                                                               cgcgTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtgg                                                                                                   >  5:37503/1‑118 (MQ=255)
                                                                               cgcgTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtgg                                                                                                   <  6:37503/118‑1 (MQ=255)
                                                                               cgcgGCACCGGCTCTGACTCCAACGTCTATACCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGc                                                                      <  6:259890/147‑1 (MQ=255)
                                                                                   tCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCa                                                                <  1:3232/149‑1 (MQ=255)
                                                                                       cGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCgg                                                            >  2:318714/1‑149 (MQ=255)
                                                                                         gCTCTGACGGCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGtt                                                          <  2:88606/149‑1 (MQ=255)
                                                                                         gCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGtt                                                          <  8:281953/149‑1 (MQ=255)
                                                                                         gCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGtt                                                          <  1:170893/149‑1 (MQ=255)
                                                                                         gCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGtt                                                          <  4:371482/149‑1 (MQ=255)
                                                                                                aCGCCAACGTTTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAAt                                                      <  1:320669/147‑1 (MQ=255)
                                                                                                aCGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGATATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAAt                                                      <  4:227540/147‑1 (MQ=255)
                                                                                                  gCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGt                                                  >  1:332677/1‑149 (MQ=255)
                                                                                                  gCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGt                                                  >  7:111282/1‑149 (MQ=255)
                                                                                                    cAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGc                                                >  5:162762/1‑149 (MQ=255)
                                                                                                    cAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGc                                                >  7:202369/1‑149 (MQ=255)
                                                                                                    cAACGTCTATCCCCAGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGc                                                >  7:319235/1‑149 (MQ=255)
                                                                                                     aaCGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCt                                               >  7:135636/1‑149 (MQ=255)
                                                                                                         tCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTgca                                            >  7:258614/1‑146 (MQ=255)
                                                                                                              cccccGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAggg                                      >  2:32734/1‑149 (MQ=255)
                                                                                                              cccccGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAggg                                      >  1:189522/1‑149 (MQ=255)
                                                                                                               ccccGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGTCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGt                                                  <  3:280504/136‑1 (MQ=255)
                                                                                                               ccccGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCgg                                                            >  7:117281/1‑126 (MQ=255)
                                                                                                               ccccGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCgg                                                            <  8:117281/126‑1 (MQ=255)
                                                                                                               ccccGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGt                                                  >  4:280504/1‑136 (MQ=255)
                                                                                                               ccccGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGc                                     >  7:154032/1‑149 (MQ=255)
                                                                                                                    ctccccACCCAGCTCCAAGCCCAGGCTATCGAACTCACCCAAGCCGATCACCCCGCCCAGCTGGATCCGGCTCCAGACCTTGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGTCCTGcc                                <  6:288562/144‑1 (MQ=255)
                                                                                                                        gCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTAc                            <  1:360351/149‑1 (MQ=255)
                                                                                                                            ccAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTg                                  >  7:47861/1‑139 (MQ=255)
                                                                                                                            ccAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTg                                  <  8:47861/139‑1 (MQ=255)
                                                                                                                               gCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCg                     <  4:115152/149‑1 (MQ=255)
                                                                                                                               gCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGGAGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCg                     <  1:292650/149‑1 (MQ=255)
                                                                                                                                 tCGAAGCCCAGGGTATCGAACTGACACAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCCACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGc                   >  7:312759/1‑149 (MQ=255)
                                                                                                                                   gAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGccc                 >  4:188459/1‑149 (MQ=255)
                                                                                                                                        ccAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAg                                                               >  4:36995/1‑98 (MQ=255)
                                                                                                                                        ccAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAg                                                               <  3:36995/98‑1 (MQ=255)
                                                                                                                                          aGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGCCCCCAGtgg          >  4:110254/1‑149 (MQ=255)
                                                                                                                                                tCGTACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGa                                             <  8:60186/108‑1 (MQ=255)
                                                                                                                                                tCGTACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGa                                             <  8:65384/108‑1 (MQ=255)
                                                                                                                                                tCGTACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGa                                             >  7:65384/1‑108 (MQ=255)
                                                                                                                                                tCGTACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGa                                             >  7:60186/1‑108 (MQ=255)
                                                                                                                                                tCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGCCCCCAGTGGCTGGcc    <  1:120635/149‑1 (MQ=255)
                                                                                                                                                  gAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGCCCCCAGTGGCTGGCCGa  >  6:380604/1‑149 (MQ=255)
                                                                                                                                                  gAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGCCCCCAGTGGCTGGCCGa  >  1:219203/1‑149 (MQ=255)
                                                                                                                                                   aaCTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCg                     >  8:238354/1‑129 (MQ=255)
                                                                                                                                                   aaCTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCg                     <  7:238354/129‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTCGCGCGGCAACCCAGCGGTTGAATACGTGCTGAACAAGGGCCTGCCCTACGTCTCCGGCCCCCAGTGGCTGGCCGA  >  NC_002947/635773‑636066

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