Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,285,028 G→A 6.7% G87S (GGT→AGT) ‡ PP_1122 → OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,285,0280GA 6.7% 219.7 / 4.0 75G87S (GGT→AGT) ‡PP_1122OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base A (3/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

GGTGGCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAA  >  NC_002947/1284880‑1285174
                                                                                                                                                    |                                                                                                                                                   
ggtggCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCggtg                                                                                                                                                     >  6:339607/1‑149 (MQ=255)
  tggCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCggtggt                                                                                                                                                   >  2:30729/1‑149 (MQ=255)
      tGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGtggtgggg                                                                                                                                               >  3:185316/1‑147 (MQ=255)
        ccAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGgcgc                                                                                                                                              <  3:47390/148‑1 (MQ=255)
        ccAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGcc                                                                                                                                             >  1:55367/1‑149 (MQ=255)
         cAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTAGCCCATCTGCAGCCTGGTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCt                                                                                                                                           >  1:111056/1‑149 (MQ=255)
           gCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTg                                                                                                                                         >  4:34348/1‑149 (MQ=255)
                cTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGGTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgc                                                                                                                                      <  3:291449/147‑1 (MQ=255)
                 tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGgtggtg                                                                                                                                                  >  1:292400/1‑135 (MQ=255)
                 tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGgtggtg                                                                                                                                                  <  2:292400/135‑1 (MQ=255)
                 tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgcgc                                                                                                                                    >  7:228542/1‑148 (MQ=255)
                   aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAt                                                                                                                                 >  7:210102/1‑149 (MQ=255)
                   aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAt                                                                                                                                 <  3:224222/149‑1 (MQ=255)
                    cGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAtt                                                                                                                                <  3:186344/149‑1 (MQ=255)
                     gggCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCGTTGGTGCCATCGAGAGTTCTTCATGGGCCGCCGGTGGTGGCGCCGTGGGCGCCATTg                                                                                                                               >  8:135020/1‑149 (MQ=255)
                      ggCTGCGCTACGCACAGCTATGGCTGCGCGCCCCTCAATCAAAGGCCCTGGCCCTTCTGCCGCCTTCTGGGGGGATTGGTGGGTGTTGGCCTTGTTGCAATGGATATTTCTTCCTGGGCCGCCTGTGGTGGCGCCCAAGGCGCCaccac                                                                                                                              <  4:185316/149‑5 (MQ=255)
                         tGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGCCCGCCGGCAATGGCGCCCAAGGCGCCaccaccgg                                                                                                                           >  1:130655/1‑142 (MQ=38)
                            gcTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCcggcgg                                                                                                                        <  8:210102/149‑1 (MQ=255)
                             cTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGcc                                                                                                                                   <  6:368326/137‑1 (MQ=255)
                             cTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGcc                                                                                                                                   >  5:368326/1‑137 (MQ=255)
                              tACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGGGGAATATTACATGGCCCCCGGGTGTTGCGCCGTTGGCGGCCATTGCcgccggc                                                                                                                       >  1:308335/1‑148 (MQ=255)
                                cGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTg                                                                                                                    <  2:111056/149‑1 (MQ=255)
                                cGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGCCCGCCGGCAATGGCGCCCAAGGCGCCACCACcggcggccca                                                                                                                    <  2:130655/149‑5 (MQ=37)
                                  cacaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCCCCTTGGGCTCCATTTCCCGCCGGCAGAt                                                                                                                  >  8:305464/1‑149 (MQ=255)
                                   acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAtt                                                                                                                 >  1:5024/1‑149 (MQ=255)
                                   acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAtt                                                                                                                 <  4:27740/149‑1 (MQ=255)
                                   acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAtt                                                                                                                 >  5:26575/1‑149 (MQ=255)
                                   acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGCAATGGCGCCCAAGGCGCCACCACcggcggcccagga                                                                                                                 <  8:297639/149‑8 (MQ=37)
                                   acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGCAATGGCGCCCAAGGCGCCACCACcggcggcccagg                                                                                                                  >  7:297639/1‑142 (MQ=37)
                                    caGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                >  3:104679/1‑149 (MQ=255)
                                        gATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTa                                                                                                            >  3:173619/1‑149 (MQ=255)
                                            gCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGcc                                                                                                        <  4:161034/149‑1 (MQ=255)
                                               gcgcgcCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAg                                                                                                     <  4:6783/149‑1 (MQ=255)
                                                   gcCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTGGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACg                                                                                                 <  3:34348/149‑1 (MQ=255)
                                                    ccccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACgg                                                                                                >  6:317714/1‑149 (MQ=255)
                                                    ccccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACgg                                                                                                <  5:339607/149‑1 (MQ=255)
                                                     cccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGt                                                                                               <  3:55148/149‑1 (MQ=255)
                                                       cTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                             >  8:227144/1‑149 (MQ=255)
                                                       cTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                             >  5:165181/1‑149 (MQ=255)
                                                          aatcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACg                                                                                           >  1:347923/1‑148 (MQ=255)
                                                            tcaaAGGACCTGGCCCATCTGCAGCCTGCTCGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCAGGGGCCCCAATTCCGGCCGGCTGAATTGGTTCGCCCACGCGGGTGACg                                                                                           >  4:698/1‑146 (MQ=255)
                                                               aaGGACCTGGCCCATCTGCAGGCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                     <  4:104679/149‑1 (MQ=255)
                                                               aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                     <  7:316007/149‑1 (MQ=255)
                                                               aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                     <  3:156606/149‑1 (MQ=255)
                                                                aGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATgg                                                                                    <  8:228542/149‑1 (MQ=255)
                                                                aGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATgg                                                                                    <  6:26575/149‑1 (MQ=255)
                                                                  gACCTGGCCCATCTGCAGCCTGCTGGGGGGCTTTGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTg                                                                                  <  7:305464/149‑1 (MQ=255)
                                                                   aCCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGa                                                                                 <  5:317714/149‑1 (MQ=255)
                                                                          ccATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                <  1:122975/143‑1 (MQ=255)
                                                                          ccATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                >  2:122975/1‑143 (MQ=255)
                                                                           cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAttt                                                                         <  2:347923/149‑1 (MQ=255)
                                                                              cTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTc                                                                       >  7:200119/1‑148 (MQ=255)
                                                                               tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGGGCCATTGCCGGCGGGGGGATTTGTTACGCCCAGGACGGGGACGAAGAGTGTGACGGGATTTTCGa                                                                     >  8:212374/1‑149 (MQ=255)
                                                                               tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGa                                                                                          >  8:93984/1‑128 (MQ=255)
                                                                               tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGa                                                                                          <  7:93984/128‑1 (MQ=255)
                                                                                gCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAc                                                                    >  6:374773/1‑149 (MQ=255)
                                                                                 cAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAcc                                                                   >  7:334458/1‑149 (MQ=255)
                                                                                 cAGCCGGCTGGGCGGCTTGGTCGGGGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAc                                                                    <  2:308335/148‑1 (MQ=255)
                                                                                          gggCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                <  7:258960/127‑1 (MQ=255)
                                                                                          gggCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                >  8:258960/1‑127 (MQ=255)
                                                                                           ggcggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCa                                                         <  8:200119/149‑1 (MQ=255)
                                                                                           ggcggcTTGGGCGGTGGTGGCCATGGGGCCATCGAGAGGTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCa                                                         <  7:135020/149‑1 (MQ=255)
                                                                                                 ttGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGccc                                                   >  2:75265/1‑149 (MQ=255)
                                                                                                     tCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCgag                                               <  5:374773/149‑1 (MQ=255)
                                                                                                              gCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGc                                       <  8:334458/148‑1 (MQ=255)
                                                                                                                    gTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcg                                >  2:214931/1‑149 (MQ=255)
                                                                                                                       ccATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcggcg                             <  6:165181/149‑1 (MQ=255)
                                                                                                                           cGAGAGTTCTTCCGGGGGCGGCGGTGGTGGCGGCTTGGGCGCCTTTGCCGGCGGGCTGATTTGTTTCGCGGAGGACGGTGACGAAGATGGTGACGGCTTTTTCGACCGGGGCGATCAGTGCCCCGAGACCCCTGCCTACACGGCGGt                           <  4:35508/147‑1 (MQ=255)
                                                                                                                            gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCgg                                        <  8:321316/133‑1 (MQ=255)
                                                                                                                            gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCgg                                        >  7:321316/1‑133 (MQ=255)
                                                                                                                             agagTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTgcc                       >  5:91206/1‑147 (MQ=255)
                                                                                                                              gagTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGAcc                      <  1:70136/149‑1 (MQ=255)
                                                                                                                                   cttcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGAGGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACAt                  <  4:173619/148‑1 (MQ=255)
                                                                                                                                   cttcCTGGGCCGCCGGTGGCGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATg                 >  6:361334/1‑149 (MQ=255)
                                                                                                                                      ccTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATggg               <  2:5024/148‑1 (MQ=255)
                                                                                                                                          ggCCGCCGGTGGCGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTg            <  5:361334/147‑1 (MQ=255)
                                                                                                                                            ccgccgGCAATGGCGCCCAAGGCGCCACCACCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCa        >  7:58479/1‑149 (MQ=38)
                                                                                                                                             cgccgGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCAc       >  5:276489/1‑149 (MQ=255)
                                                                                                                                              gccgGTGGTGGCGCCTTGGGGGCCATTGCCGGCGGGGTGATTTGTTTCGCCCAGGACGGGGACGAATATGGTGCCGGCATTTTAGACCCGGGCGCTCACAGCCCCGCGACCCCGGCCAACACGGGGGTTTACCACATTGGGTGgccccg      >  6:42109/1‑146 (MQ=255)
                                                                                                                                              gccgGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACt      >  7:78524/1‑149 (MQ=255)
                                                                                                                                                 ggtggtGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaa   <  6:276489/149‑1 (MQ=255)
                                                                                                                                                  gtggtgGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCCCTGaaa  >  1:44943/1‑149 (MQ=255)
                                                                                                                                                  gtggtgGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa  <  1:75265/149‑1 (MQ=255)
                                                                                                                                                  gtggtgGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa  <  2:55367/149‑1 (MQ=255)
                                                                                                                                                  gtggtgGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa  <  5:42109/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTGGCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAA  >  NC_002947/1284880‑1285174

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