Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,955,819 C→A 6.2% intergenic (+45/+159) 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,8190CA 6.2% 208.2 / 3.6 65intergenic (+45/+159)elbB/PP_5221putative hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base A (2/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

TCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCC  >  NC_002947/5955674‑5955965
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                                   <  4:85226/118‑1 (MQ=255)
                                                             cAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCtgtg                                                                                                                   >  3:85226/1‑118 (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                                                                        <  3:92153/149‑1 (MQ=255)
                                                                         ccGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCtt                                                                        >  8:256379/1‑149 (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                                                                            <  3:176500/114‑1 (MQ=255)
                                                                                                        aCTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATc                                                                            >  4:176500/1‑114 (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                               >  3:157982/1‑148 (MQ=255)
                                                                                                                   ttcttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGc                               >  6:197819/1‑148 (MQ=255)
                                                                                                                   ttcttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGc                               >  2:270598/1‑148 (MQ=255)
                                                                                                                        cgcgGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCTCTGCCGGTgttcgc                           <  7:256379/147‑6 (MQ=255)
                                                                                                                              cATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCgcgc                   >  1:205157/1‑149 (MQ=255)
                                                                                                                                 gCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAg                >  4:198101/1‑149 (MQ=255)
                                                                                                                                 gCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAg                >  4:202155/1‑149 (MQ=255)
                                                                                                                                 gCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAg                >  5:252053/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)
                                                                                                                                                 |                                                                                                                                                  
TCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCC  >  NC_002947/5955674‑5955965

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: