Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,198,338 C→T 5.6% 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 5.6% 291.1 / 4.1 99A214A (GCC→GCTgspEprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (45/48);  new base T (2/4);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 105 total aligned reads displayed.

ATGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198189‑1198487
                                                                                                                                                      |                                                                                                                                                     
aTGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGcc                                                                                                                                                       >  1:1073421/1‑150 (MQ=255)
 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:138306/150‑1 (MQ=255)
                         aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:1784029/150‑1 (MQ=255)
                         aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:1633809/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)
                                              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                                                                                                    >  2:73673/1‑150 (MQ=255)
                                                   aCCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGc                                                                                                    <  1:73673/150‑1 (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)
                                                                               aCGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGAGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACg                                                                        <  1:788596/150‑1 (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)
                                                                                                                               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                                                   >  2:366259/1‑114 (MQ=255)
                                                                                                                                         tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGc                                                   <  1:366259/114‑1 (MQ=255)
                                                                                                                                         tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               >  2:642584/1‑150 (MQ=255)
                                                                                                                                         tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               <  2:1282563/150‑1 (MQ=255)
                                                                                                                                           gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         >  1:1280916/1‑122 (MQ=255)
                                                                                                                                           gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         >  1:1280879/1‑122 (MQ=255)
                                                                                                                                           gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTc                                         <  2:1280879/122‑1 (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)
                                                                                                                                             ccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCg                                        <  2:1861763/119‑1 (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:1016512/150‑1 (MQ=255)
                                                                                                                                                gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        >  1:1527780/1‑150 (MQ=255)
                                                                                                                                                 cAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            <  2:1492104/145‑1 (MQ=255)
                                                                                                                                                 cAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            >  1:1492104/1‑145 (MQ=255)
                                                                                                                                                 cAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGCTGCCTGGCGCGGCCCAACGACATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGCCCACGACCTTCTACGCCAGCCTGAAc       <  2:384669/150‑1 (MQ=255)
                                                                                                                                                    gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    >  2:248022/1‑150 (MQ=255)
                                                                                                                                                    gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    >  1:4089/1‑150 (MQ=255)
                                                                                                                                                    gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    <  2:383871/150‑1 (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)
                                                                                                                                                      |                                                                                                                                                     
ATGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198189‑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: