Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,627,499 C→T 6.5% intergenic (+30/+36) PP_1421 → / ← PP_1422 sensor histidine kinase TctE/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,627,4990CT 6.5% 819.8 / 31.4 261intergenic (+30/+36)PP_1421/PP_1422sensor histidine kinase TctE/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (135/109);  new base T (6/11);  total (141/120)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 267 total aligned reads displayed.

GGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCG  >  NC_002947/1627350‑1627648
                                                                                                                                                     |                                                                                                                                                     
        gTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCcgc                                                                                                                                               <  2:1865880/150‑1 (MQ=255)
        gTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCcgc                                                                                                                                               <  2:5448299/150‑1 (MQ=255)
         ttAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCcgcg                                                                                                                                              >  2:4583645/1‑150 (MQ=255)
         ttAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCcgcg                                                                                                                                              >  1:3283990/1‑150 (MQ=255)
             ggCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAg                                                                                                                                          <  2:143657/150‑1 (MQ=255)
              gCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAgg                                                                                                                                         <  2:234775/150‑1 (MQ=255)
               cTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCACGTGACCCTGCACGATGGCGAGAAGGTCGGGTTGCGGGTGCGGGTGAGCTTGAGCGAGGATTGATGGTTTGCCTGTGTCTCCCCTTTCATGGGCTCCCCCGCGAAAggg                                                                                                                                        <  2:4673622/150‑1 (MQ=255)
               cTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAggg                                                                                                                                        <  1:1173379/150‑1 (MQ=255)
               cTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAggg                                                                                                                                        <  2:2831806/150‑1 (MQ=255)
               cTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAggg                                                                                                                                        <  2:7181194/150‑1 (MQ=255)
                tcgacaTCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAgggg                                                                                                                                       >  2:2177540/6‑150 (MQ=255)
                 ggCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGc                                                                                                                                      <  2:6604890/150‑1 (MQ=255)
                 ggCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGc                                                                                                                                      >  2:7335803/1‑150 (MQ=255)
                 ggCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGc                                                                                                                                      <  1:292396/150‑1 (MQ=255)
                 ggCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGGGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGc                                                                                                                                      <  2:2638594/150‑1 (MQ=255)
                  gCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCg                                                                                                                                     <  2:3625254/150‑1 (MQ=38)
                   cGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGa                                                                                                                                    <  2:1838776/150‑1 (MQ=38)
                    gATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCggc                                                                                                                                   >  2:4449886/1‑148 (MQ=255)
                    gATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGac                                                                                                                                   >  1:6783137/1‑150 (MQ=255)
                       cGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCAc                                                                                                                                                   >  2:6863110/1‑131 (MQ=255)
                       cGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGAcaca                                                                                                                                <  1:967975/150‑1 (MQ=255)
                         ttGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAgggg                                                                                                                                       <  1:15808/141‑1 (MQ=255)
                           ggCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGAGGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCt                                                                                                                            <  1:279319/150‑1 (MQ=255)
                            gCGAGATCTGTCTGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTg                                                                                                                           <  1:4263075/150‑1 (MQ=255)
                            gCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTg                                                                                                                           <  1:7422378/150‑1 (MQ=255)
                             cGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAGGcaaa                                                                                                                          <  1:5608661/150‑4 (MQ=37)
                             cGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAGGcaaa                                                                                                                          <  2:3741632/150‑4 (MQ=37)
                              gagaTCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGt                                                                                                                         >  1:6845421/1‑150 (MQ=255)
                               agaTCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTa                                                                                                                        <  1:6367907/150‑1 (MQ=255)
                               agaTCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTa                                                                                                                        <  2:6433849/150‑1 (MQ=255)
                                    tGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGacac                                                                                                                                 >  1:5027465/1‑136 (MQ=37)
                                    tGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCt                                                                                                                   >  1:7347367/1‑150 (MQ=255)
                                      tCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCa                                                                                                                 <  1:7071228/150‑1 (MQ=255)
                                        gggCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAg                                                                                                                               <  2:4863250/134‑1 (MQ=37)
                                            gcACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCaa                                                                                                           <  2:3165427/150‑1 (MQ=255)
                                            gcACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCaa                                                                                                           <  1:6863844/150‑1 (MQ=255)
                                              aCCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAACGGGGCGACACAGGCTGGTACAGCTCAGCGCAACa                                                                                                         >  1:5049269/1‑150 (MQ=255)
                                                 tGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGc                                                                                                      >  2:1560453/1‑150 (MQ=255)
                                                       aaGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACCGGCTGGTACAGCTCAGCGCAACATGCTCcgcg                                                                                                >  1:205714/1‑150 (MQ=255)
                                                              ccTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcag                                                                                         <  1:4918673/150‑1 (MQ=255)
                                                              ccTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcag                                                                                         <  1:6664852/150‑1 (MQ=255)
                                                               cTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAgggg                                                                                                                                       <  1:2608337/103‑1 (MQ=255)
                                                               cTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagc                                                                                        >  1:5762028/1‑150 (MQ=255)
                                                                tGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                       >  1:3763321/1‑150 (MQ=255)
                                                                tGCACGATGGCGAGAAGGGCGGGGTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCGGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                       <  2:276860/150‑1 (MQ=255)
                                                                  cACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGt                                                                                     <  1:7301493/150‑1 (MQ=255)
                                                                   aCGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGtt                                                                                    >  2:5514093/1‑150 (MQ=255)
                                                                     gagggCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCa                                                                                  <  2:6905704/147‑1 (MQ=255)
                                                                     gATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCa                                                                                  <  1:284968/150‑1 (MQ=255)
                                                                        ggCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCg                                                                               <  2:2497937/150‑1 (MQ=255)
                                                                           gagaAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGAcc                                                                            >  1:1091761/1‑150 (MQ=255)
                                                                             gaAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCCTGGGCTCACCCGCGCAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCCGTTCATCGACCTc                                                                          >  2:5669047/1‑150 (MQ=255)
                                                                              aaGGGCGGGTTGCGGGTGCGGGTGAGCTTTCTCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTcc                                                                         <  2:29305/150‑1 (MQ=255)
                                                                              aaGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTcc                                                                         <  1:2518825/150‑1 (MQ=255)
                                                                                 ggCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGc                                                                      >  1:537619/1‑150 (MQ=255)
                                                                                 ggCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGc                                                                      >  2:6318794/1‑150 (MQ=255)
                                                                                      gTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGt                                                                 <  1:3360027/150‑1 (MQ=255)
                                                                                         gcgggtgcgggtGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGcc                                                              >  2:3430196/1‑150 (MQ=255)
                                                                                            ggtgcgggtgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCaaa                                                           >  2:5641839/1‑150 (MQ=255)
                                                                                              tgcgggtgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACg                                                         <  2:178187/150‑1 (MQ=255)
                                                                                                    tgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAAcc                                                   >  2:787354/1‑150 (MQ=255)
                                                                                                    tgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAAcc                                                   >  1:6047757/1‑150 (MQ=255)
                                                                                                      aGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCAAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAAcccc                                                 >  1:2825204/1‑150 (MQ=255)
                                                                                                         tttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGa                                              >  2:6359110/1‑150 (MQ=255)
                                                                                                         tttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGa                                              <  2:5996127/150‑1 (MQ=255)
                                                                                                          ttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                       <  1:1323447/108‑1 (MQ=255)
                                                                                                          ttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                       >  2:1323447/1‑108 (MQ=255)
                                                                                                           tATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATg                                            <  1:2097528/150‑1 (MQ=255)
                                                                                                                 ggATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTa                                      >  2:1148920/1‑150 (MQ=255)
                                                                                                                      gATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTGAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCg                                 <  1:6899868/150‑1 (MQ=255)
                                                                                                                       aTGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACa                                                                                                         >  1:3529077/1‑77 (MQ=255)
                                                                                                                       aTGGTTTGCCTGGGTCGCCCCTTTCATGGGCTCACCCGCGACAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCgg                                <  2:3870539/150‑1 (MQ=255)
                                                                                                                           tttGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGgcagc                            >  2:423987/1‑150 (MQ=255)
                                                                                                                                cTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGgcagcgcagc                       >  2:2127856/1‑150 (MQ=255)
                                                                                                                                cTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGgcagcgcagc                       <  2:1598751/150‑1 (MQ=255)
                                                                                                                                  gtgtCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCa                     >  1:2749412/1‑150 (MQ=255)
                                                                                                                                    gtCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAAt                   <  1:2309429/150‑1 (MQ=255)
                                                                                                                                    gtCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCAGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAAt                   >  1:2502965/1‑150 (MQ=255)
                                                                                                                                     tCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGcc                                  <  2:3039601/134‑1 (MQ=37)
                                                                                                                                     tCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATg                  >  2:448470/1‑150 (MQ=38)
                                                                                                                                     tCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATg                  <  1:6275702/150‑1 (MQ=255)
                                                                                                                                     tCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATg                  >  2:4197154/1‑150 (MQ=255)
                                                                                                                                       gCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCGCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCCGCCAATGGc                >  1:149747/1‑150 (MQ=255)
                                                                                                                                       gCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCGCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGGGCCAAACGGCCACCCCAGATGGCGGTCCGCCGgcggcgcggcggctggg                >  2:5372057/1‑143 (MQ=255)
                                                                                                                                        ccccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGcc               >  1:3748462/1‑150 (MQ=255)
                                                                                                                                        ccccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGcc               >  1:1832283/1‑150 (MQ=255)
                                                                                                                                         cccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCCGCGCAGCCAATGGCCg              >  2:6675561/1‑150 (MQ=255)
                                                                                                                                          ccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCgg             >  2:5322637/1‑150 (MQ=255)
                                                                                                                                          ccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCgg             >  1:6404313/1‑150 (MQ=255)
                                                                                                                                          ccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCgg             >  2:50892/1‑150 (MQ=255)
                                                                                                                                          ccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCgg             >  1:2677153/1‑150 (MQ=255)
                                                                                                                                           cTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGt            >  2:5222514/1‑150 (MQ=255)
                                                                                                                                            tttCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTg           >  1:4397822/1‑150 (MQ=255)
                                                                                                                                            tttCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTg           >  1:5082295/1‑150 (MQ=255)
                                                                                                                                                aTGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGc       <  2:5193380/150‑1 (MQ=255)
                                                                                                                                                aTGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGc       <  2:4067500/150‑1 (MQ=255)
                                                                                                                                                 tGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGcc      >  1:4629094/1‑150 (MQ=255)
                                                                                                                                                   ggCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCcg    <  2:249046/150‑1 (MQ=255)
                                                                                                                                                   ggCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCcg    >  1:4076328/1‑150 (MQ=255)
                                                                                                                                                     cTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCcgcg  >  1:614169/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCG  >  NC_002947/1627350‑1627648

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