Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 508,774 Δ1 bp 7.2% coding (240/1890 nt) estP ← esterase EstP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947508,7740C. 7.2% 218.9 / 11.0 69coding (240/1890 nt)estPesterase EstP
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base . (3/2);  total (39/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 = 9.64e-01

AACCCGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGA  >  NC_002947/508625‑508919
                                                                                                                                                     |                                                                                                                                                  
aaCCCGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGgc                                                                                                                                                    >  2:955132/1‑150 (MQ=255)
aaCCCGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGgc                                                                                                                                                    <  2:116409/150‑1 (MQ=255)
              gTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggcg                                                                                                                                     <  2:1099282/150‑1 (MQ=255)
               tttCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggcgg                                                                                                                                    >  1:1197241/1‑150 (MQ=255)
                ttCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggcggc                                                                                                                                   <  1:1157459/150‑1 (MQ=255)
                   cTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCAc                                                                                                                                >  1:889258/1‑150 (MQ=255)
                   cTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCAc                                                                                                                                >  2:541203/1‑150 (MQ=255)
                      gggggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCCGCTGCACATCCTTTGTGGACCCGCCGAGTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAGGCGAccttaac                                                                                                                             <  1:870956/150‑6 (MQ=255)
                       ggggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAg                                                                                                                            <  1:955132/150‑1 (MQ=255)
                         ggATGAAAGTCTCAGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGAGGTGTACCCGCCGACTGCCCAGTTATTAACATCGGGCTGGATACCTACGGATGTCGACGGGTTCAGATCGCCTGGCGCGACGC‑CAAGACGGCAACCGAGGa                                                                                                                          >  2:1152157/1‑150 (MQ=255)
                          gATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGCCGCCTGG‑GCGTCGCGCCAggcgac                                                                                                                                   >  1:1003283/1‑138 (MQ=255)
                          gATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGCCGCCTGG‑GCGTCGCGCCAggcgac                                                                                                                                   <  2:1003283/140‑3 (MQ=255)
                            tGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGatca                                                                                                                       >  2:502883/1‑150 (MQ=255)
                                       gACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGGGGTGTACCCGCCGACTGCCCAGGTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACAcc                                                                                                            >  1:696755/1‑150 (MQ=255)
                                        aCGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCg                                                                                                           >  1:830210/1‑150 (MQ=255)
                                              tCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGa                                                                                                     >  2:857249/1‑150 (MQ=255)
                                               cGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGaa                                                                                                    <  1:291943/150‑1 (MQ=255)
                                                   aTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGtt                                                                                                 <  1:660405/149‑1 (MQ=255)
                                                   aTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGtt                                                                                                 >  2:660405/1‑149 (MQ=255)
                                                       aGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGt                                                                                            >  2:121233/1‑150 (MQ=255)
                                                         tCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCg                                                                                          >  1:387303/1‑150 (MQ=255)
                                                           cAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGc                                                                                        >  2:815876/1‑150 (MQ=255)
                                                               aTCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCa                                                                                    <  2:387303/150‑1 (MQ=255)
                                                                tCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGcc                                                                                              <  1:62058/139‑1 (MQ=255)
                                                                tCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGcc                                                                                              >  2:62058/1‑139 (MQ=255)
                                                                tCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCAc                                                                                   >  2:1001314/1‑150 (MQ=255)
                                                                    ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGATGAGCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGtt                                                                               <  1:23620/150‑1 (MQ=255)
                                                                     tgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCAcc                                                                                                                               >  2:784197/1‑101 (MQ=255)
                                                                     tgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCAcc                                                                                                                               <  1:784197/101‑1 (MQ=255)
                                                                     tgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTg                                                                              <  2:78339/150‑1 (MQ=255)
                                                                       ggtggtGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                            <  2:101159/150‑1 (MQ=255)
                                                                         tggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGa                                                                          >  1:164422/1‑150 (MQ=255)
                                                                          ggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGaa                                                                         <  1:857249/150‑1 (MQ=255)
                                                                                 ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgcg                                                                  >  2:1194958/1‑150 (MQ=255)
                                                                                 ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgcg                                                                  >  1:963037/1‑150 (MQ=255)
                                                                                    ccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCg                                                                                       >  1:961567/1‑126 (MQ=255)
                                                                                    ccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCg                                                                                       <  2:961567/126‑1 (MQ=255)
                                                                                         tGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTgg                                                                             <  2:129678/131‑1 (MQ=255)
                                                                                         tGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTgg                                                                             >  1:129678/1‑131 (MQ=255)
                                                                                            ccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACAGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAcc                                                       >  1:820399/1‑150 (MQ=255)
                                                                                                  attaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaac                                                 <  2:963037/150‑1 (MQ=255)
                                                                                                       cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGAt                                            >  2:1051678/1‑150 (MQ=255)
                                                                                                       cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGAt                                            >  2:1060701/1‑150 (MQ=255)
                                                                                                        aTCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                           >  2:290667/1‑150 (MQ=255)
                                                                                                          cGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGATGAAACGCGCGCCTGGGGTACCGCCAACAAGATCgg                                         <  1:1001314/150‑1 (MQ=255)
                                                                                                            ggCTTGATACCTACGCATGTCGACGGTTTAAGGTCGCCTGGCGCGACGC‑CCAGCCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTCGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGa                                       <  1:1152157/150‑1 (MQ=255)
                                                                                                              cTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTggg                                                            >  2:547982/1‑127 (MQ=255)
                                                                                                              cTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTggg                                                            <  1:547982/127‑1 (MQ=255)
                                                                                                               tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACgg                                                                                 >  2:860459/1‑105 (MQ=255)
                                                                                                               tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACgg                                                                                 <  1:860459/105‑1 (MQ=255)
                                                                                                                       tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGa                                             >  2:931259/1‑133 (MQ=255)
                                                                                                                       tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGa                                             <  1:931259/133‑1 (MQ=255)
                                                                                                                        aCGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGcc                                >  1:1074040/1‑145 (MQ=255)
                                                                                                                        aCGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGcc                                <  2:1074040/145‑1 (MQ=255)
                                                                                                                            aTGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCgg                       >  1:1161382/1‑150 (MQ=255)
                                                                                                                                 gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                           <  2:803691/125‑1 (MQ=255)
                                                                                                                                 gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                           >  1:803691/1‑125 (MQ=255)
                                                                                                                                  aCGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGc                 >  2:142070/1‑150 (MQ=255)
                                                                                                                                  aCGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGc                 <  1:502883/150‑1 (MQ=255)
                                                                                                                                      gTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTc             >  1:154786/1‑150 (MQ=255)
                                                                                                                                      gTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTc             <  2:889258/150‑1 (MQ=255)
                                                                                                                                          aGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                  <  1:982788/109‑1 (MQ=255)
                                                                                                                                          aGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                  >  2:982788/1‑109 (MQ=255)
                                                                                                                                           ggTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGGAACTGCCCGGCATCGGCGAGGCTGTCGCCaa        >  2:190641/1‑150 (MQ=255)
                                                                                                                                            gTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCg                                          >  2:406574/1‑115 (MQ=255)
                                                                                                                                            gTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCg                                          <  1:406574/115‑1 (MQ=255)
                                                                                                                                              cGCCTGG‑GCGTCGCGCCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAaca     >  2:1153907/1‑150 (MQ=255)
                                                                                                                                              cGCCTGG‑GCGTCGCGCCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAaca     >  1:735917/1‑150 (MQ=255)
                                                                                                                                                 cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAg                   <  2:630882/133‑1 (MQ=255)
                                                                                                                                                 cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAg                   >  1:630882/1‑133 (MQ=255)
                                                                                                                                                 cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGa  <  1:290667/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AACCCGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGA  >  NC_002947/508625‑508919

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