Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,199,573 C→G 9.3% A141G (GCC→GGC)  gspF → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,199,5730CG 9.3% 286.3 / 17.7 97A141G (GCC→GGC) gspFprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (40/48);  new base G (5/4);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

CAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGCCCAAA  >  NC_002947/1199425‑1199718
                                                                                                                                                    |                                                                                                                                                 
cAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGcc                                                                                                                                                  <  2:648809/150‑1 (MQ=255)
   cctgccctgcACATTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAg                                                                                                                                               >  2:286214/1‑150 (MQ=255)
   cctgccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAg                                                                                                                                               >  2:2271823/1‑150 (MQ=255)
        cctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCt                                                                                                                                          >  2:1056476/1‑150 (MQ=255)
         ctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGTGccc                                                                                                                                         <  1:778167/150‑3 (MQ=255)
         ctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGAGCACCTGGGCCAGGTGccc                                                                                                                                         <  2:1776362/150‑3 (MQ=255)
          tgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGa                                                                                                                                        <  1:2188863/150‑1 (MQ=255)
          tgcACAGTGTGCTGGTGGCCTTGCGCGGCTAACTCGCNGAAGGCCTGNGCATGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGTGcccg                                                                                                                                        >  1:1068750/1‑147 (MQ=255)
          tgcACAGTGTGCTGGTGGACTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGa                                                                                                                                        <  2:1124833/150‑1 (MQ=255)
                gtgtGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCt                                                                                                                                  <  1:1056476/150‑1 (MQ=255)
                 tgtgCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCtt                                                                                                                                 <  2:958672/150‑1 (MQ=255)
                 tgtgCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCtt                                                                                                                                 <  1:1462753/150‑1 (MQ=255)
                 tgtgCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCCTTGCTTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCtt                                                                                                                                 <  2:1873143/150‑1 (MQ=255)
                      tggtggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGTCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCg                                                                                                                            <  1:2102487/150‑1 (MQ=255)
                        gtggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTg                                                                                                                                <  1:697283/144‑1 (MQ=255)
                        gtggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTg                                                                                                                                >  2:697283/1‑144 (MQ=255)
                        gtggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGac                                                                                                                          >  2:738271/1‑150 (MQ=255)
                         gggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGacc                                                                                                                         <  2:1439016/149‑1 (MQ=255)
                              ttGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGTAGGGTGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGTGGCCGGTCCTGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTg                                                                                                                    <  1:2244086/150‑1 (MQ=255)
                              ttGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTg                                                                                                                    <  1:1900930/150‑1 (MQ=255)
                                  gcgGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGTGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGc                                                                                                                >  1:1043970/1‑150 (MQ=255)
                                    ggCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAg                                                                                                              <  1:1798673/150‑1 (MQ=255)
                                      cTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGTGcccg                                                                                                                                        <  1:1308617/122‑4 (MQ=255)
                                      cTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGTGcccg                                                                                                                                        >  2:1308617/1‑119 (MQ=255)
                                      cTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGt                                                                                                            <  1:510914/150‑1 (MQ=255)
                                       tCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTg                                                                                                           >  1:842088/1‑150 (MQ=255)
                                         aCTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCa                                                                                                         >  2:41119/1‑150 (MQ=255)
                                         aCTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCCTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGCACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCGGACCACCCGGAGCAGGGGCa                                                                                                         >  1:1088429/1‑150 (MQ=255)
                                           tCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAg                                                                                                                                               <  2:1447341/110‑1 (MQ=255)
                                           tCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAg                                                                                                                                               >  1:1447341/1‑110 (MQ=255)
                                           tCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAAc                                                                                                       >  2:874792/1‑150 (MQ=255)
                                           tCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAAc                                                                                                       >  1:1471846/1‑150 (MQ=255)
                                           tCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAAc                                                                                                       <  2:492145/150‑1 (MQ=255)
                                           tCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAAc                                                                                                       >  1:991847/1‑150 (MQ=255)
                                             gCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGc                                                                                                     >  1:2183776/1‑150 (MQ=255)
                                               cGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCa                                                                                                   <  1:1517564/150‑1 (MQ=255)
                                               cGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCa                                                                                                   <  1:1822759/150‑1 (MQ=255)
                                                gAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGACCGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAg                                                                                                  >  1:494537/1‑150 (MQ=255)
                                                gAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAg                                                                                                  >  1:286418/1‑150 (MQ=255)
                                                 aaGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGc                                                                                                 >  2:2118782/1‑150 (MQ=255)
                                                 aaGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGc                                                                                                 <  1:738271/150‑1 (MQ=255)
                                                      gtgggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAacac                                                                                            >  2:1658635/2‑150 (MQ=255)
                                                      ctgggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGcc                                                                                                                                                  >  2:2225993/1‑96 (MQ=255)
                                                      ctgggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGcc                                                                                                                                                  <  1:2225993/96‑1 (MQ=255)
                                                      ctgggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAacac                                                                                            >  1:2261892/1‑150 (MQ=255)
                                                       tgggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAAcaca                                                                                           >  2:1064566/1‑150 (MQ=255)
                                                            ctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGcc                                                                                      >  2:271386/1‑150 (MQ=255)
                                                                cccGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGcc                                                                                      >  1:189789/1‑146 (MQ=255)
                                                                cccGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGcc                                                                                      <  2:189789/146‑1 (MQ=255)
                                                                cccGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGca                                                                                  <  1:1064566/150‑1 (MQ=255)
                                                                 ccGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGcac                                                                                 >  1:275319/1‑150 (MQ=255)
                                                                   gCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAg                                                                               >  2:818200/1‑150 (MQ=255)
                                                                     ggccaggCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAgcg                                                                             <  1:694073/145‑1 (MQ=255)
                                                                       ccTGGCCCGGCAGGGCGCACCCTTCACTGGCTTGTACTGCGCGCTGGTCGAGGGCGGCGAACGTTCCGGGCACCGGGCGCAGGTGCTGACCCCCTTAGCAGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCAAAgcaat                                                                           <  2:1011192/150‑4 (MQ=255)
                                                                         tGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGt                                                                                                            <  2:608456/115‑1 (MQ=255)
                                                                         tGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGt                                                                                                            >  1:608456/1‑115 (MQ=255)
                                                                         tGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGa                                                                         <  1:874792/150‑1 (MQ=255)
                                                                          ggCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAAAATTCCTGGCACCGGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGAt                                                                        <  1:286214/150‑1 (MQ=255)
                                                                              cGGTAGGGCGCCCCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTAc                                                                    <  2:137338/150‑1 (MQ=255)
                                                                                        cACCGTTCACTGGCTTGTACTGCTCGCTGGTCTAGGCCGGCGAACGTTCCGGGCACTTGGCTCAGGTGCTGACCCGCTTGGCCGATCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGc                                                          <  2:323680/150‑1 (MQ=255)
                                                                                        cACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGc                                                          >  2:948738/1‑150 (MQ=255)
                                                                                              tCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATgg                                                    <  2:95826/150‑1 (MQ=255)
                                                                                              tCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCACGCTGTGCTGATgg                                                    <  2:842088/150‑1 (MQ=255)
                                                                                                      ttGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCg                                            >  1:1807295/1‑150 (MQ=255)
                                                                                                      ttGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCg                                            <  1:338360/150‑1 (MQ=255)
                                                                                                        gTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCt                                          <  1:2271823/150‑1 (MQ=255)
                                                                                                        gTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCt                                          <  1:2118782/150‑1 (MQ=255)
                                                                                                         tACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTg                                         <  2:991847/150‑1 (MQ=255)
                                                                                                            tGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGcc                                      >  2:1879951/1‑150 (MQ=255)
                                                                                                                cgcTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTTACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCGGATCTACCCCGCTGTGTTGATGGGGGTGTCGCTGGCCgtgg                                  <  2:707764/150‑1 (MQ=255)
                                                                                                                 gcTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTg                                                                                                           <  1:2305321/76‑1 (MQ=255)
                                                                                                                 gcTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTg                                                                                                           >  2:2305321/1‑76 (MQ=255)
                                                                                                                 gcTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGa                                                                         >  2:776362/1‑110 (MQ=255)
                                                                                                                 gcTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGa                                                                         <  1:776362/110‑1 (MQ=255)
                                                                                                                    ggTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACACGGCCCGCTTCGCGCTTATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGAt                              <  1:1268600/150‑1 (MQ=255)
                                                                                                                     gTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGAtt                             >  1:240550/1‑150 (MQ=255)
                                                                                                                         aGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGcc                         <  1:1687230/150‑1 (MQ=255)
                                                                                                                         aGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGcc                         >  2:1687230/1‑150 (MQ=255)
                                                                                                                         aGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGcc                         <  1:1658635/150‑1 (MQ=255)
                                                                                                                          gggcGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAg                                                                                                              >  1:2405643/4‑64 (MQ=255)
                                                                                                                          ggCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAg                                                                                                              <  2:2405643/64‑1 (MQ=255)
                                                                                                                             cGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCtgat                     >  2:976795/1‑150 (MQ=255)
                                                                                                                                   aCGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGAtt                             >  1:496306/1‑136 (MQ=255)
                                                                                                                                   aCGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGAtt                             <  2:496306/136‑1 (MQ=255)
                                                                                                                                   aCGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGtt               >  2:2361437/1‑150 (MQ=255)
                                                                                                                                    cGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCgtggt                                 <  1:227902/131‑1 (MQ=255)
                                                                                                                                    cGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCgtggt                                 >  2:227902/1‑131 (MQ=255)
                                                                                                                                    cGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTtc              <  1:976795/150‑1 (MQ=255)
                                                                                                                                     gTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTtcg             <  1:41119/150‑1 (MQ=255)
                                                                                                                                     gTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTtcg             >  1:1042885/1‑150 (MQ=255)
                                                                                                                                        gtcagCACCTGGGCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTtcgtcg          >  2:1638383/5‑150 (MQ=255)
                                                                                                                                         tcagCACCTGGGCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTcgtcgt         >  2:1029895/4‑150 (MQ=255)
                                                                                                                                         cGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGTTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTcgtcgt         <  2:1043970/150‑1 (MQ=255)
                                                                                                                                           ggCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGc       <  1:818200/150‑1 (MQ=255)
                                                                                                                                             cACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGccc     <  2:1362877/150‑1 (MQ=255)
                                                                                                                                              aCCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCtgat                     >  1:1307462/1‑133 (MQ=255)
                                                                                                                                              aCCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCtgat                     <  2:1307462/133‑1 (MQ=255)
                                                                                                                                                cTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAac                                                                                              <  1:2070279/58‑1 (MQ=255)
                                                                                                                                                cTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAac                                                                                              >  2:2070279/1‑58 (MQ=255)
                                                                                                                                                cTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGCCCaaa  >  1:2205576/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGCCCAAA  >  NC_002947/1199425‑1199718

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