Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,993,647 Δ1 bp 7.8% intergenic (+42/+56) PP_1780 → / ← PP_1781 mannosyltransferase/acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,993,6470A. 7.8% 248.3 / 17.3 90intergenic (+42/+56)PP_1780/PP_1781mannosyltransferase/acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (37/41);  new base . (4/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CTACATTGGAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAGGCTTAGCGGGA  >  NC_002947/1993507‑1993794
                                                                                                                                             |                                                                                                                                                   
cTACATTGGAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTg                                                                                                                                            <  2:1830845/150‑1 (MQ=255)
 tACATTGGAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGc                                                                                                                                           >  2:950253/1‑150 (MQ=255)
  aCATTGGAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCTGCGCCGCACTGCa                                                                                                                                          <  2:1221260/150‑1 (MQ=255)
  aCATTGGAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCa                                                                                                                                          <  2:1211043/150‑1 (MQ=255)
  aCATTGGAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCa                                                                                                                                          <  2:2129353/150‑1 (MQ=255)
        gAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCtt                                                                                                                                    >  1:2200536/1‑150 (MQ=255)
             aaaGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCTGCGCCGCACTGCAgg                                                                                                                                        <  1:1167163/141‑1 (MQ=255)
             aaaGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCTGCGCCGCACTGCAgg                                                                                                                                        >  2:1167163/1‑141 (MQ=255)
               aGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCat                                                                                                                             <  1:2129729/150‑1 (MQ=255)
               aGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCat                                                                                                                             >  2:63295/1‑150 (MQ=255)
                 tGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTGAAGGGCACGGATTGCGTGCGCTTGGCCGTCCATGAACCCTGCAGTGCGCCGCGGACGCACTGCAAGGCTTTTGCCatat                                                                                                                           >  2:620313/1‑150 (MQ=255)
                  gCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATg                                                                                                                          >  1:654781/1‑150 (MQ=255)
                      cTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACtt                                                                                                                      >  2:1655504/1‑150 (MQ=255)
                        gAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGt                                                                                                                    <  1:1843184/150‑1 (MQ=255)
                          aCTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTAt                                                                                                                  <  2:2177126/150‑1 (MQ=255)
                          aCTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTCAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACCCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTAt                                                                                                                  <  1:453452/150‑1 (MQ=255)
                           cTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGACATATGACTTGtggg                                                                                                                 >  1:1788907/1‑147 (MQ=255)
                                    cTGGTGTTGTATGTTACAGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGAACTGATTGCGTACGCTTGGACATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTAGGAGGCTCGCa                                                                                                        <  2:674011/150‑1 (MQ=255)
                                               tGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCt                                                                                             >  2:1291071/1‑150 (MQ=255)
                                                     ggAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCat                                                                                                                             >  2:1833323/1‑112 (MQ=255)
                                                     ggAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCat                                                                                                                             <  1:1833323/112‑1 (MQ=255)
                                                     ggAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCg                                                                                       <  1:1291071/150‑1 (MQ=255)
                                                     ggAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCg                                                                                       >  2:1643673/1‑150 (MQ=255)
                                                         aaTACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTa                                                                                   >  1:1378863/1‑150 (MQ=255)
                                                             cTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGAACGCAGTAAACGCGCTTATCCGGCTACggg                                                                               <  1:63295/150‑1 (MQ=255)
                                                              tCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGc                                                                              <  2:586999/150‑1 (MQ=255)
                                                               cAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGca                                                                             <  2:654781/150‑1 (MQ=255)
                                                                aGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGcac                                                                            >  2:2147705/1‑150 (MQ=255)
                                                                      aGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTg                                                                      <  2:2104483/150‑1 (MQ=255)
                                                                             ggCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGc                                                               <  2:1212431/150‑1 (MQ=255)
                                                                             ggCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGc                                                               <  1:1655504/150‑1 (MQ=255)
                                                                              gcgcCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCa                                                                                                        <  2:1849112/108‑1 (MQ=255)
                                                                              gcgcCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCa                                                                                                        >  1:1849112/1‑108 (MQ=255)
                                                                                gcCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACggg                                                                               <  2:2186994/131‑1 (MQ=255)
                                                                                gcCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACggg                                                                               >  1:2186994/1‑131 (MQ=255)
                                                                                 ccTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGgc                                                           <  2:861134/150‑1 (MQ=255)
                                                                                   tGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTTTTTTTTTTTTTGTATGAGGCACGCAGTAAACCCTCTAAGCCGGCCACGTTCACAGCCTGGTGGCGATGgcgt                                                         >  2:967147/1‑149 (MQ=255)
                                                                                    gCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCa                                                        >  2:947388/1‑150 (MQ=255)
                                                                                     cTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCa                                                                                                                                          >  2:1613032/1‑67 (MQ=255)
                                                                                     cTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCa                                                                                                                                          <  1:1613032/67‑1 (MQ=255)
                                                                                      tCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGa                                                      >  1:2064207/1‑150 (MQ=255)
                                                                                      tCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGa                                                      >  1:1994676/1‑150 (MQ=255)
                                                                                      tCAAGCGCACAGATCGCGTAAGCTTGGCCATCAATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACATGTATGAGGAACGCAGTAAACGCGATTATCCGGCTAAGGGCACAGCATGGCGACGCTGGCGCAGa                                                      >  2:1646578/1‑150 (MQ=255)
                                                                                        aaGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCa                                                                                                        >  2:294636/1‑98 (MQ=255)
                                                                                        aaGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCa                                                                                                        <  1:294636/98‑1 (MQ=255)
                                                                                        aaGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGt                                                                                                      <  1:446786/100‑1 (MQ=255)
                                                                                        aaGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGt                                                                                                      >  2:446786/1‑100 (MQ=255)
                                                                                        aaGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATa                                                    >  1:2123419/1‑150 (MQ=255)
                                                                                            gcACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGt                                                >  1:1664316/1‑150 (MQ=255)
                                                                                             taCTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCACAGCAGCTTTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGtt                                               <  1:1310782/149‑1 (MQ=255)
                                                                                             cACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGtt                                               >  2:878163/1‑150 (MQ=255)
                                                                                             cACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGtt                                               <  1:2092098/150‑1 (MQ=255)
                                                                                              aCTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTg                                              >  1:1880129/1‑150 (MQ=255)
                                                                                               cTGATTGCGTACGCTTGGCCATCCACGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGa                                             <  2:1994676/150‑1 (MQ=255)
                                                                                                    tGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCgg                                        <  2:2200536/150‑1 (MQ=255)
                                                                                                       gTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTg                                                                      <  2:1627225/117‑1 (MQ=255)
                                                                                                       gTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTg                                                                      >  1:1627225/1‑117 (MQ=255)
                                                                                                       gTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCg                                     >  1:357834/1‑150 (MQ=255)
                                                                                                       gTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCg                                     >  2:1043036/1‑150 (MQ=255)
                                                                                                                 ccatccatGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTccagccag                           >  1:934700/1‑150 (MQ=255)
                                                                                                                  catccatGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCcagccagc                          >  1:1704302/1‑150 (MQ=255)
                                                                                                                   atccatGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCg                         <  1:723263/150‑1 (MQ=255)
                                                                                                                   atcaatGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGa                                             <  1:902613/130‑1 (MQ=255)
                                                                                                                   atcaatGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGa                                             >  2:902613/1‑130 (MQ=255)
                                                                                                                       atGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGc                                                                         <  2:899219/98‑1 (MQ=255)
                                                                                                                       atGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGc                                                                         >  1:899219/1‑98 (MQ=255)
                                                                                                                       atGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTgg                     <  2:2123419/150‑1 (MQ=255)
                                                                                                                        tGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTggg                    >  1:822815/1‑150 (MQ=255)
                                                                                                                        tGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTggg                    <  2:934700/150‑1 (MQ=255)
                                                                                                                         gAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAgg                                                                                                              >  2:1363118/1‑59 (MQ=255)
                                                                                                                         gAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAgg                                                                                                              <  1:1363118/59‑1 (MQ=255)
                                                                                                                           aCCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCa                 >  1:848884/1‑150 (MQ=255)
                                                                                                                             ccTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAAc               <  2:848884/150‑1 (MQ=255)
                                                                                                                              cTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACa              <  2:822815/150‑1 (MQ=255)
                                                                                                                               tGCAGTGCGGCTGCGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGc                                                                                     >  1:1558168/1‑78 (MQ=37)
                                                                                                                               tGCAGTGCGGCTGCGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGa                                             <  2:1482865/118‑1 (MQ=255)
                                                                                                                               tGCAGTGCGGCTGCGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGa                                             >  1:1482865/1‑118 (MQ=255)
                                                                                                                               tGCAGTGCGGCTGCGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAg             >  1:2163565/1‑150 (MQ=255)
                                                                                                                               tGCAGTGCGGCTGCGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAg             >  2:2196763/1‑150 (MQ=255)
                                                                                                                               tGCAGTGCGGCTGCGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGc                                                                                     <  2:1558168/78‑1 (MQ=37)
                                                                                                                               tGCAGTGCGGCTGCGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAg             >  2:223682/1‑150 (MQ=255)
                                                                                                                               tGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTgg                     >  2:1693451/1‑142 (MQ=255)
                                                                                                                               tGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTgg                     <  1:1693451/142‑1 (MQ=255)
                                                                                                                               tGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAg             <  1:620313/150‑1 (MQ=255)
                                                                                                                                  aGTGCGGCTGCGCCGCACTGCAGGGTTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGc                                                  <  1:1218274/110‑1 (MQ=255)
                                                                                                                                    tGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAGGCTTa        <  2:2064207/150‑1 (MQ=255)
                                                                                                                                        gcgcAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAGGCTTAGCgg    <  2:1378863/150‑1 (MQ=255)
                                                                                                                                        gcgcAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAGGCTTAGCgg    <  1:950253/150‑1 (MQ=255)
                                                                                                                                         cgcAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAGGCTTAGCggg   >  2:2083333/1‑150 (MQ=255)
                                                                                                                                           gcAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAGGCTTAGCGGGa  <  2:690606/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CTACATTGGAGTGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCGCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACAGGCTTAGCGGGA  >  NC_002947/1993507‑1993794

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: