Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,955,821 A→C 7.4% intergenic (+47/+157) elbB → / ← PP_5221 putative hydrolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,955,8210AC 7.4% 191.5 / 3.5 65intergenic (+47/+157)elbB/PP_5221putative hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/24);  new base C (2/3);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACA  >  NC_002947/5955674‑5955969
                                                                                                                                                   |                                                                                                                                                    
tCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAg                                                                                                                                                     >  8:117900/1‑149 (MQ=255)
  gTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAgc                                                                                                                                                    <  6:292382/148‑1 (MQ=255)
   tGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAgcg                                                                                                                                                   >  3:125298/1‑148 (MQ=255)
         ccccGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCtt                                                                                                                                            >  8:274118/1‑149 (MQ=255)
             ggCGTACATGGAGGCCAAGTCGATCAGCGAGGCGTCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAg                                                                                                                                        <  4:125298/149‑1 (MQ=255)
             ggCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAg                                                                                                                                        <  7:5443/149‑1 (MQ=255)
                        aGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTa                                                                                                                             >  6:250867/1‑149 (MQ=255)
                             aaGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCtgt                                                                                                                        >  1:209281/1‑149 (MQ=37)
                                  gATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTggga                                                                                                                    <  2:209281/148‑2 (MQ=37)
                                   aTCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCAcgc                                                                                                                                 >  4:254659/1‑134 (MQ=255)
                                   aTCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCAcgc                                                                                                                                 <  3:254659/134‑1 (MQ=255)
                                   aTCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGc                                                                                                                   >  1:252297/1‑148 (MQ=255)
                                     cAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCAcgcg                                                                                                                                <  5:323113/133‑1 (MQ=255)
                                     cAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCAcgcg                                                                                                                                >  6:323113/1‑133 (MQ=255)
                                      aGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACt                                                                                                               >  8:239427/1‑149 (MQ=255)
                                      aGCGAGGCGGCGGGCGGGAGCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACt                                                                                                               <  3:177193/149‑1 (MQ=255)
                                       gCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTg                                                                                                              >  5:133070/1‑149 (MQ=255)
                                          aGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTgcgc                                                                                                           <  5:290722/149‑1 (MQ=255)
                                              ggcgggcgGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATa                                                                                                       <  6:219169/149‑1 (MQ=255)
                                              ggcgggcgGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGCGCAACTGCGCGATa                                                                                                       <  7:274118/149‑1 (MQ=255)
                                                  ggcgggATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTCCCTGTGGGGCAACTGCGCGATAGCCt                                                                                                   >  2:102787/1‑149 (MQ=255)
                                                            aCAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCaa                                                                                          >  7:108173/1‑148 (MQ=255)
                                                            aCAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAAc                                                                                         >  3:224502/1‑149 (MQ=255)
                                                            aCAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAAc                                                                                         >  7:17658/1‑149 (MQ=255)
                                                             cAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCtgtg                                                                                                                       >  3:85226/1‑118 (MQ=255)
                                                             cAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCtgtg                                                                                                                       <  4:85226/118‑1 (MQ=255)
                                                             cAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACt                                                                                        <  6:133070/149‑1 (MQ=255)
                                                              aaGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTg                                                                                       >  7:20320/1‑149 (MQ=255)
                                                                 cTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTgcgc                                                                                    <  2:252297/149‑1 (MQ=255)
                                                                   ggtggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGa                                                                                  >  7:51891/1‑149 (MQ=255)
                                                                   ggtggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGa                                                                                  >  8:270521/1‑149 (MQ=255)
                                                                     tggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATc                                                                                <  7:270521/149‑1 (MQ=255)
                                                                         ccGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCtt                                                                            >  8:256379/1‑149 (MQ=255)
                                                                         ccGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCtt                                                                            <  3:92153/149‑1 (MQ=255)
                                                                         ccGAGTGCTGCAGCTGACGAACGAGTAACTTACTGCGTTGGCTTCTTCTCGTGCATGCCCGCTCCCACAGGTCAATCGCACGCTTCAACCCCACCCGGTACCTGTGCGGAAACTGAGCTATAGCCTGTGGGGCAACTGAGCGAACCCtt                                                                            >  2:166110/1‑149 (MQ=255)
                                                                          cGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTg                                                                           <  8:51891/149‑1 (MQ=255)
                                                                              ttcTGGAGCTGACGCACGAGTAAATTACTGCGTTGGCTTCTTCGGGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTa                                                                         <  1:166110/145‑1 (MQ=255)
                                                                                   gAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCgg                                                                  <  8:17658/149‑1 (MQ=255)
                                                                                        gACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTa                                                             >  7:74916/1‑149 (MQ=255)
                                                                                         aCGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTAc                                                            <  7:117900/149‑1 (MQ=255)
                                                                                               gAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCcgcg                                                       >  8:75737/1‑148 (MQ=255)
                                                                                                      ttACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCgg                                               >  6:6577/1‑149 (MQ=255)
                                                                                                        aCTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATc                                                                                >  4:176500/1‑114 (MQ=255)
                                                                                                        aCTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATc                                                                                <  3:176500/114‑1 (MQ=255)
                                                                                                            cGTTGGCTTGCTCCCGGGCCTGCCCGGGCCCCCCTGTCCCTCGCCTGCTTCCACGCCCCGGGGTAACCGTGNGGCCACCGCGCGCTCGCCGGGGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGc                                         <  1:102787/149‑1 (MQ=16)
                                                                                                            cGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGctctgc                                         <  4:224502/149‑6 (MQ=255)
                                                                                                            cGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGc                                         <  7:239427/149‑1 (MQ=255)
                                                                                                                 gCTTCTTAGCGGGCATGCCCGCTCCCACAGGTCCAGCGAAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGAAGGAGCGGGTTTACCCGCCAACACCGGCAGAGCCggag                                    >  4:185333/1‑147 (MQ=255)
                                                                                                                   ttcttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGc                                   >  6:197819/1‑148 (MQ=255)
                                                                                                                   ttcttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGc                                   >  2:270598/1‑148 (MQ=255)
                                                                                                                   ttcttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGc                                   >  3:157982/1‑148 (MQ=255)
                                                                                                                        cgcgGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCTCTGCCGGTgttcgc                               <  7:256379/147‑6 (MQ=255)
                                                                                                                              cATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCgcgc                       >  1:205157/1‑149 (MQ=255)
                                                                                                                                 gCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAg                    >  5:252053/1‑149 (MQ=255)
                                                                                                                                 gCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAg                    >  4:198101/1‑149 (MQ=255)
                                                                                                                                 gCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAg                    >  4:202155/1‑149 (MQ=255)
                                                                                                                                   ccGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCCGCGCCAGcc                  >  6:139160/1‑149 (MQ=255)
                                                                                                                                   ccGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCGGTGGGGCAACTGCGCGATAGCCTGTGGCGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCAGCGAACACCGCCAGCGCCGGTGCCAACCACCGAGCCAccc                  >  1:228582/1‑149 (MQ=38)
                                                                                                                                         ccACAGGTCCAGGGCAGGCTTCTAGGCCACGCGGTAGCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGc             <  8:108173/148‑1 (MQ=255)
                                                                                                                                          cacaGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGc                   <  3:291057/141‑1 (MQ=25)
                                                                                                                                          cacaGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGc                   >  4:291057/1‑141 (MQ=25)
                                                                                                                                               gTCCTGCGCAGGCTTCAAGGCCTCGCGGTACCTGTGGGGCAAGTGCGGGATAGCGTGTGGGGCAACTGCGCGATCGGTTGTAGGAGGGGGTTTACCGGCGATCAGGGGCTGAGCCGGTGCCATCCAGCGCGCCAGCCTGGGCTGTCGcc      <  2:205157/149‑1 (MQ=255)
                                                                                                                                               gTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGcc      >  6:138911/1‑149 (MQ=255)
                                                                                                                                                tCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGt                                     >  2:296873/1‑117 (MQ=255)
                                                                                                                                                tCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGt                                     <  1:296873/117‑1 (MQ=255)
                                                                                                                                                   aGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACa  <  1:270598/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACA  >  NC_002947/5955674‑5955969

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: