Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,285,029 G→A 6.6% G87D (GGT→GAT) ‡ PP_1122 → OmpA family protein

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

TGGCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAA  >  NC_002947/1284882‑1285174
                                                                                                                                                   |                                                                                                                                                  
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)
       cAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTAGCCCATCTGCAGCCTGGTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑t                                                                                                                                           >  1:111056/1‑149 (MQ=255)
         gCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTg                                                                                                                                         >  4:34348/1‑149 (MQ=255)
              cTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGGTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGgc                                                                                                                                      <  3:291449/147‑1 (MQ=255)
               tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGgtggtg                                                                                                                                                  >  1:292400/1‑135 (MQ=255)
               tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGgtggtg                                                                                                                                                  <  2:292400/135‑1 (MQ=255)
               tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGgcgc                                                                                                                                    >  7:228542/1‑148 (MQ=255)
                 aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCAt                                                                                                                                 >  7:210102/1‑149 (MQ=255)
                 aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCAt                                                                                                                                 <  3:224222/149‑1 (MQ=255)
                  cGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCAtt                                                                                                                                <  3:186344/149‑1 (MQ=255)
                   gggCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCGTTGGTGCCATCGAGAGTTCTTCATGGGCCGCCGGTGGTGGCGCC‑GTGGGCGCCATTg                                                                                                                               >  8:135020/1‑149 (MQ=255)
                    ggCTGCGCTACGCACAGCTATGGCTGCGCGCCCCTCAATCAAAGGCCCTGGCCCTTCTGCCGCCTTCTGGGGGGATTGGTGGGTGTTGGCCTTGTTGCAATGGATATTTCTTCCTGGGCCGCCTGTGGTGGCGCC‑CAAGGCGCCaccac                                                                                                                              <  4:185316/149‑5 (MQ=255)
                       tGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGCCCGCCGGCAATGGCGCC‑CAAGGCGCCaccaccgg                                                                                                                           >  1:130655/1‑142 (MQ=38)
                          gcTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCcggcgg                                                                                                                        <  8:210102/149‑1 (MQ=255)
                           cTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGcc                                                                                                                                   <  6:368326/137‑1 (MQ=255)
                           cTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGcc                                                                                                                                   >  5:368326/1‑137 (MQ=255)
                            tACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGGGGAATATTACATGGCCCCCGGGTGTTGCGCCGTTGGCGGCCATTGCcgccggc                                                                                                                       >  1:308335/1‑148 (MQ=255)
                              cGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTg                                                                                                                    <  2:111056/149‑1 (MQ=255)
                              cGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGCCCGCCGGCAATGGCGCC‑CAAGGCGCCACCACcggcggccca                                                                                                                    <  2:130655/149‑5 (MQ=37)
                                cacaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCCCC‑TTGGGCTCCATTTCCCGCCGGCAGAt                                                                                                                  >  8:305464/1‑149 (MQ=255)
                                 acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGAtt                                                                                                                 >  5:26575/1‑149 (MQ=255)
                                 acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGAtt                                                                                                                 <  4:27740/149‑1 (MQ=255)
                                 acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGAtt                                                                                                                 >  1:5024/1‑149 (MQ=255)
                                 acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGCAATGGCGCC‑CAAGGCGCCACCACcggcggcccagga                                                                                                                 <  8:297639/149‑8 (MQ=37)
                                 acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGCAATGGCGCC‑CAAGGCGCCACCACcggcggcccagg                                                                                                                  >  7:297639/1‑142 (MQ=37)
                                  caGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                >  3:104679/1‑149 (MQ=255)
                                      gATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTa                                                                                                            >  3:173619/1‑149 (MQ=255)
                                          gCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGcc                                                                                                        <  4:161034/149‑1 (MQ=255)
                                             gcgcgcCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAg                                                                                                     <  4:6783/149‑1 (MQ=255)
                                                 gcCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTGGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACg                                                                                                 <  3:34348/149‑1 (MQ=255)
                                                  ccccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACgg                                                                                                <  5:339607/149‑1 (MQ=255)
                                                  ccccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACgg                                                                                                >  6:317714/1‑149 (MQ=255)
                                                   cccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGt                                                                                               <  3:55148/149‑1 (MQ=255)
                                                     cTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                             >  8:227144/1‑149 (MQ=255)
                                                     cTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                             >  5:165181/1‑149 (MQ=255)
                                                        aatcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACg                                                                                           >  1:347923/1‑148 (MQ=255)
                                                          tcaaAGGACCTGGCCCATCTGCAGCCTGCTCGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑AGGGGCCCCAATTCCGGCCGGCTGAATTGGTTCGCCCACGCGGGTGACg                                                                                           >  4:698/1‑146 (MQ=255)
                                                             aaGGACCTGGCCCATCTGCAGGCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                     <  4:104679/149‑1 (MQ=255)
                                                             aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                     <  3:156606/149‑1 (MQ=255)
                                                             aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                     <  7:316007/149‑1 (MQ=255)
                                                              aGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATgg                                                                                    <  8:228542/149‑1 (MQ=255)
                                                              aGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATgg                                                                                    <  6:26575/149‑1 (MQ=255)
                                                                gACCTGGCCCATCTGCAGCCTGCTGGGGGGCTTTGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTg                                                                                  <  7:305464/149‑1 (MQ=255)
                                                                 aCCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGa                                                                                 <  5:317714/149‑1 (MQ=255)
                                                                        ccATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                <  1:122975/143‑1 (MQ=255)
                                                                        ccATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                >  2:122975/1‑143 (MQ=255)
                                                                         cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAttt                                                                         <  2:347923/149‑1 (MQ=255)
                                                                            cTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTc                                                                       >  7:200119/1‑148 (MQ=255)
                                                                             tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGGGCCATTGCCGGCGGGGGGATTTGTTACGCCCAGGACGGGGACGAAGAGTGTGACGGGATTTTCGa                                                                     >  8:212374/1‑149 (MQ=255)
                                                                             tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGa                                                                                          >  8:93984/1‑128 (MQ=255)
                                                                             tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGa                                                                                          <  7:93984/128‑1 (MQ=255)
                                                                              gCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAc                                                                    >  6:374773/1‑149 (MQ=255)
                                                                               cAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAcc                                                                   >  7:334458/1‑149 (MQ=255)
                                                                               cAGCCGGCTGGGCGGCTTGGTCGGGGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAc                                                                    <  2:308335/148‑1 (MQ=255)
                                                                                        gggCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                <  7:258960/127‑1 (MQ=255)
                                                                                        gggCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGAc                                                                                >  8:258960/1‑127 (MQ=255)
                                                                                         ggcggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCa                                                         <  8:200119/149‑1 (MQ=255)
                                                                                         ggcggcTTGGGCGGTGGTGGCCATGGGGCCATCGAGAGGTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCa                                                         <  7:135020/149‑1 (MQ=255)
                                                                                               ttGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGccc                                                   >  2:75265/1‑149 (MQ=255)
                                                                                                   tCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCgag                                               <  5:374773/149‑1 (MQ=255)
                                                                                                            gCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGc                                       <  8:334458/148‑1 (MQ=255)
                                                                                                                  gTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcg                                >  2:214931/1‑149 (MQ=255)
                                                                                                                     ccATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcggcg                             <  6:165181/149‑1 (MQ=255)
                                                                                                                         cGAGAGTTCTTCCGGGGGCGGCGGTGGTGGCGGC‑TTGGGCGCCTTTGCCGGCGGGCTGATTTGTTTCGCGGAGGACGGTGACGAAGATGGTGACGGCTTTTTCGACCGGGGCGATCAGTGCCCCGAGACCCCTGCCTACACGGCGGt                           <  4:35508/147‑1 (MQ=255)
                                                                                                                          gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCgg                                        <  8:321316/133‑1 (MQ=255)
                                                                                                                          gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCgg                                        >  7:321316/1‑133 (MQ=255)
                                                                                                                           agagTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTgcc                       >  5:91206/1‑147 (MQ=255)
                                                                                                                            gagTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGAcc                      <  1:70136/149‑1 (MQ=255)
                                                                                                                                 cttcCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGAGGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACAt                  <  4:173619/148‑1 (MQ=255)
                                                                                                                                 cttcCTGGGCCGCCGGTGGCGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATg                 >  6:361334/1‑149 (MQ=255)
                                                                                                                                    ccTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATggg               <  2:5024/148‑1 (MQ=255)
                                                                                                                                        ggCCGCCGGTGGCGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTg            <  5:361334/147‑1 (MQ=255)
                                                                                                                                          ccgccgGCAATGGCGCC‑CAAGGCGCCACCACCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCa        >  7:58479/1‑149 (MQ=38)
                                                                                                                                           cgccgGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCAc       >  5:276489/1‑149 (MQ=255)
                                                                                                                                            gccgGTGGTGGCGCC‑TTGGGGGCCATTGCCGGCGGGGTGATTTGTTTCGCCCAGGACGGGGACGAATATGGTGCCGGCATTTTAGACCCGGGCGCTCACAGCCCCGCGACCCCGGCCAACACGGGGGTTTACCACATTGGGTGgccccg      >  6:42109/1‑146 (MQ=255)
                                                                                                                                            gccgGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACt      >  7:78524/1‑149 (MQ=255)
                                                                                                                                               ggtggtGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaa   <  6:276489/149‑1 (MQ=255)
                                                                                                                                                gtggtgGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCCCTGaaa  >  1:44943/1‑149 (MQ=255)
                                                                                                                                                gtggtgGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa  <  5:42109/149‑1 (MQ=255)
                                                                                                                                                gtggtgGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa  <  1:75265/149‑1 (MQ=255)
                                                                                                                                                gtggtgGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa  <  2:55367/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCC‑TTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAA  >  NC_002947/1284882‑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: ‑