Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,798,379	T→C	7.7%	D495G (GAC→GGC) ‡	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,798,379	0	T	C	7.7%	215.7 / 7.2	79	D495G (GAC→GGC) ‡	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base C (3/3);  total (36/43)
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.25e-01

GCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCACATAATT  >  NC_002947/5798232‑5798527                                                                                                                                                    |                                                                                                                                                     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                                                                                                                               >  2:256639/1‑149 (MQ=255)                       tGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                                                                                                                               <  5:33492/149‑1 (MQ=255)                       tGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                                                                                                                               <  4:251385/149‑1 (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                                                                                      <  2:394751/148‑1 (MQ=255)                                                                 gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                                      <  8:264979/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                                                      <  5:145619/138‑1 (MQ=255)                                                                                                           ggACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGTAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGt                                                      >  6:145619/1‑138 (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:139060/1‑148 (MQ=255)                                                                                                                            ttCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTgg                           >  4:205901/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)                                                                                                                                                    tCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCACATAAtt  >  7:94495/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCTTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCACATAATT  >  NC_002947/5798232‑5798527

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