Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,198,339 A→G 5.0% S215G (AGC→GGC)  gspE → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,198,3390AG 5.0% 311.1 / 3.8 103S215G (AGC→GGC) gspEprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base A (48/48);  new base G (2/3);  total (51/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 104 total aligned reads displayed.

TGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198190‑1198487
                                                                                                                                                      |                                                                                                                                                    
tGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGCCGTACGACTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGAGGAGATGCGCGTGCTCGACAAGCAGGCCa                                                                                                                                                      <  1:1789867/150‑1 (MQ=255)
    ccGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  >  1:1770608/1‑150 (MQ=255)
    ccGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  >  2:55585/1‑150 (MQ=255)
      gCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgc                                                                                                                                                >  1:1299403/1‑150 (MQ=255)
      ccgacgACCTTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgc                                                                                                                                                <  1:1418581/145‑1 (MQ=255)
        cTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgctg                                                                                                                                              >  2:1640777/1‑150 (MQ=255)
        cTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgctg                                                                                                                                              >  1:1282563/1‑150 (MQ=255)
         tGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgcttg                                                                                                                                             <  1:1725429/150‑3 (MQ=255)
          ggACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgc                                                                                                                                            >  2:1633809/1‑150 (MQ=255)
              aTTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgctt                                                                                                                                              <  1:616763/144‑2 (MQ=255)
              aTTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgctt                                                                                                                                              >  2:616763/1‑143 (MQ=255)
              aTTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTg                                                                                                                                        >  2:1045723/1‑150 (MQ=255)
              aTTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTg                                                                                                                                        <  2:605478/150‑1 (MQ=255)
                        aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:1784029/150‑1 (MQ=255)
                        aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:138306/150‑1 (MQ=255)
                              aGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGc                                                                                                                        <  1:1455833/150‑1 (MQ=255)
                              aGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGc                                                                                                                        <  2:1830217/150‑1 (MQ=255)
                                ccccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCg                                                                                                                      >  1:383871/1‑150 (MQ=255)
                                 cccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGcc                                                                                                                       >  1:1490939/1‑148 (MQ=255)
                                 cccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGcc                                                                                                                       <  2:1490939/148‑1 (MQ=255)
                                 cccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCgg                                                                                                                     >  1:384669/1‑150 (MQ=255)
                                  ccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCCGCCTACTCTCGCGGGATAACCTGGGCATGCCGGc                                                                                                                    >  2:441765/1‑150 (MQ=255)
                                  ccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGc                                                                                                                    >  1:1319269/1‑150 (MQ=255)
                                      gACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCg                                                                                                                <  2:1523422/150‑1 (MQ=255)
                                          gCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAt                                                                                                            >  1:945204/1‑150 (MQ=255)
                                           ccGCCTCACCTTGCGCGCGGCGCCGCGCCAATTGGCCGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAtt                                                                                                           <  2:714203/150‑1 (MQ=255)
                                            cGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTAGATAACCTGGGCATGCCGGCAGCGGTATTg                                                                                                          >  2:1016512/1‑150 (MQ=255)
                                            cGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCACGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTg                                                                                                          >  2:1205183/1‑150 (MQ=255)
                                             gCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACATTGACCCGCATACACGGCGAGCTCTAATGTATGCGCTTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGc                                                                                                         >  1:1177912/1‑150 (MQ=255)
                                              cATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCa                                                                                                        <  1:1816501/150‑1 (MQ=255)
                                               aTCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAg                                                                                                       <  2:795629/150‑1 (MQ=255)
                                               aTCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTCGAGATTTAGAGGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAg                                                                                                       >  1:1553413/1‑150 (MQ=255)
                                                tCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAgg                                                                                                      <  2:1317330/150‑1 (MQ=255)
                                                tCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAgg                                                                                                      <  2:765328/150‑1 (MQ=255)
                                                tCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAgg                                                                                                      <  2:1770608/150‑1 (MQ=255)
                                                  aCCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGc                                                                                                    <  1:73673/150‑1 (MQ=255)
                                                  aCCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGc                                                                                                    >  2:73673/1‑150 (MQ=255)
                                                    tttGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGGTGCGCGTGCACGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCt                                                                                                  <  2:359657/149‑1 (MQ=255)
                                                     ttGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTg                                                                                                 >  2:14256/1‑150 (MQ=255)
                                                            cggcggGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGtt                                                                                          <  1:1640777/150‑1 (MQ=255)
                                                                  ggCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTgg                                                                                    <  1:55585/150‑1 (MQ=255)
                                                                       gAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCgg                                                                               <  1:1158950/150‑1 (MQ=255)
                                                                       gAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCgg                                                                               <  1:1186755/150‑1 (MQ=255)
                                                                             gACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCAGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAAc                                                                         <  1:441765/150‑1 (MQ=255)
                                                                                gTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGc                                                                      >  1:380921/1‑150 (MQ=255)
                                                                                  aCGAGTCTCGACCGTGCCCGTCATCCACTGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCTTTCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAt                                                                    <  1:594710/150‑1 (MQ=255)
                                                                                  aCGAGGCTCGACTTTGTCCGGCATCCATTGTGAGCGGGTGGTGATGCGGGTGCTTGACAAGCGGGCTAGTCTGCTGGCGCTGGTTAACCTGGGGATGGCGGCAGCGGTATTGCAGGGCCTGCGCCGTTGCCTGGCGCGGCCCAACGGCAt                                                                    <  2:629379/150‑1 (MQ=255)
                                                                                   ggAGTCTCGACCTTGCCCGGCATCCACGGCTAGCGGGTGGTGATGCGCGTGCTCGCCAATCGGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGNAtt                                                                   <  2:1203819/149‑1 (MQ=255)
                                                                                    gAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtg                                                                  >  2:360036/1‑150 (MQ=255)
                                                                                    gAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtg                                                                  >  2:1158378/1‑150 (MQ=255)
                                                                                      gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  1:12403/150‑1 (MQ=255)
                                                                                      gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  1:360036/150‑1 (MQ=255)
                                                                                        ctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACTGCATTGTGCt                                                              <  1:1205183/150‑1 (MQ=255)
                                                                                        ctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAGTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              <  2:503154/150‑1 (MQ=255)
                                                                                         tcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCACCGGCATTGTGCTg                                                             >  2:1295642/1‑150 (MQ=255)
                                                                                            aCCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTcc                                                          >  1:1093954/1‑150 (MQ=255)
                                                                                                 cccccGCATCCACGGCGCGCGGGTGGGGATGCGCGCGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATACCCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCTTGGCGCGGCCCAACGGCATTGTGCTGTCCACCgg                                                     <  1:1272287/149‑1 (MQ=255)
                                                                                                   ccGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTTGGCATGCCGGCAGCGGTATTGCAGGGCCTCCGCAGTTGCCTGGGGCGGCCCCGCGGGATGGTGCGGTCCACCgggg                                                   >  2:286508/1‑149 (MQ=255)
                                                                                                      gCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGTCAGCCTGCTGGCGCTGGATAACCTGGGCGTGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGa                                                <  1:1562926/150‑1 (MQ=255)
                                                                                                             cGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         <  2:1914534/150‑1 (MQ=255)
                                                                                                               gCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCgg                                       <  2:945204/150‑1 (MQ=255)
                                                                                                                 gAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCt                                                                                      >  2:1136323/1‑101 (MQ=255)
                                                                                                                 gAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCt                                                                                      <  1:1136323/101‑1 (MQ=255)
                                                                                                                   gCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGCGCCTGCGCAGTTGTCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGGCCGGTTCGGGCaa                                   >  2:1031667/1‑150 (MQ=255)
                                                                                                                       gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGAcc                               >  2:1708325/1‑150 (MQ=255)
                                                                                                                           tgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGa                           >  1:788395/1‑150 (MQ=255)
                                                                                                                            gATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAc                          >  2:295598/1‑150 (MQ=255)
                                                                                                                              aTGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAcc                         <  2:1123134/150‑1 (MQ=255)
                                                                                                                              aTGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAcc                         >  2:1555999/1‑150 (MQ=255)
                                                                                                                                gcgcGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCtt                       <  1:1728377/150‑1 (MQ=255)
                                                                                                                                  gcgTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGt                     >  1:1831206/1‑150 (MQ=255)
                                                                                                                                   cgTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTa                    <  1:1045723/150‑1 (MQ=255)
                                                                                                                                      gCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGc                 >  1:1312826/1‑150 (MQ=255)
                                                                                                                                        tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGc                                                   <  1:366259/114‑1 (MQ=255)
                                                                                                                                        tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGc                                                   >  2:366259/1‑114 (MQ=255)
                                                                                                                                        tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               >  2:642584/1‑150 (MQ=255)
                                                                                                                                        tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               <  2:1282563/150‑1 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         >  1:1280879/1‑122 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         <  2:1280879/122‑1 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         >  1:1280916/1‑122 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         <  2:1280916/122‑1 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGc             >  1:310507/1‑150 (MQ=255)
                                                                                                                                           aCAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            >  2:718956/1‑150 (MQ=255)
                                                                                                                                            ccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCg                                        >  1:1861763/3‑121 (MQ=255)
                                                                                                                                              aGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGa         >  2:630649/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCgg                                                                               >  1:285474/1‑79 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCgg                                                                               <  2:285474/79‑1 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        >  1:1527780/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        <  1:1016512/150‑1 (MQ=255)
                                                                                                                                                cAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            >  1:1492104/1‑145 (MQ=255)
                                                                                                                                                cAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            <  2:1492104/145‑1 (MQ=255)
                                                                                                                                                cAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGCTGCCTGGCGCGGCCCAACGACATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGCCCACGACCTTCTACGCCAGCCTGAAc       <  2:384669/150‑1 (MQ=255)
                                                                                                                                                   gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    <  2:383871/150‑1 (MQ=255)
                                                                                                                                                   gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    >  2:248022/1‑150 (MQ=255)
                                                                                                                                                   gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    >  1:4089/1‑150 (MQ=255)
                                                                                                                                                    ccAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCCGCCTGAACAGcc   >  1:176866/1‑150 (MQ=255)
                                                                                                                                                    ccAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGcc   >  1:953281/1‑150 (MQ=255)
                                                                                                                                                     cAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGAc                                               <  2:583387/105‑1 (MQ=255)
                                                                                                                                                     cAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGAc                                               >  1:583387/1‑105 (MQ=255)
                                                                                                                                                     cAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCt  >  1:1301633/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: