Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,307 C→G 5.0% P316R (CCG→CGG)  PP_3812 → lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,3070CG 5.0% 256.0 / 2.7 79P316R (CCG→CGG) PP_3812lipase
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base G (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GCTTGCCACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAA  >  NC_002947/4338160‑4338454
                                                                                                                                                   |                                                                                                                                                   
gCTTGCCACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCg                                                                                                                                                    >  3:61268/1‑149 (MQ=255)
    gCCACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGtt                                                                                                                                                  <  1:208016/147‑1 (MQ=255)
    gCCACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGtt                                                                                                                                                  <  8:254376/147‑1 (MQ=255)
    gCCACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTg                                                                                                                                                 <  5:249161/148‑1 (MQ=255)
     ccACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAg                                                                                                                                               >  1:157222/1‑149 (MQ=255)
      cACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGc                                                                                                                                              >  4:50637/1‑149 (MQ=255)
         cGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACGGGTTGAGcagc                                                                                                                                           <  2:134527/149‑4 (MQ=255)
         cGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGcc                                                                                                                                             >  1:24046/1‑147 (MQ=255)
           gCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGt                                                                                                                                         <  8:179603/149‑1 (MQ=255)
                 gTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCgg                                                                                                                                    <  5:125378/148‑1 (MQ=255)
                   gCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAg                                                                                                                                 >  3:242249/1‑149 (MQ=255)
                       aCGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTaaa                                                                                                                             <  7:193325/149‑1 (MQ=255)
                       aCGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTa                                                                                                                               >  5:374656/1‑147 (MQ=255)
                         ggCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCACGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAgg                                                                                                                           <  4:345920/149‑1 (MQ=255)
                         ggCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAgg                                                                                                                           <  1:193134/149‑1 (MQ=255)
                                 tCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTg                                                                                                                   >  1:164377/1‑149 (MQ=255)
                                   gATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGAGTGACAGTAACCTCGGTGCGCTACAACGGAAAGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCg                                                                                                                 >  7:151787/1‑149 (MQ=255)
                                             gggATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCggg                                                                                                        >  1:7708/1‑148 (MQ=255)
                                             gggATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCggg                                                                                                        <  3:166422/148‑1 (MQ=255)
                                              ggATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCTGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTaaa                                                                                                                             <  5:379244/126‑1 (MQ=255)
                                              ggATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTaaa                                                                                                                             >  6:379244/1‑126 (MQ=255)
                                              ggATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGc                                                                                                       >  8:304762/1‑148 (MQ=255)
                                                aTGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAAc                                                                                                    <  6:374656/149‑1 (MQ=255)
                                                   aaGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGaa                                                                                                 <  1:77193/149‑1 (MQ=255)
                                                    aGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAAc                                                                                                >  7:259354/1‑149 (MQ=255)
                                                    aGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAAc                                                                                                <  3:169732/149‑1 (MQ=255)
                                                      gTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACtt                                                                                              <  2:24046/149‑1 (MQ=255)
                                                       tGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACtt                                                                                              >  4:259414/1‑148 (MQ=255)
                                                       tGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACtt                                                                                              <  3:259414/148‑1 (MQ=255)
                                                        gAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAgg                                                                                                                           >  1:60407/1‑118 (MQ=255)
                                                        gAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAgg                                                                                                                           <  2:60407/118‑1 (MQ=255)
                                                         aGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTaaa                                                                                           <  7:304762/149‑1 (MQ=255)
                                                             gTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGc                                                                                       <  3:50637/149‑1 (MQ=255)
                                                                      cAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGt                                                                              >  7:269455/1‑149 (MQ=255)
                                                                         gTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTAATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAg                                                                           >  1:309352/1‑149 (MQ=255)
                                                                           tGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAgg                                                                          <  2:164377/148‑1 (MQ=255)
                                                                                 cggcggGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGa                                                                    <  4:242249/148‑1 (MQ=255)
                                                                                  ggcgggcgTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGc                                                                  <  4:61268/149‑1 (MQ=255)
                                                                                  ggcgggcgTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGc                                                                  <  4:104102/149‑1 (MQ=255)
                                                                                         gTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCATGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGCGCAATacac                                                           >  6:370577/1‑149 (MQ=255)
                                                                                            aCAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCg                                                        >  6:384228/1‑149 (MQ=255)
                                                                                             cAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGc                                                       >  4:6424/1‑149 (MQ=255)
                                                                                             cAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGc                                                       <  2:309352/149‑1 (MQ=255)
                                                                                              aGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGcc                                                      >  1:379883/1‑149 (MQ=255)
                                                                                                tAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTg                                                    >  8:114582/1‑149 (MQ=255)
                                                                                                  aCCTCGGTGCGCTACAACGGATTGATTCATGACTATGGGCTGCTCAACCCGTTGCGCCAAGTGCCGGCAGTAAAGGCGGCGCTGGGGCAGGCGGGAAACGAAATTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGcagggg                                                  >  6:368756/1‑144 (MQ=255)
                                                                                                  aCCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTgg                                                   <  2:157222/148‑1 (MQ=255)
                                                                                                    cTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCg                                                >  6:105812/1‑149 (MQ=255)
                                                                                                     tCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                               >  4:350096/1‑149 (MQ=255)
                                                                                                       ggTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCggg                                              <  3:6424/148‑1 (MQ=255)
                                                                                                          gcgcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtt                                          >  8:259468/1‑149 (MQ=255)
                                                                                                           ctcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtt                                          <  5:368756/146‑1 (MQ=255)
                                                                                                           cgcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCCAGGCGGGTCtat                                         <  2:7708/149‑3 (MQ=255)
                                                                                                            gcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtttt                                        <  5:384228/149‑1 (MQ=255)
                                                                                                            gcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtttt                                        <  3:350096/149‑1 (MQ=255)
                                                                                                                cAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCtt                                    <  5:105812/149‑1 (MQ=255)
                                                                                                                cAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCt                                     >  1:337993/1‑148 (MQ=255)
                                                                                                                 aaCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGt                                             >  4:157457/1‑139 (MQ=255)
                                                                                                                 aaCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGt                                             <  3:157457/139‑1 (MQ=255)
                                                                                                                    ggAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAt                                >  8:279897/1‑149 (MQ=255)
                                                                                                                        tGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                            <  8:259354/149‑1 (MQ=255)
                                                                                                                          aTTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAg                                                                   >  4:134109/1‑108 (MQ=255)
                                                                                                                          aTTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAg                                                                   <  3:134109/108‑1 (MQ=255)
                                                                                                                          aTTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCgg                          <  7:279897/149‑1 (MQ=255)
                                                                                                                           ttCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGa                         <  7:114582/149‑1 (MQ=255)
                                                                                                                              ctGACTATCGGCTGCTCAACCCGTTGAGCCAAGTGCAGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGaaa                       <  8:269455/147‑1 (MQ=255)
                                                                                                                                gACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAg                    >  1:368001/1‑149 (MQ=255)
                                                                                                                                       ggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACGCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCa              <  7:259468/148‑1 (MQ=255)
                                                                                                                                        ttggctCAACGGGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAgg            >  1:91693/4‑149 (MQ=255)
                                                                                                                                        ttggctCAACGGGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAgg            >  6:22631/4‑149 (MQ=255)
                                                                                                                                           gctCAACGGGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCtt         <  5:22631/149‑1 (MQ=255)
                                                                                                                                            ctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCttt                                         >  2:288871/1‑116 (MQ=255)
                                                                                                                                            ctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCttt                                         <  1:288871/116‑1 (MQ=255)
                                                                                                                                             tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                               <  8:86078/109‑1 (MQ=255)
                                                                                                                                             tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                               >  7:86078/1‑109 (MQ=255)
                                                                                                                                             tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTc        >  8:176411/1‑148 (MQ=255)
                                                                                                                                                aCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGAt                                                                      <  3:99400/83‑1 (MQ=255)
                                                                                                                                                aCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGAt                                                                      >  4:99400/1‑83 (MQ=255)
                                                                                                                                                aCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAAt     >  4:66849/1‑148 (MQ=255)
                                                                                                                                                 cccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATg    >  1:225663/1‑148 (MQ=255)
                                                                                                                                                  ccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGaa  <  3:321804/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTTGCCACCGGCCCTGGTGCAAACGGCTGAGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAA  >  NC_002947/4338160‑4338454

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