Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,381,269 G→A 10.0% intergenic (+34/+31) 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,2690GA10.0% 186.7 / 8.2 69intergenic (+34/+31)bamE/PP_4732outer membrane protein assembly factor BamE/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base A (4/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAG  >  NC_002947/5381121‑5381416
                                                                                                                                                     |                                                                                                                                                   
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                                                                                                                                         <  5:137640/131‑1 (MQ=255)
                              ttCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtttt                                                                                                                                         >  6:137640/1‑131 (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                                                                                              >  3:140723/1‑149 (MQ=255)
                                                       gagaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGc                                                                                              >  5:193688/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                 >  6:251822/1‑149 (MQ=255)
                                                                                                                                     aaaaaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc                 <  2:216550/149‑1 (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)
                                                                                                                                                     ggTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAg  <  4:140723/148‑1 (MQ=255)
                                                                                                                                                     ggTACGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGACTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAg  <  6:189091/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAG  >  NC_002947/5381121‑5381416

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: