Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,538,532	T→G	5.4%	intergenic (+22/+99)	PP_4871 → / ← PP_4872	lysine decarboxylase family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,538,532	0	T	G	5.4%	210.0 / 2.0	67	intergenic (+22/+99)	PP_4871/PP_4872	lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base G (2/2);  total (33/34)
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.66e-01

AGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCC  >  NC_002947/5538384‑5538679                                                                                                                                                      |                                                                                                                                                      agGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCt                                                                                                                                                       >  2:103360/1‑149 (MQ=255) agGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCt                                                                                                                                                       >  8:140745/1‑149 (MQ=255)   gTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTc                                                                                                                                                     <  1:108784/149‑1 (MQ=255)    tttCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTcg                                                                                                                                                    >  3:77239/1‑149 (MQ=255)        gTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCgg                                                                                                                                                 <  2:248146/148‑1 (MQ=255)        gTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGt                                                                                                                                                <  1:114181/149‑1 (MQ=255)              ccGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTc                                                                                                                                               >  1:367869/1‑144 (MQ=255)              ccGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTc                                                                                                                                               <  2:367869/144‑1 (MQ=255)              ccGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTc                                                                                                                                          <  1:68215/149‑1 (MQ=255)               cGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGGCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCGAAGCGGTcccaaaa                                                                                                                                         <  5:265812/149‑6 (MQ=255)                  aaCACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGa                                                                                                                                      <  1:103360/149‑1 (MQ=255)                        gggCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGGCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAg                                                                                                                                <  4:111817/149‑1 (MQ=255)                         ggCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAgg                                                                                                                               <  4:77239/149‑1 (MQ=255)                         ggCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAgg                                                                                                                               >  2:254695/1‑149 (MQ=255)                          gCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGATGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGc                                                                                                                              <  5:335465/149‑1 (MQ=255)                            gATGTTGCTGCTGGGGCAGCAGCCGGACGAGGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGcc                                                                                                                             >  1:240748/1‑149 (MQ=255)                            gATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCg                                                                                                                            <  1:39787/149‑1 (MQ=255)                                ttGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑cc                                                                                                                      <  3:113551/149‑1 (MQ=255)                                aagcgactgGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑cc                                                                                                                      <  7:346619/143‑1 (MQ=255)                                  gctgctGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑cccc                                                                                                                    >  7:340000/1‑148 (MQ=255)                                  gctgctGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCca                                                                                                                    >  5:70142/1‑149 (MQ=255)                                    tgctgGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGCCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCcac                                                                                                                   >  8:38209/1‑148 (MQ=255)                                         gggCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCGTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAgg                                                                                                              <  4:268177/148‑1 (MQ=255)                                           gcagcaGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCT‑AGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAAt                                                                                                          >  6:188524/1‑149 (MQ=255)                                            cagcagCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAAt                                                                                                          >  6:40805/1‑149 (MQ=255)                                                  ccGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATgggg                                                                                                      >  3:311576/1‑147 (MQ=255)                                                     gACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCg                                                                                                 <  2:145411/149‑1 (MQ=255)                                                      aCGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCCCATTGCCTTTGTGGGAGCgg                                                                                                <  2:240748/149‑1 (MQ=38)                                                       cGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCCCATTGCCTTTGTGGGAGCGGc                                                                                               <  8:340000/149‑1 (MQ=38)                                                        gA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCACAAAGGCAATGGGGGAGCgg                                                                                                <  2:379411/149‑1 (MQ=255)                                                        gA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCCCATTGCCTTTGTGGGAGCGGcc                                                                                              <  7:140745/149‑1 (MQ=38)                                                           gctgctGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                            <  5:40805/149‑1 (MQ=255)                                                                 ggAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGcc                                                                                              <  4:344912/141‑1 (MQ=255)                                                                 ggAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGcc                                                                                              >  3:344912/1‑141 (MQ=255)                                                                 ggAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcg                                                                                      >  6:197088/1‑149 (MQ=255)                                                                  gAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTATCGCGACACAAGGCCGCT‑‑CCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTcgc                                                                                     <  8:6251/149‑1 (MQ=37)                                                                                 tttGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAagc                                                                      >  8:235250/1‑149 (MQ=255)                                                                                   tgtgGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAagcag                                                                    >  5:152136/1‑149 (MQ=255)                                                                                               cgcgCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTc                                                                                                                                          <  3:74077/69‑1 (MQ=255)                                                                                               cgcgCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTc                                                                                                                                          >  4:74077/1‑69 (MQ=255)                                                                                                gcgcCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGa                                                       >  6:334656/1‑149 (MQ=255)                                                                                                   ccTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGt                                                          <  3:252088/143‑1 (MQ=255)                                                                                                   ccTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGt                                                          >  4:252088/1‑143 (MQ=255)                                                                                                    cTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAg                                                   >  4:131183/1‑149 (MQ=255)                                                                                                    cTAAGTGGGGCGACAAGCCGCCTGCCTGAAGCTGAGACTTTTGGGACCGCGTCGCGGTCCTTTCGCGACCCAAGGCCGCT‑‑CCCACACAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATc                                                     <  6:70142/147‑1 (MQ=38)                                                                                                    cTAAGTGGGCCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAg                                                   <  7:38209/149‑1 (MQ=255)                                                                                                      aaGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGt                                                  >  5:171542/1‑148 (MQ=255)                                                                                                                   agcagcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGAc                                    <  1:254695/149‑1 (MQ=255)                                                                                                                      agcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAgg                                 <  5:188524/149‑1 (MQ=255)                                                                                                                      agcCTGACGGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAgg                                 <  2:281999/149‑1 (MQ=255)                                                                                                                       gcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGc                                >  8:157180/1‑149 (MQ=255)                                                                                                                        ccTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAATCACACCCAATAGTCGATCAGTTAAAGCGGAATGa                                     >  1:372153/1‑143 (MQ=255)                                                                                                                           gactgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTg                              >  3:121471/1‑147 (MQ=255)                                                                                                                                aaGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGc                       <  7:349677/149‑1 (MQ=255)                                                                                                                                   cTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttc                    <  5:174886/149‑1 (MQ=255)                                                                                                                                   cTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttc                    >  4:350989/1‑149 (MQ=255)                                                                                                                                    tGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGAAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttct                   >  5:325389/1‑149 (MQ=255)                                                                                                                                         cTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttgtt               <  7:235250/148‑1 (MQ=255)                                                                                                                                          ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAGGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGc             >  3:129431/1‑149 (MQ=255)                                                                                                                                           gaaaggACCGCGAAGCGGTCCTTTCGCGACACAAGGCCGCA‑‑CCCACAAAGGCAATGGGGGAGCGGCCGTGTGACGCGATGGGCTGCAAAGCACACCCAATAGTCGATCAGTTAAACCGGGATGACAGGCTGAAGCGGCGTCatgttggc            >  4:390824/5‑147 (MQ=255)                                                                                                                                             aaggACCGCGAAGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCt          >  5:354213/3‑149 (MQ=255)                                                                                                                                               ggACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTg        >  5:109512/1‑149 (MQ=255)                                                                                                                                               ggACCGCTTCGCGGTCCTTTCGCGACACAAGGCCACT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGCCGCGGCTTCTTGTTGCCCTTg        >  2:324634/1‑149 (MQ=255)                                                                                                                                                   cGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCgg                                          >  6:167210/1‑111 (MQ=255)                                                                                                                                                   cGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCgg                                          <  7:274718/111‑1 (MQ=255)                                                                                                                                                   cGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCgg                                          <  5:167210/111‑1 (MQ=255)                                                                                                                                                   cGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCgg                                          >  8:274718/1‑111 (MQ=255)                                                                                                                                                   cGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCAc    >  4:347547/1‑149 (MQ=255)                                                                                                                                                     cTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGCTGGGCTGCAAAGCAGCCCCAATAGTGGCTCAGTTAACGCGGCATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGGCAccc  >  7:150210/1‑149 (MQ=39)                                                                                                                                                      |                                                                                                                                                      AGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGA‑GCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCT‑‑CCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCC  >  NC_002947/5538384‑5538679

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: ‑