Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,198,338 C→T 7.7% A214A (GCC→GCT gspE → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,198,3380CT 7.7% 366.7 / 14.9 130A214A (GCC→GCTgspEprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (54/65);  new base T (3/7);  total (57/73)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 132 total aligned reads displayed.

TGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198190‑1198487
                                                                                                                                                    |                                                                                                                                                     
 ggCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAg                                                                                                                                                     >  2:2038699/1‑150 (MQ=255)
   cccGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGcc                                                                                                                                                   <  2:583760/150‑1 (MQ=255)
        cTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgctt                                                                                                                                              <  1:1879762/150‑2 (MQ=255)
         tGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgcttg                                                                                                                                             <  1:968799/150‑3 (MQ=255)
            aCATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcgc                                                                                                                                          <  1:746196/150‑1 (MQ=255)
            aCATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcgc                                                                                                                                          <  2:2042822/150‑1 (MQ=255)
             cATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgcttgtcga                                                                                                                                         <  2:488830/150‑7 (MQ=255)
              aTTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTg                                                                                                                                        >  1:1549959/1‑150 (MQ=255)
               ttGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTgg                                                                                                                                       >  2:2233330/1‑150 (MQ=255)
                  cGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATa                                                                                                                                    <  1:2382078/150‑1 (MQ=255)
                   ggAAAAGCGTCAGCCCCAGGACGGACGCATCACCTTGCGCGCGGCGGGGCGCGACGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCACGCGCTCGACAAGCAGGCCAGCCTGCTGGCGCTCGATac                                                                                                                                   >  2:2301007/1‑149 (MQ=255)
                    gAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGCATAGCAGTCCAGCCTGCTGGCGCTGGATAAc                                                                                                                                  <  1:1871512/150‑1 (MQ=255)
                     aaaaGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  <  1:319653/133‑1 (MQ=255)
                     aaaaGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  >  2:319653/1‑133 (MQ=255)
                        agcctgggCCCCAGGACGGCCGTCTCACCGTGCGTGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGTCCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAACGAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:924089/143‑1 (MQ=255)
                            tCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCt                                                                                                                                         >  2:2386687/1‑135 (MQ=255)
                            tCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCt                                                                                                                                         <  1:2386687/135‑1 (MQ=255)
                            tCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCAt                                                                                                                          >  1:1675225/1‑150 (MQ=255)
                            tCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCAGGCTGgcctgc                                                                                                                                        <  1:1169767/136‑5 (MQ=255)
                            tCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCAGGCTGgcctgc                                                                                                                                        >  2:1169767/1‑132 (MQ=255)
                                 cccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCgg                                                                                                                     <  2:1039656/150‑1 (MQ=255)
                                   cAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCa                                                                                                                   <  2:1387638/150‑1 (MQ=255)
                                   cAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCa                                                                                                                   >  2:844572/1‑150 (MQ=255)
                                   cAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCa                                                                                                                   <  1:2038699/150‑1 (MQ=255)
                                        cGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGt                                                                                                              <  1:67016/150‑1 (MQ=255)
                                            cGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTg                                                                                                          >  2:1598908/1‑150 (MQ=255)
                                                tCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAgg                                                                                                      <  1:588371/150‑1 (MQ=255)
                                                tCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGATCAACAAGAAGGCCAGGCTGTTGGCGCTCGATATCCTAGGCATGTTGGTAGCGGTATTGCAgg                                                                                                      <  2:636911/150‑1 (MQ=255)
                                                  aCCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGc                                                                                                    >  2:2002831/1‑150 (MQ=255)
                                                             ggcggGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTg                                                                                         >  1:1633128/1‑150 (MQ=255)
                                                              gcggGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGc                                                                                        >  1:2150729/1‑150 (MQ=255)
                                                               cggGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATCACCTAGAACTGCGGACAAGGGTAGTGCAGGGCCTGCGCAGTTGcc                                                                                       >  2:633988/1‑150 (MQ=255)
                                                                ggggCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCt                                                                                      <  2:2124189/150‑1 (MQ=255)
                                                                ggggCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGAATTGCAGGGCCTGCGCAGTTGCCt                                                                                      >  2:916551/1‑150 (MQ=255)
                                                                   gcgcGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgc                                                                                   >  2:1897812/1‑150 (MQ=255)
                                                                   gcgcGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgc                                                                                   <  1:1534073/150‑1 (MQ=255)
                                                                     gcgAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAg                                                                                                       <  2:1300103/128‑1 (MQ=255)
                                                                        aaGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGc                                                                              <  1:624700/150‑1 (MQ=255)
                                                                           tGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCa                                                                           <  1:810898/150‑1 (MQ=255)
                                                                             gACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCt                                                                                                                                         >  1:2028977/1‑86 (MQ=255)
                                                                             gACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAAc                                                                         >  1:1028305/1‑150 (MQ=255)
                                                                             gACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAAc                                                                         <  1:1598908/150‑1 (MQ=255)
                                                                             gACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGAGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAAc                                                                         >  2:2393061/1‑150 (MQ=255)
                                                                               cGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACgg                                                                       <  2:1633128/150‑1 (MQ=255)
                                                                                  aCGAGTCTCGACCTTGCCCGGCATCAACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAt                                                                    >  2:2124484/1‑150 (MQ=255)
                                                                                   cGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAtt                                                                   <  2:2284221/150‑1 (MQ=255)
                                                                                    gAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtg                                                                  >  1:2262530/1‑150 (MQ=255)
                                                                                      gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  1:689657/150‑1 (MQ=255)
                                                                                       tctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCgg                                                                                                                     >  2:1567699/1‑96 (MQ=255)
                                                                                       tctcGACCTTGCCCGGCATCCACGGCGAGCGGGCGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCGGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGc                                                               >  1:803369/1‑150 (MQ=255)
                                                                                        ctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAtt                                                                   <  2:2302062/145‑1 (MQ=255)
                                                                                        ctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              <  1:10641/150‑1 (MQ=255)
                                                                                        ctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCCTGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              <  1:2069242/150‑1 (MQ=255)
                                                                                          cGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGt                                                            <  1:1079249/150‑1 (MQ=255)
                                                                                           gACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgcgc                                                                                 <  1:1692689/128‑1 (MQ=255)
                                                                                           gACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgcgc                                                                                 >  2:1692689/1‑128 (MQ=255)
                                                                                            aCCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTcc                                                          >  2:1279120/1‑150 (MQ=255)
                                                                                             ccTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAAc                                                                                                                                  >  1:87572/1‑77 (MQ=255)
                                                                                              cTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCt                                                                                                                                <  2:621008/78‑1 (MQ=255)
                                                                                              cTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCt                                                                                                                                >  1:621008/1‑78 (MQ=255)
                                                                                               ttGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCAcc                                                       >  1:1449845/1‑150 (MQ=255)
                                                                                                   ccGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcgc                                                                                                                                          >  2:1666781/1‑63 (MQ=255)
                                                                                                   ccGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGc                                                   <  2:1325593/150‑1 (MQ=255)
                                                                                                      gCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGa                                                >  2:185232/1‑150 (MQ=255)
                                                                                                       cATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGAc                                               >  1:1388199/1‑150 (MQ=255)
                                                                                                        aTCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGAcc                                              >  2:51519/1‑150 (MQ=255)
                                                                                                           cACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGt                                           <  2:104016/150‑1 (MQ=255)
                                                                                                             cGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         <  1:2233330/150‑1 (MQ=255)
                                                                                                              ggCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCg                                        <  2:2030771/150‑1 (MQ=255)
                                                                                                               gCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCgg                                       <  2:1582232/150‑1 (MQ=255)
                                                                                                               ccGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCCCTGCATAACCTGGGCACGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCTAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCgg                                       <  2:772256/149‑1 (MQ=255)
                                                                                                                cGAGCCGGTTGTGAGGCGCTTGCTCGACAAGCCGGCAAGCCTGCTGGCGCTGGATAACCTGGGCATGTCGGCAGCGGTATTGCAGGGCCTGCGCAGTGGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCCCCGGGCCGACCGGTTCggg                                      <  1:410545/150‑1 (MQ=255)
                                                                                                                    cGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAg                                  <  2:1028305/150‑1 (MQ=255)
                                                                                                                       gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGtt                                                                                          >  1:2193217/1‑91 (MQ=255)
                                                                                                                       gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGtt                                                                                          <  2:2193217/91‑1 (MQ=255)
                                                                                                                       gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              <  1:2369196/119‑1 (MQ=255)
                                                                                                                       gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              >  2:2369196/1‑119 (MQ=255)
                                                                                                                       gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGAcc                               >  1:425467/1‑150 (MQ=255)
                                                                                                                        tggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCNCTGGATNACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTTCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCa                              >  1:1222077/1‑150 (MQ=255)
                                                                                                                         ggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCAc                             <  1:1279120/150‑1 (MQ=255)
                                                                                                                         ggtgATGCGCGCGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCAc                             <  1:130106/150‑1 (MQ=255)
                                                                                                                            gATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCg                                                                                                                >  1:1743400/1‑64 (MQ=255)
                                                                                                                            gATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCg                                                                                                                <  2:1743400/64‑1 (MQ=255)
                                                                                                                            gATGCGCGTGCCCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGGCCGCCGTTTCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAc                          <  2:2019591/150‑1 (MQ=255)
                                                                                                                             aTGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAcc                         <  2:1034665/150‑1 (MQ=255)
                                                                                                                               gcgcGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCtt                       >  2:896576/1‑150 (MQ=255)
                                                                                                                               gcgcGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCtt                       >  2:896735/1‑150 (MQ=255)
                                                                                                                               gcgcGTGCTCGACAAGCAGGCCAGCCTCATCGAGATCTATAACCTTGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCtt                       <  2:803369/150‑1 (MQ=255)
                                                                                                                                   gTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTAc                   <  1:1547061/150‑1 (MQ=255)
                                                                                                                                   gTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTAc                   <  1:1105639/150‑1 (MQ=255)
                                                                                                                                        gcgccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGa                                 >  1:1647598/6‑131 (MQ=255)
                                                                                                                                        gcgccaGCAGGCTGGCCTGCTGCCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGa                                 <  2:1647598/126‑1 (MQ=255)
                                                                                                                                         gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGc             <  2:612152/150‑1 (MQ=255)
                                                                                                                                           ccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCgg                                                                               >  2:793232/3‑82 (MQ=255)
                                                                                                                                           ccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCgg                                                                               <  1:793232/80‑1 (MQ=255)
                                                                                                                                           ccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGc                                                   >  1:792972/3‑110 (MQ=255)
                                                                                                                                           cAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCt           <  1:844572/150‑1 (MQ=255)
                                                                                                                                              gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        <  2:1449845/150‑1 (MQ=255)
                                                                                                                                               cAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAAc       >  1:893107/1‑150 (MQ=255)
                                                                                                                                                    cAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCt  >  1:1851348/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198190‑1198487

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: