Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,199,573 C→G 5.4% 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 5.4% 441.2 / 11.0 140A141G (GCC→GGC) gspFprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (59/73);  new base G (5/3);  total (64/76)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

CAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGCCCAAACT  >  NC_002947/1199425‑1199720
                                                                                                                                                    |                                                                                                                                                   
cAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGcc                                                                                                                                                    <  1:1274121/150‑1 (MQ=255)
   cctgccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAg                                                                                                                                                 <  2:834412/150‑1 (MQ=255)
   cctgccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAg                                                                                                                                                 >  1:317879/1‑150 (MQ=255)
    ctgccctgcACCGCGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCCGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAgg                                                                                                                                                >  2:938084/1‑150 (MQ=255)
    ctgccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAgg                                                                                                                                                <  1:1436480/150‑1 (MQ=255)
     tgccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGt                                                                                                                                               >  2:2025292/1‑150 (MQ=255)
     tgccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGt                                                                                                                                               >  2:12961/1‑150 (MQ=255)
      gccctgcACAGTGTGCTTGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTg                                                                                                                                              >  1:1078071/1‑150 (MQ=255)
       ccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGc                                                                                                                                             <  1:859217/150‑1 (MQ=255)
       ccctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTCTCCTGCGCGCTTGTCGATGCCGGCAAACGTTCCGGGCACCTGGCCCAGGTGc                                                                                                                                             <  2:273162/150‑1 (MQ=255)
        cctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGTGcc                                                                                                                                            <  1:1815939/150‑2 (MQ=255)
         ctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGTgatc                                                                                                                                           <  1:83199/150‑4 (MQ=255)
         ctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGGCCAGGTGccc                                                                                                                                           <  1:2025292/150‑3 (MQ=255)
         ctgcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTg                                                                                                                                           >  2:586910/1‑150 (MQ=255)
           gcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGAc                                                                                                                                         <  1:333932/150‑1 (MQ=255)
           gcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGAc                                                                                                                                         <  2:1639474/150‑1 (MQ=255)
           gcACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCNGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTCGTCGCGGCCGGCGCACGTTTCGGGCACCTGGCCCCGGTGCTGAc                                                                                                                                         >  1:2134117/1‑150 (MQ=255)
            cacaGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGa                                                                                                                                          <  2:1443462/148‑1 (MQ=255)
            cacaGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGa                                                                                                                                          >  1:1443462/1‑148 (MQ=255)
             acaGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGAccc                                                                                                                                       <  2:317879/150‑1 (MQ=255)
             acaGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGAccc                                                                                                                                       <  1:2360979/150‑1 (MQ=255)
               aGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGc                                                                                                                                     <  2:1431686/150‑1 (MQ=255)
                      tggtggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCg                                                                                                                              >  2:2198643/1‑150 (MQ=255)
                       ggtggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGa                                                                                                                             <  1:692997/150‑1 (MQ=255)
                        gtggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGac                                                                                                                            <  1:2198643/150‑1 (MQ=255)
                         tggCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGacc                                                                                                                           >  1:1709353/1‑150 (MQ=255)
                           gCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGACTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGaccac                                                                                                                         >  1:1133735/1‑150 (MQ=255)
                            ccTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGaccacc                                                                                                                        >  1:640877/1‑150 (MQ=255)
                             cTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCt                                                                                                                       >  1:670284/1‑150 (MQ=255)
                              ttGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGaccac                                                                                                                         >  2:1897481/1‑147 (MQ=255)
                              ttGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGaccac                                                                                                                         <  1:1897481/147‑1 (MQ=255)
                              ttGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTg                                                                                                                      >  1:1015372/1‑150 (MQ=255)
                                 cgcgGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAg                                                                                                                   <  2:1706925/150‑1 (MQ=255)
                                  gcgGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGc                                                                                                                  <  2:612021/150‑1 (MQ=255)
                                  gcgGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCGGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGc                                                                                                                  <  2:704369/150‑1 (MQ=255)
                                     gCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAgg                                                                                                               <  1:700484/150‑1 (MQ=255)
                                     gCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAgg                                                                                                               <  2:1159956/150‑1 (MQ=255)
                                         aCTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGcg                                                                                                           >  2:1739106/1‑149 (MQ=255)
                                         aCTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCa                                                                                                           >  2:156755/1‑150 (MQ=255)
                                           tCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAAc                                                                                                         >  2:654041/1‑150 (MQ=255)
                                              ccGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGcc                                                                                                      <  2:1709353/150‑1 (MQ=255)
                                               cGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGACGACCAACTGGAGCAGGTGCAACGCCa                                                                                                     <  1:1627648/150‑1 (MQ=255)
                                                gAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCCAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAg                                                                                                    <  2:1658313/150‑1 (MQ=255)
                                                      ctgggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAacac                                                                                              >  2:696487/1‑150 (MQ=255)
                                                        gggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGac                                                                                                                            >  1:680486/1‑118 (MQ=255)
                                                        gggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGac                                                                                                                            <  2:680486/118‑1 (MQ=255)
                                                         ggcctggCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAg                                                                                           <  2:2178104/150‑1 (MQ=255)
                                                                 ccGCAGCCTGGCCCGGCCGGGCGCACCGTTCACTGGCTTGTAATGCGCGCTGGTCGCGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAATGACCGTAACACAAggacagcgc                                                                                   >  2:178013/1‑143 (MQ=255)
                                                                 ccGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGa                                                                                                                                          <  1:1662755/95‑1 (MQ=255)
                                                                 ccGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGa                                                                                                                                          >  2:1662755/1‑95 (MQ=255)
                                                                   gCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCtt                                                                                                                                   <  1:2104395/100‑1 (MQ=255)
                                                                   gCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCtt                                                                                                                                   >  2:2104395/1‑100 (MQ=255)
                                                                       ccTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCt                                                                             >  1:121889/1‑150 (MQ=255)
                                                                         tGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTAATGAACGCAGGTAGAAAAGGATAAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGa                                                                           >  1:45874/1‑150 (MQ=255)
                                                                             ccGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTa                                                                       >  2:1114492/1‑150 (MQ=255)
                                                                             ccGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTa                                                                       >  1:442851/1‑150 (MQ=255)
                                                                              cGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTAc                                                                      >  2:604903/1‑150 (MQ=255)
                                                                                      cgcACCGTTCACTGTCTTGGTCTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCAGATTTGCCGACCACCCGGAGCAGGTGCATCGCCTGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCtgt                                                              <  1:1993295/150‑1 (MQ=255)
                                                                                          ccGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTg                                                          >  2:1917323/1‑150 (MQ=255)
                                                                                          ccGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTg                                                          >  2:1314409/1‑150 (MQ=255)
                                                                                           cGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGa                                                         >  2:1830403/1‑150 (MQ=255)
                                                                                           cGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGa                                                         <  1:586910/150‑1 (MQ=255)
                                                                                            gTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGAt                                                        <  1:1314409/150‑1 (MQ=255)
                                                                                             ttCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACATAGGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATg                                                       >  1:20916/1‑150 (MQ=255)
                                                                                             ttCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATg                                                       >  2:696463/1‑150 (MQ=255)
                                                                                              tCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTCCCGGGCACCGGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATgg                                                      <  2:964549/150‑1 (MQ=255)
                                                                                              gcACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATgg                                                      <  1:267112/149‑1 (MQ=255)
                                                                                                aCTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGgt                                                    <  1:1516628/150‑1 (MQ=255)
                                                                                                 cTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGgtg                                                   >  1:550703/1‑150 (MQ=255)
                                                                                                  tGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCa                                                                                                     <  2:1424499/99‑1 (MQ=255)
                                                                                                  tGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCa                                                                                                     >  1:1424499/1‑99 (MQ=255)
                                                                                                      ttGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAgg                                                                                                               <  2:595229/85‑1 (MQ=255)
                                                                                                      ttGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAgg                                                                                                               >  1:595229/1‑85 (MQ=255)
                                                                                                        gTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCt                                            <  2:284395/150‑1 (MQ=255)
                                                                                                         tACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTg                                           <  1:920912/150‑1 (MQ=255)
                                                                                                            tGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGa                                                                                                                    >  1:1251071/1‑74 (MQ=255)
                                                                                                            tGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGcc                                        <  2:2164261/150‑1 (MQ=255)
                                                                                                            tGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGATGCTGACCCGCTTGGCCGACCACCTGGa                                                                                                                    <  2:1251071/74‑1 (MQ=255)
                                                                                                            tGCGCGCTGGTCGAGACCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGcc                                        <  2:1037565/150‑1 (MQ=255)
                                                                                                              cgcgcTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCgt                                      <  2:121889/150‑1 (MQ=255)
                                                                                                                 gcTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCgtggt                                   <  1:12961/150‑1 (MQ=255)
                                                                                                                    ggTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGAt                                <  1:696487/150‑1 (MQ=255)
                                                                                                                    ggTCGAGGCCGGCGAACGTTCCGGCCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGTGTGTCGCTGGCCGTGGTGAtt                               <  1:938084/150‑1 (MQ=255)
                                                                                                                       cGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTgg                             <  1:518117/150‑1 (MQ=255)
                                                                                                                        gAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGc                            <  2:2280662/150‑1 (MQ=255)
                                                                                                                        gAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGc                            <  1:920194/150‑1 (MQ=255)
                                                                                                                          ggCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACTCAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGgtgt                                                  >  1:1973412/1‑126 (MQ=255)
                                                                                                                          ggCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGgtgt                                                  <  2:1973412/126‑1 (MQ=255)
                                                                                                                            ccGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCtga                        <  2:1015757/150‑1 (MQ=255)
                                                                                                                            ccGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGAGCCGCGTCAGGGATACAAGGCCCGCACAGCGCTGATCTACCCCACTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCtga                        <  2:1380612/150‑1 (MQ=255)
                                                                                                                             cGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATGTACCCCGCTGTGCTGATGCGTGTGTCGCTGGCCGTGGTGATTGGCCtgat                       >  2:1250315/1‑150 (MQ=255)
                                                                                                                                  aaCGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCACCACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGt                  <  1:156755/150‑1 (MQ=255)
                                                                                                                                     gTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTtcg               >  1:2355934/1‑150 (MQ=255)
                                                                                                                                     gTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTtcg               >  2:867781/1‑150 (MQ=255)
                                                                                                                                           ggCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGc         <  1:1582171/150‑1 (MQ=255)
                                                                                                                                            gCACCTGGGCCAGGTGCTGACCCGCTTGGCCGACCACCTGGCGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGcc        >  2:617118/1‑150 (MQ=255)
                                                                                                                                              aCCTGGCCCAGGTGCTGACCTGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGc                                                            >  1:1669566/1‑96 (MQ=255)
                                                                                                                                                cTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGCCCaaa    >  1:1692049/1‑150 (MQ=255)
                                                                                                                                                  ggCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGCCCAAACt  <  1:604903/150‑1 (MQ=255)
                                                                                                                                                   gCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGgt                                                    <  1:1321373/99‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCTGGCCCGGCAGGGCGCACCGTTCACTGGCTTGTACTGCGCGCTGGTCGAGGCCGGCGAACGTTCCGGGCACCTGGCCCAGGTGCTGACCCGCTTGGCCGACCACCTGGAGCAGGTGCAACGCCAGCAACACAAGGCCCGCACAGCGCTGATCTACCCCGCTGTGCTGATGGGTGTGTCGCTGGCCGTGGTGATTGGCCTGATGACGTTCGTCGTGCCCAAACT  >  NC_002947/1199425‑1199720

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: