Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,538,381 A→G 8.1% E153E (GAA→GAG PP_4871 → lysine decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,538,3810AG 8.1% 198.0 / 6.8 74E153E (GAA→GAGPP_4871lysine decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base A (40/28);  new base G (4/2);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

GCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGAC  >  NC_002947/5538235‑5538528
                                                                                                                                                  |                                                                                                                                                   
gccttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCgaa                                                                                                                                                     >  8:31749/1‑147 (MQ=255)
gccttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCgaa                                                                                                                                                     <  7:31749/147‑1 (MQ=255)
 ccttgccGGGTTGGCTGGGTACGCTGGAGGAACTGTTCGAATTGTTGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACTGTTTCTACGACAAGCTGGGCGTGTTCCTCGACCATATCGTCgaagaa                                                                                                                                                  <  7:310477/149‑1 (MQ=255)
   ttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAgg                                                                                                                                                >  2:88150/1‑149 (MQ=255)
        gggTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCg                                                                                                                                           >  7:22559/1‑149 (MQ=255)
           tGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTgc                                                                                                                                        <  1:88150/149‑1 (MQ=255)
             ggCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTgcgc                                                                                                                                      >  1:314707/1‑149 (MQ=255)
              gCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGcc                                                                                                                                     >  1:326866/1‑149 (MQ=255)
                         tGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTc                                                                                                                                            <  5:62590/131‑1 (MQ=255)
                         tGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTc                                                                                                                                            >  6:62590/1‑131 (MQ=255)
                           gaggaACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCggg                                                                                                                         <  8:22559/148‑1 (MQ=255)
                               aaCTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGAt                                                                                                                     <  2:326866/148‑1 (MQ=255)
                                aCTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATAGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGt                                                                                                                   >  3:113032/1‑149 (MQ=255)
                                    ttCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTtgct                                                                                                               >  3:172190/1‑149 (MQ=255)
                                       gAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTgcgc                                                                                                                                      >  1:5493/1‑123 (MQ=255)
                                       gAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTgcgc                                                                                                                                      <  2:5493/123‑1 (MQ=255)
                                           tgtgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTgggg                                                                                                        >  3:156293/1‑149 (MQ=255)
                                            gtgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGgc                                                                                                       <  3:355850/149‑1 (MQ=255)
                                             tgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTgggg                                                                                                        >  3:264838/1‑147 (MQ=255)
                                                    gggggCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACTCCGGGCGATGTTGCTGCTGGGgcagatcgg                                                                                                >  7:174447/1‑143 (MQ=255)
                                                    gggggCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTGTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTgggg                                                                                                        <  8:174447/140‑1 (MQ=255)
                                                     ggggCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCg                                                                                               >  8:60160/1‑148 (MQ=255)
                                                          aGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAg                                                                                         >  3:309740/1‑149 (MQ=255)
                                                           gCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAgc                                                                                        <  4:156293/149‑1 (MQ=255)
                                                                gCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGctgctg                                                                                   >  8:248893/1‑149 (MQ=255)
                                                                 cTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTgg                                                                                  <  7:60160/149‑1 (MQ=255)
                                                                 cTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTgg                                                                                  <  3:254860/149‑1 (MQ=255)
                                                                 cTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGATgg                                                                                  >  1:315660/1‑149 (MQ=255)
                                                                  tATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTgg                                                                                  >  7:185550/1‑148 (MQ=255)
                                                                   aTCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTgg                                                                                  >  7:25299/1‑147 (MQ=255)
                                                                        gCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGAt                                                                           >  2:25786/1‑149 (MQ=255)
                                                                         ccAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATg                                                                          <  8:25299/149‑1 (MQ=255)
                                                                               cGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGc                                                                    >  4:290759/1‑149 (MQ=255)
                                                                                gCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCt                                                                   >  7:145987/1‑149 (MQ=255)
                                                                                    gggCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCTACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTtgt                                                               <  6:374567/149‑1 (MQ=255)
                                                                                    gggCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCttt                                                                 <  4:172190/147‑1 (MQ=255)
                                                                                    gggCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTtgt                                                               >  7:35442/1‑149 (MQ=255)
                                                                                    gggCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTtgt                                                               <  7:248893/149‑1 (MQ=255)
                                                                                     ggCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTtgtg                                                              <  4:80467/149‑1 (MQ=255)
                                                                                      gCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTgg                                                             <  7:291410/149‑1 (MQ=255)
                                                                                        tACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGgcg                                                           >  6:10338/1‑149 (MQ=255)
                                                                                           tGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCg                                                        >  8:246373/1‑149 (MQ=255)
                                                                                             gATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCg                                                        <  2:314707/147‑1 (MQ=255)
                                                                                             gATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATAGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGAGGGACAGCTTTGTGGCGCCGGt                                                      >  5:344121/1‑149 (MQ=255)
                                                                                                gtgAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGa                                                                        >  4:88710/1‑128 (MQ=255)
                                                                                                gtgAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGa                                                                        <  3:88710/128‑1 (MQ=255)
                                                                                                 tgAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTcgc                                                   <  4:113032/148‑1 (MQ=255)
                                                                                                  gAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGa                                                                        >  8:153593/1‑126 (MQ=255)
                                                                                                  gAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGa                                                                        <  7:153593/126‑1 (MQ=255)
                                                                                                  gAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCgcgc                                                 >  4:32195/1‑149 (MQ=255)
                                                                                                    aCGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGAt                                                                           >  7:277417/1‑121 (MQ=255)
                                                                                                    aCGGTTTCTACGACAAGCTGGGCGGGGTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGAt                                                                           <  8:277417/121‑1 (MQ=255)
                                                                                                        tttCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGCCGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGt                                           >  8:49598/1‑149 (MQ=255)
                                                                                                        tttCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGGCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGt                                           <  8:22272/149‑1 (MQ=255)
                                                                                                            tACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGt                                       >  3:360312/1‑149 (MQ=255)
                                                                                                            tACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGt                                       >  3:254952/1‑149 (MQ=255)
                                                                                                            tACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCCGGCGATGTTGCAGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGCTCGCGCCTAAGTGGGt                                       >  7:277190/1‑149 (MQ=255)
                                                                                                               gACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGa                                    >  7:60698/1‑149 (MQ=255)
                                                                                                                aCAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGAc                                   <  3:82091/149‑1 (MQ=255)
                                                                                                                 cAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGAc                                   >  8:315912/1‑148 (MQ=255)
                                                                                                                 cAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGAc                                   >  2:355901/1‑148 (MQ=255)
                                                                                                                  aaGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACaa                                 <  7:246373/149‑1 (MQ=255)
                                                                                                                   aGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAag                                >  1:13133/1‑149 (MQ=255)
                                                                                                                    gCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAag                                >  1:158692/1‑148 (MQ=255)
                                                                                                                     cTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAagca                              <  8:145987/149‑1 (MQ=255)
                                                                                                                       gggCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAagcagc                            >  8:200641/1‑149 (MQ=255)
                                                                                                                         gCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCct                          <  8:185550/149‑1 (MQ=255)
                                                                                                                           accttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctga                        <  5:300021/146‑1 (MQ=255)
                                                                                                                            ccttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGGAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctgac                       >  4:353388/3‑149 (MQ=255)
                                                                                                                             cttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctgact                      >  4:263597/2‑149 (MQ=255)
                                                                                                                             cttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctgact                      >  3:81521/2‑149 (MQ=255)
                                                                                                                              ttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctgact                      >  4:289132/1‑148 (MQ=255)
                                                                                                                                ccTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctgactga                    <  3:32195/148‑1 (MQ=255)
                                                                                                                                    gACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTg               <  3:290759/149‑1 (MQ=255)
                                                                                                                                      ccATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCt                >  5:254009/1‑146 (MQ=255)
                                                                                                                                             gtcgAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACtttt       <  3:353388/148‑1 (MQ=255)
                                                                                                                                              tcgAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTgg     >  7:115188/1‑149 (MQ=255)
                                                                                                                                                 aagaagGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGAc  <  1:355901/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGAC  >  NC_002947/5538235‑5538528

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: