Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,381,267 T→C 9.0% 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.0% 211.0 / 9.3 78intergenic (+32/+33)bamE/PP_4732outer membrane protein assembly factor BamE/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base C (4/3);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

GCGCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGA  >  NC_002947/5381119‑5381413
                                                                                                                                                     |                                                                                                                                                   
gcggAACCTGAAGAAAAACCGGCCAAGCCTGGTTCTACCGAAGACTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACt                                                                                                                                                     <  5:220347/149‑1 (MQ=255)
 cgcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACtt                                                                                                                                                    >  6:175205/1‑149 (MQ=255)
 cgcAACCTGAAGAAAAACCGGCAAAGCCTTTTTAGAAAGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACtt                                                                                                                                                    <  6:295598/149‑1 (MQ=255)
  gcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTg                                                                                                                                                   <  3:80901/149‑1 (MQ=255)
  gcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCACGGACTTg                                                                                                                                                   >  6:259714/1‑149 (MQ=255)
       cTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCg                                                                                                                                              <  1:135282/149‑1 (MQ=255)
         gaagaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCggg                                                                                                                                            >  4:44483/1‑149 (MQ=255)
           agaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGc                                                                                                                                           <  3:147811/148‑1 (MQ=255)
            gaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCACGGACTTGGTCCGGGCtt                                                                                                                                         >  5:69773/1‑149 (MQ=255)
             aaaaaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttt                                                                                                                                        >  4:177572/1‑149 (MQ=255)
              aaaaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCtttt                                                                                                                                       <  5:175205/149‑1 (MQ=255)
                aaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCtttttt                                                                                                                                     >  2:160495/1‑149 (MQ=255)
                   cGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgtt                                                                                                                                  <  1:160495/149‑1 (MQ=255)
                     gCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgtttg                                                                                                                                >  5:209955/1‑149 (MQ=255)
                      ccAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGGCCCGCCCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCCGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTgc                                                                                                                               <  1:227358/149‑1 (MQ=255)
                        aaGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTgcgc                                                                                                                             <  8:240864/149‑1 (MQ=255)
                          gCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCt                                                                                                                           <  3:99235/149‑1 (MQ=255)
                            cTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTta                                                                                                                         >  1:135580/1‑149 (MQ=255)
                            cTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtttttt                                                                                                                                     <  7:120221/137‑1 (MQ=255)
                            cTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtttttt                                                                                                                                     >  8:120221/1‑137 (MQ=255)
                              ggTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGt                                                                                                                                                 >  3:125526/1‑123 (MQ=255)
                              ggTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGt                                                                                                                                                 <  4:125526/123‑1 (MQ=255)
                                  cGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGt                                                                                                                    <  7:120562/148‑1 (MQ=255)
                                    aCCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAg                                                                                                                 >  3:74540/1‑149 (MQ=255)
                                      cGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGc                                                                                                               <  8:194582/149‑1 (MQ=255)
                                         agagTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCAtt                                                                                                            <  3:76638/149‑1 (MQ=255)
                                             tccatccaGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAAt                                                                                                        >  1:304252/1‑149 (MQ=255)
                                                   caGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTg                                                                                                  <  3:148853/149‑1 (MQ=255)
                                                     gcgcGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGc                                                                                                <  3:241547/149‑1 (MQ=255)
                                                       gcgAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCAc                                                                                              >  4:76044/1‑149 (MQ=255)
                                                           gaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTg                                                                                         >  7:53669/1‑149 (MQ=255)
                                                             tCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCTCTGGAAACCCGCTAAGCGTGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGc                                                                                       <  5:259714/149‑1 (MQ=255)
                                                             tCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCTGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGc                                                                                       <  3:96433/149‑1 (MQ=255)
                                                             tCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGc                                                                                       <  2:222231/149‑1 (MQ=255)
                                                               gACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTgg                                                                                        >  5:104271/1‑146 (MQ=255)
                                                               gACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAg                                                                                     <  4:74540/149‑1 (MQ=255)
                                                                     aTCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCa                                                                               >  3:108492/1‑149 (MQ=255)
                                                                      tCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCAt                                                                              >  2:275466/1‑149 (MQ=255)
                                                                             tACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAg                                                                       >  7:45591/1‑149 (MQ=255)
                                                                               ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                      <  7:212798/148‑1 (MQ=255)
                                                                               ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                      >  8:212798/1‑148 (MQ=255)
                                                                               ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                      <  3:239983/148‑1 (MQ=255)
                                                                                cccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCg                                                                     >  8:29536/1‑148 (MQ=255)
                                                                                 ccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTa                                                                   >  2:220613/1‑149 (MQ=255)
                                                                                 ccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTa                                                                   >  7:216149/1‑149 (MQ=255)
                                                                                  cGGTCCCGACCCCGGCACCGCTGGAAACCTGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTAt                                                                  >  1:189280/1‑149 (MQ=255)
                                                                                  cGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGACCGCATCGTCCAGCGTAt                                                                  >  7:130307/1‑149 (MQ=255)
                                                                                    gTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACg                                                                                 >  3:255795/1‑132 (MQ=255)
                                                                                    gTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACg                                                                                 <  4:255795/132‑1 (MQ=255)
                                                                                    gTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTc                                                                <  8:216149/149‑1 (MQ=255)
                                                                                        cGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTcgccg                                                            >  7:35280/1‑149 (MQ=255)
                                                                                             ccGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCATGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGa                                                       >  7:197754/1‑149 (MQ=255)
                                                                                              cGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGAt                                                      >  1:126948/1‑149 (MQ=255)
                                                                                                    cGCTGGAAACCCGCTAAGCGGGAACATCAGATTCAAAAAAGACCGGACTTGGTCCTCGCTTTTATGTTTGCGCCTTATACGTCACGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCATGAtacaat                                                 >  2:268505/1‑143 (MQ=255)
                                                                                                     gCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAg                                               <  6:104271/149‑1 (MQ=255)
                                                                                                     gCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAg                                               <  3:14363/149‑1 (MQ=255)
                                                                                                      cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGc                                              <  1:220613/149‑1 (MQ=255)
                                                                                                             cccGCTAAGCGGGAACATCAGATGCAAAAAAGCTCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGcc                                         >  1:39368/1‑148 (MQ=255)
                                                                                                             cccGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGccc                                        <  1:80709/149‑1 (MQ=255)
                                                                                                                    aGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTg                                  <  2:304252/147‑1 (MQ=255)
                                                                                                                       gggAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGgcgcg                             >  6:144187/1‑149 (MQ=255)
                                                                                                                          aaCATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGcac                          <  3:177572/149‑1 (MQ=255)
                                                                                                                            cATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACaa                        <  5:271958/149‑1 (MQ=255)
                                                                                                                              tCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATc                      <  8:45591/149‑1 (MQ=255)
                                                                                                                                 gATGCAAAAAAGCCCGGACTTGGTCCGGGGGTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGATCGCATAGACCAGCGTATTCGCCGACTGATTGAATAGCGGGCCCAACGTGGCGAGCACAAACGGc                   >  3:17738/1‑149 (MQ=255)
                                                                                                                                 gATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGc                   <  2:126948/149‑1 (MQ=255)
                                                                                                                                     cAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCgg               >  3:111244/1‑149 (MQ=255)
                                                                                                                                      aaaaaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAg                                                                       <  4:91112/92‑1 (MQ=255)
                                                                                                                                      aaaaaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGc                                              <  1:158030/117‑1 (MQ=255)
                                                                                                                                      aaaaaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGc                                              >  2:158030/1‑117 (MQ=255)
                                                                                                                                       aaaaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTAGGCACGCTGGCAGAACGCATCGAACAGCGAATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc              >  6:3207/1‑148 (MQ=255)
                                                                                                                                          aaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTaa          <  8:53669/149‑1 (MQ=255)
                                                                                                                                             cccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATc        >  4:212072/1‑148 (MQ=255)
                                                                                                                                             cccGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATc        >  4:165618/1‑148 (MQ=255)
                                                                                                                                               ccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc              <  7:78142/141‑1 (MQ=255)
                                                                                                                                               ccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc              >  8:78142/1‑141 (MQ=255)
                                                                                                                                                cGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGGCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGaa     <  2:135580/149‑1 (MQ=255)
                                                                                                                                                cGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGaa     <  6:69773/149‑1 (MQ=255)
                                                                                                                                                cGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGaa     <  2:39368/149‑1 (MQ=255)
                                                                                                                                                 ggACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAg    <  6:209955/149‑1 (MQ=255)
                                                                                                                                                 ggACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAg    <  3:44483/149‑1 (MQ=255)
                                                                                                                                                     ttGGTCCGGGCTTTTTTGAATGCGCCTTATACGTCAGCAAATCAATTACTTGGCACGCTGGCAGAACGCAACGACAAGCTTATTCTCAGACTGAATCAATAGCGGGCCCAACGTGGCGCGCCCAATCGGCCCGCCGAAATCGAAAGa  >  6:97473/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGA  >  NC_002947/5381119‑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: