Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,798,375	T→G	7.2%	A496A (GCA→GCC) ‡	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,798,375	0	T	G	7.2%	247.5 / 7.3	80	A496A (GCA→GCC) ‡	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/41);  new base G (3/3);  total (36/44)
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.73e-01

TTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCA  >  NC_002947/5798229‑5798520                                                                                                                                                   |                                                                                                                                                  ttCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                                                                                                                 <  3:151921/149‑1 (MQ=255)    gCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTc                                                                                                                                              >  5:132214/1‑149 (MQ=255)             cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTTGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCGACGCAGCGGTTg                                                                                                                                       <  5:243540/147‑1 (MQ=255)             cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                                                                                                     >  6:66464/1‑149 (MQ=255)             cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCGACGCAGCGGTTgcg                                                                                                                                     <  7:50686/149‑3 (MQ=255)             cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCGACGCAGCGGTTgcg                                                                                                                                     <  1:426946/149‑3 (MQ=255)                 gcggcgGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGaaa                                                                                                                                 <  2:327996/149‑1 (MQ=255)                  cggcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAAc                                                                                                                                <  4:227975/149‑1 (MQ=255)                   ggcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgcg                                                                                                                                           >  4:122926/1‑137 (MQ=255)                   ggcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgcg                                                                                                                                           <  3:122926/137‑1 (MQ=255)                    gcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTg                                                                                                                              <  3:30670/149‑1 (MQ=255)                    gcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTg                                                                                                                              <  8:232281/149‑1 (MQ=255)                        ggTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGt                                                                                                                          <  3:397226/149‑1 (MQ=255)                          tGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                                                                                                                        <  5:33492/149‑1 (MQ=255)                          tGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                                                                                                                        <  4:251385/149‑1 (MQ=255)                          tGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                                                                                                                        >  2:256639/1‑149 (MQ=255)                            aTGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGt                                                                                                                      >  5:220313/1‑149 (MQ=255)                             tGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTc                                                                                                                     <  6:132214/149‑1 (MQ=255)                               ccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                                                                                                                   >  4:427981/1‑149 (MQ=255)                                 gccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCacga                                                                                                                 >  8:245661/1‑149 (MQ=255)                                          gTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACagag                                                                                                        >  6:46419/1‑149 (MQ=255)                                               aGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCtt                                                                                                    <  5:66464/148‑1 (MQ=255)                                               aGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTg                                                                                                   <  2:249433/149‑1 (MQ=255)                                               aGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTg                                                                                                   >  5:167172/1‑149 (MQ=255)                                                  gggCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGACTTGccc                                                                                                >  6:409224/1‑149 (MQ=255)                                                    gCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTg                                                                                              <  1:117688/149‑1 (MQ=255)                                                     ccGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGc                                                                                             <  7:14221/149‑1 (MQ=255)                                                     acggggCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTTCGGCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGc                                                                                             <  4:22907/145‑1 (MQ=255)                                                      cGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTg                                                                                              >  1:406603/1‑147 (MQ=255)                                                       gTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAg                                                                                           <  3:95056/149‑1 (MQ=255)                                                          gCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACACCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGtt                                                                                        >  6:406356/1‑149 (MQ=255)                                                              gCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTc                                                                                    >  2:5011/1‑149 (MQ=255)                                                                    gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                                 <  2:159232/146‑1 (MQ=255)                                                                    gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                               <  8:264979/148‑1 (MQ=255)                                                                    gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                               <  2:394751/148‑1 (MQ=255)                                                                    gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGa                                                                                <  6:236081/147‑1 (MQ=255)                                                                     ttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGt                                                                             <  1:5011/149‑1 (MQ=255)                                                                      tgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAg                                                                              >  8:123090/1‑147 (MQ=255)                                                                        gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGc                                                                          >  3:318074/1‑149 (MQ=255)                                                                          cacTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGa                                                                        <  7:262214/147‑1 (MQ=255)                                                                          cGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGa                                                                        >  2:231164/1‑149 (MQ=255)                                                                            cTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATa                                                                      <  1:231164/149‑1 (MQ=255)                                                                             tGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgcg                                                                                                                                           >  8:194851/1‑79 (MQ=255)                                                                               ccGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGACTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAg                                                                     <  5:409224/147‑1 (MQ=255)                                                                                    ggTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCAc                                                              <  6:220313/149‑1 (MQ=255)                                                                                    ggTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCAc                                                              >  6:220697/1‑149 (MQ=255)                                                                                        cGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTg                                                          <  1:256639/149‑1 (MQ=255)                                                                                        cGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTg                                                          >  4:31485/1‑149 (MQ=255)                                                                                          tGGCCAAGGTTTACCAGCCTGGAGCTTCGCTTGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGCGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTggg                                                        <  5:406356/149‑1 (MQ=255)                                                                                             ccAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCt                                                     >  3:81759/1‑149 (MQ=255)                                                                                                      aCCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTgg                                                         >  7:132626/1‑136 (MQ=255)                                                                                                      aCCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTgg                                                         <  8:132626/136‑1 (MQ=255)                                                                                                       ccAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGc                                           <  3:31485/149‑1 (MQ=255)                                                                                                          gCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCTGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTg                                        >  6:303133/1‑149 (MQ=255)                                                                                                           ccTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTgg                                       <  8:106966/149‑1 (MQ=255)                                                                                                              ggACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGTAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGt                                               >  6:145619/1‑138 (MQ=255)                                                                                                              ggACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGTAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGt                                               <  5:145619/138‑1 (MQ=255)                                                                                                              ggACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGc                                                                 >  7:88629/1‑120 (MQ=255)                                                                                                              ggACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGc                                                                 <  8:88629/120‑1 (MQ=255)                                                                                                              ggACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGcc                                    >  7:375839/1‑149 (MQ=255)                                                                                                                 ccTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGAc                                 >  7:118868/1‑149 (MQ=255)                                                                                                                  cTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACg                                >  5:427593/1‑149 (MQ=255)                                                                                                                          gCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGTCCGACGCCTCGGAg                        >  4:285668/1‑149 (MQ=255)                                                                                                                            aaGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCt                      <  7:245661/149‑1 (MQ=255)                                                                                                                             aGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTg                     <  6:167172/149‑1 (MQ=255)                                                                                                                               ttCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTgg                    >  4:205901/1‑148 (MQ=255)                                                                                                                               ttCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTgg                    >  4:139060/1‑148 (MQ=255)                                                                                                                                    gATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCt                      <  3:429299/141‑1 (MQ=255)                                                                                                                                    gATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCt                      >  4:429299/1‑141 (MQ=255)                                                                                                                                    gATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTAcc                <  4:318074/147‑1 (MQ=255)                                                                                                                                    gATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTAc                 <  8:375839/146‑1 (MQ=255)                                                                                                                                    gATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCa               <  7:123090/148‑1 (MQ=255)                                                                                                                                     aTTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAg             <  5:46419/149‑1 (MQ=255)                                                                                                                                     aTTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAg             >  1:319431/1‑149 (MQ=255)                                                                                                                                     aTTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGAAAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCaa              <  5:303133/148‑1 (MQ=255)                                                                                                                                           ccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTg       >  5:154772/2‑149 (MQ=255)                                                                                                                                           ccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTATTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTg       >  7:389189/2‑149 (MQ=255)                                                                                                                                            cAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGc      <  5:220697/149‑1 (MQ=255)                                                                                                                                            cAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTg       >  1:394456/1‑148 (MQ=255)                                                                                                                                                cGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGcacca  >  5:147389/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCA  >  NC_002947/5798229‑5798520

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑