Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,381,267 T→C 9.8% intergenic (+32/+33) bamE → / ← PP_4732 outer membrane protein assembly factor BamE/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,381,2670TC 9.8% 191.6 / 9.7 71intergenic (+32/+33)bamE/PP_4732outer membrane protein assembly factor BamE/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/29);  new base C (4/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGA  >  NC_002947/5381121‑5381413
                                                                                                                                                   |                                                                                                                                                  
gcAACCTGAAGAAAAACCGGCCAAGCCTGGTTGTACAGAAGAGTCCTTCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTg                                                                                                                                                  <  6:263141/149‑1 (MQ=255)
gcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTg                                                                                                                                                  >  7:184770/1‑149 (MQ=255)
      tGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCAGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCgg                                                                                                                                            >  7:216379/1‑149 (MQ=255)
       gaagaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCggg                                                                                                                                           >  5:201026/1‑149 (MQ=255)
       gaagaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCggg                                                                                                                                           >  6:68063/1‑149 (MQ=255)
        aagaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGc                                                                                                                                          <  1:1201/149‑1 (MQ=255)
        aagaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATACAAAAAAGCCCGGACCAAGTCCGGGc                                                                                                                                          <  6:278579/149‑1 (MQ=255)
          gaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATTCAAAAAAGCCCGGACTTGGTCCGGGCtt                                                                                                                                        <  1:193167/149‑1 (MQ=255)
             aaaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttttt                                                                                                                                     >  8:26517/1‑149 (MQ=255)
             aaaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttttt                                                                                                                                     >  6:288514/1‑149 (MQ=255)
              aaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGACAGGACTTGGTCCGGGCtttttt                                                                                                                                    >  5:324085/1‑149 (MQ=255)
                ccGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgt                                                                                                                                  <  8:262716/149‑1 (MQ=255)
                ccGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgt                                                                                                                                  <  4:239909/149‑1 (MQ=255)
                  ggCCAAGCCTGGTTCGACGGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgttt                                                                                                                                <  8:303263/149‑1 (MQ=255)
                    ccAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCCGAAGCGGGGACAAAACAATCAAAAAAACCCCGCAATTGGGCCGGCGttttttttttcg                                                                                                                              >  5:11915/1‑147 (MQ=255)
                        gCCTGGTTCTAACGAAGAGTCCATCCAGCGCGAGATCTACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCt                                                                                                                          <  6:91243/149‑1 (MQ=255)
                             gTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATAc                                                                                                                     <  8:184770/149‑1 (MQ=255)
                              ttCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtttt                                                                                                                                      >  6:137640/1‑131 (MQ=255)
                              ttCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtttt                                                                                                                                      <  5:137640/131‑1 (MQ=255)
                                 gACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCa                                                                                                                 <  6:157427/149‑1 (MQ=255)
                                 gACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtt                                                                                                                                        <  2:275025/126‑1 (MQ=255)
                                 gACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCAGGGCtt                                                                                                                                        >  1:275025/1‑126 (MQ=255)
                                  aaataaGAGTCCTTCAAGAGCGAGATCTACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCa                                                                                                                 <  6:294808/144‑1 (MQ=255)
                                     gAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCa                                                                                                             >  6:28744/1‑149 (MQ=255)
                                           tccatccaGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAAt                                                                                                       >  1:216550/1‑149 (MQ=255)
                                           tccatccaGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAAt                                                                                                       >  5:228316/1‑149 (MQ=255)
                                             catccaGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTg                                                                                                     >  4:132345/1‑149 (MQ=255)
                                              atccaGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGc                                                                                                    >  4:111990/1‑149 (MQ=255)
                                                 caGCGCGATATCGACACCATCGACACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCtt                                                                                                  <  6:246974/148‑1 (MQ=255)
                                                 caGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCtt                                                                                                  <  7:71432/148‑1 (MQ=255)
                                                 caGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTg                                                                                                 <  8:216379/149‑1 (MQ=255)
                                                   gcgcGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGc                                                                                               <  5:68063/149‑1 (MQ=255)
                                                       gagaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGc                                                                                           >  5:193688/1‑149 (MQ=255)
                                                       gagaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGc                                                                                           >  3:140723/1‑149 (MQ=255)
                                                       gagaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCAGTCAATTGATTGGCACGc                                                                                           >  5:189091/1‑149 (MQ=255)
                                                         gaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTg                                                                                         >  7:47001/1‑149 (MQ=255)
                                                            cGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCa                                                                                      >  4:208497/1‑149 (MQ=255)
                                                             gACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAATCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCa                                                                                      <  1:71568/148‑1 (MQ=255)
                                                              acacCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTc                                                                                                          >  2:192537/1‑127 (MQ=255)
                                                              acacCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTATTGTTTGCGCCTTATACGTCAGGCATTc                                                                                                          <  1:192537/127‑1 (MQ=255)
                                                               cacCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGaa                                                                                   >  8:43186/1‑149 (MQ=255)
                                                                    tCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCAt                                                                              >  6:142552/1‑149 (MQ=255)
                                                                       agaCTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGa                                                                           >  6:310909/1‑149 (MQ=255)
                                                                             ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                      <  5:288514/148‑1 (MQ=255)
                                                                             ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                      <  5:28744/148‑1 (MQ=255)
                                                                             ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCg                                                                     >  4:80296/1‑149 (MQ=255)
                                                                               ccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTa                                                                   >  1:86362/1‑149 (MQ=255)
                                                                                        aCCCCGGCACCGCTGGAAACCCGCTAAGCGGGAAGATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCgccgcc                                                          <  6:325176/149‑1 (MQ=255)
                                                                                          cccGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCt                                                         >  4:236445/1‑148 (MQ=255)
                                                                                                  cGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGaa                                                  <  4:99696/147‑1 (MQ=255)
                                                                                                  cGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATa                                                >  3:99696/1‑149 (MQ=255)
                                                                                                   gCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAg                                               >  2:227340/1‑149 (MQ=255)
                                                                                                    cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGaa                                                  >  3:100607/1‑145 (MQ=255)
                                                                                                    cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGaa                                                  <  4:100607/145‑1 (MQ=255)
                                                                                                    cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGc                                              <  5:310909/149‑1 (MQ=255)
                                                                                                       gAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCgg                                            <  3:203130/148‑1 (MQ=255)
                                                                                                       gAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCg                                             <  3:111990/147‑1 (MQ=255)
                                                                                                             cGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGAtt                                                     >  8:54414/1‑133 (MQ=255)
                                                                                                             cGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAAc                                     >  8:174118/1‑149 (MQ=255)
                                                                                                              gCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACg                                    >  2:22533/1‑149 (MQ=255)
                                                                                                                                     aaaaaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc              <  2:216550/149‑1 (MQ=255)
                                                                                                                                     aaaaaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc              >  6:251822/1‑149 (MQ=255)
                                                                                                                                     aaaaaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc              >  2:293736/1‑149 (MQ=255)
                                                                                                                                        aaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTAAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGt            >  5:85167/1‑148 (MQ=255)
                                                                                                                                         aaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTaa          >  5:209242/1‑149 (MQ=255)
                                                                                                                                           gCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTAAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATc        <  6:85167/149‑1 (MQ=255)
                                                                                                                                            cccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCgg               >  7:226884/1‑141 (MQ=255)
                                                                                                                                            cccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCgg               <  8:226884/141‑1 (MQ=255)
                                                                                                                                              cGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCg       >  8:31014/1‑147 (MQ=255)
                                                                                                                                               ggACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAg    <  3:80296/149‑1 (MQ=255)
                                                                                                                                                 aCTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGa  <  6:228316/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGA  >  NC_002947/5381121‑5381413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: