Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,393,863 G→A 5.7% E798K (GAG→AAG)  PP_2097 → sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,393,8630GA 5.7% 268.1 / 3.4 88E798K (GAG→AAG) PP_2097sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (41/42);  new base A (3/2);  total (44/44)
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.35e-01

CAATGTACCGGTGGAAGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTGCCAGCCAGCTTC  >  NC_002947/2393715‑2394010
                                                                                                                                                    |                                                                                                                                                    
cAATGTACCGGTGGAAGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAggggg                                                                                                                                                      <  2:46640/149‑1 (MQ=255)
cAATGTACCGGGGGAAGGGCGGGGCCAGGGGACTGGAAGGTGCGGCGGGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAggggg                                                                                                                                                      <  5:186606/149‑1 (MQ=255)
         ggTGGAAGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGtat                                                                                                                                             <  5:69845/149‑1 (MQ=255)
         ggTGGAAGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGta                                                                                                                                              >  5:312698/1‑148 (MQ=255)
           tGGAAGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGtata                                                                                                                                            >  5:160973/1‑148 (MQ=255)
             gAAGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATaaaa                                                                                                                                         <  3:85032/149‑1 (MQ=255)
              aaGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAAc                                                                                                                                        <  4:217193/149‑1 (MQ=255)
                 gCGTGGGGCAGGTGACTGCAAGCTGCGGGGTGGCGGGGTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAACCAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAAcc                                                                                                                                       <  8:222453/147‑1 (MQ=255)
                 gCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGt                                                                                                                                     >  5:160346/1‑149 (MQ=255)
                   gTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTgc                                                                                                                                   <  2:259367/149‑1 (MQ=255)
                       gCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTgc‑gcgc                                                                                                                              >  7:31948/1‑149 (MQ=255)
                       gCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCCGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTgc‑gcgc                                                                                                                              >  6:122349/1‑149 (MQ=255)
                            gTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCg                                                                                                                         >  4:183502/1‑149 (MQ=255)
                            gTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCg                                                                                                                         >  8:191706/1‑149 (MQ=255)
                               aCTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCg                                                                                                                      >  2:269855/1‑149 (MQ=255)
                                cTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGt                                                                                                                     >  5:163597/1‑149 (MQ=255)
                                 tGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGCGGTCCGGTCGGCGt                                                                                                                     >  1:277865/1‑149 (MQ=255)
                                 tGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTg                                                                                                                    >  3:7810/1‑149 (MQ=255)
                                 tGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTg                                                                                                                    >  8:94527/1‑149 (MQ=255)
                                 tGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGAAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTg                                                                                                                    >  7:110593/1‑149 (MQ=255)
                                   cAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTc                                                                                                                          <  5:154713/141‑1 (MQ=255)
                                   cAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTc                                                                                                                          >  6:154713/1‑141 (MQ=255)
                                   cAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTgg                                                                                                                   <  7:29437/148‑1 (MQ=255)
                                    aaGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTgggg                                                                                                                 >  6:84381/1‑149 (MQ=255)
                                     aGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTggggg                                                                                                                <  2:154764/149‑1 (MQ=255)
                                      gCTTCGGTGTGGCGGGCTGGCGGTCGTGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGc                                                                                                               <  7:140968/149‑1 (MQ=255)
                                      gCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGc                                                                                                               >  2:161078/1‑149 (MQ=255)
                                           ggTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCtt                                                                                                          <  5:84381/149‑1 (MQ=255)
                                                 gCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCg                                                                                                      >  6:80392/1‑147 (MQ=255)
                                                 gCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCg                                                                                                      <  5:80392/147‑1 (MQ=255)
                                                 gCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGc                                                                                                    <  8:53511/149‑1 (MQ=255)
                                                         gCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTg                                                                                                                    <  7:288390/125‑1 (MQ=255)
                                                         gCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTc                                                                                              <  3:276856/147‑1 (MQ=255)
                                                         gCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTc                                                                                              <  8:31948/147‑1 (MQ=255)
                                                         gCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTc                                                                                              <  1:96704/147‑1 (MQ=255)
                                                         gCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGAGGTCGGCGTg                                                                                                                    >  8:288390/1‑125 (MQ=255)
                                                           ggTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCGTTCGGGa                                                                                          >  7:261657/1‑149 (MQ=255)
                                                               gggCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTc                                                                                              <  5:72100/141‑1 (MQ=255)
                                                               gggCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTc                                                                                              >  6:72100/1‑141 (MQ=255)
                                                                 gCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGa                                                                                          <  4:124807/143‑1 (MQ=255)
                                                                 gCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGa                                                                                          >  3:124807/1‑143 (MQ=255)
                                                                    aaGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAAGAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGcc                                                                                 <  1:269855/149‑1 (MQ=255)
                                                                    aaGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGcc                                                                                 >  8:3460/1‑149 (MQ=255)
                                                                       gcgcCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCg                                                                                             >  1:150817/1‑134 (MQ=255)
                                                                       gcgcCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCg                                                                                             <  2:150817/134‑1 (MQ=255)
                                                                       gcgcCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCt                                                                              >  6:40228/1‑149 (MQ=255)
                                                                       gcgcCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCt                                                                              >  6:19553/1‑149 (MQ=255)
                                                                           cGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTg                                                                          >  7:249574/1‑149 (MQ=255)
                                                                           cGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTg                                                                          <  6:160973/149‑1 (MQ=255)
                                                                                 cTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGAAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCg                                                                                             <  3:234438/124‑1 (MQ=255)
                                                                                 cTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGAAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCg                                                                                             <  3:238719/124‑1 (MQ=255)
                                                                                 cTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGAAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCg                                                                                             >  4:238719/1‑124 (MQ=255)
                                                                                 cTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGAAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCg                                                                                             >  4:234438/1‑124 (MQ=255)
                                                                                    gCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGAcc                                                                  <  3:7817/148‑1 (MQ=255)
                                                                                            gcgcCGACGCCGGCGTGTATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCg                                                         >  8:95745/1‑149 (MQ=255)
                                                                                                cGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGcc                                                     >  6:201001/1‑149 (MQ=255)
                                                                                                  aCGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCg                                                   <  8:249574/149‑1 (MQ=255)
                                                                                                    gCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCggg                                                 >  4:186897/1‑149 (MQ=255)
                                                                                                     ccGGCGTATATGCGGCCAAACCGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGGGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGgt                                                <  3:193448/149‑1 (MQ=255)
                                                                                                         cGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGgtgtgt                                            <  5:40228/149‑1 (MQ=255)
                                                                                                           tataTGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTgg                                          <  7:94527/149‑1 (MQ=255)
                                                                                                            atatGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTgg                                          >  4:42078/1‑148 (MQ=255)
                                                                                                            atatGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTgg                                          <  8:261657/148‑1 (MQ=255)
                                                                                                             tatGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCa                                        <  3:183502/149‑1 (MQ=255)
                                                                                                             tatGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCAACGACCGCGCAAGCGGGCCCGGGTGTGTGGCa                                        <  7:191706/149‑1 (MQ=255)
                                                                                                              atGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAg                                       >  5:137528/1‑149 (MQ=255)
                                                                                                               tGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAgg                                      <  6:163597/149‑1 (MQ=255)
                                                                                                                gCGGCCAAACAGGCGGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGc                                                          <  8:276114/128‑1 (MQ=255)
                                                                                                                gCGGCCAAACAGGCGGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGc                                                          >  7:276114/1‑128 (MQ=255)
                                                                                                                 cGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGt                                     <  5:122349/148‑1 (MQ=255)
                                                                                                                    ccAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCg                                 >  8:6459/1‑149 (MQ=255)
                                                                                                                     cAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGt                                <  5:19553/149‑1 (MQ=255)
                                                                                                                          aGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGc                                                          >  2:308906/1‑118 (MQ=255)
                                                                                                                          aGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGc                                                          <  1:308906/118‑1 (MQ=255)
                                                                                                                            gCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTgg                         >  2:283561/1‑149 (MQ=255)
                                                                                                                            gCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTgg                         >  2:14587/1‑149 (MQ=255)
                                                                                                                             ccGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGc                        <  3:42078/149‑1 (MQ=255)
                                                                                                                              cgggcgggACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGcc                       <  1:14587/149‑1 (MQ=255)
                                                                                                                              cgggcgggACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATAGACCGCGCAAGCGGGCCCgggag                                               >  4:53929/1‑123 (MQ=255)
                                                                                                                                    ggACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTa                 <  2:277865/149‑1 (MQ=255)
                                                                                                                                    ggACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTa                 <  4:7810/149‑1 (MQ=255)
                                                                                                                                      aCCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACt               <  7:95745/149‑1 (MQ=255)
                                                                                                                                        cGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAACGCCGCTCCTGCATCGACCGCGCACGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTgc             >  5:158473/1‑149 (MQ=255)
                                                                                                                                        cGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTa                 <  1:7993/145‑1 (MQ=255)
                                                                                                                                        cGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTa                 >  2:7993/1‑145 (MQ=255)
                                                                                                                                         gggTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGACGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTgc             >  7:206667/1‑148 (MQ=255)
                                                                                                                                          ggTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCggg                                                       >  4:6271/1‑105 (MQ=255)
                                                                                                                                          ggTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCggg                                                       <  3:6271/105‑1 (MQ=255)
                                                                                                                                            gggAGGGGGAGCCGGTATAAACCCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTgccagc         <  3:5303/148‑1 (MQ=255)
                                                                                                                                              gAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTgccagc         >  1:38795/1‑147 (MQ=255)
                                                                                                                                                 ggggAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTGCCAGCCAGCt    >  8:94126/1‑149 (MQ=255)
                                                                                                                                                   ggAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTGCCAGCCAGCTTc  >  4:108118/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAATGTACCGGTGGAAGGCGTGGGCCAGGTGACTGCAAGCTTCGGTGTGGCGGGCTGGCGGTCGGGCGAAGGCGCCGATGCCTTGCTGTTGCGCGCCGACGCCGGCGTATATGCGGCCAAACAGGCCGGGCGGGACCGGGTGGAGGGGGAGCCGGTATAAAACCGTGC‑GCGCGGTCGGCGTGGGGGCCGCTTTGCGGCCCTTTCGGGACGCAAAGCCGCTCCTGCATCGACCGCGCAAGCGGGCCCGGGTGTGTGGCAGGTGGCGTCAGTTGGCCGCTGTACTGCCAGCCAGCTTC  >  NC_002947/2393715‑2394010

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