Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,770,444	A→T	5.9%	M143K (ATG→AAG)	PP_5060 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,770,444	0	A	T	5.9%	208.7 / 2.5	68	M143K (ATG→AAG)	PP_5060	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base T (2/2);  total (36/32)
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.93e-01

GCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGTG  >  NC_002947/5770298‑5770591                                                                                                                                                   |                                                                                                                                                       gcTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCAtt                                                                                                                                                      <  4:289885/149‑1 (MQ=255)     gCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTgg                                                                                                                                                    <  6:108863/147‑1 (MQ=255)      cGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgcag                                                                                                                                                 <  6:135967/149‑1 (MQ=255)       gttgttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcact                                                                                                                                                >  6:31496/1‑147 (MQ=255)       gttgttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgcagc                                                                                                                                                <  7:282071/149‑1 (MQ=255)        ttgttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgcagca                                                                                                                                               <  1:168952/149‑1 (MQ=255)           ttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcca                                                                                                                                            <  5:31496/149‑7 (MQ=255)           ttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAAcg                                                                                                                                            <  6:256278/149‑1 (MQ=255)               tCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGtt                                                                                                                                        <  3:146839/149‑1 (MQ=255)                  tCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAAcgc                                                                                                                                           >  8:43211/1‑143 (MQ=255)                  tCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAAcgc                                                                                                                                           <  7:43211/143‑1 (MQ=255)                    cAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcc                                                                                                                                             <  8:237214/139‑6 (MQ=255)                    cAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTACAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGAAGCATTGGCAGAAACGAGTTCAGGc                                                                                                                                   >  3:2844/1‑149 (MQ=255)                    cAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTAGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcc                                                                                                                                             >  7:237214/1‑134 (MQ=255)                     aGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGc                                                                                                                                   >  8:35310/1‑148 (MQ=255)                       cgcTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAt                                                                                                                                <  4:141880/149‑1 (MQ=255)                                cTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAt                                                                                                                                >  1:215281/1‑140 (MQ=255)                                cTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAt                                                                                                                                <  2:215281/140‑1 (MQ=255)                                 ttCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGgcg                                                                                                                     >  1:16130/1‑149 (MQ=255)                                 ttCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGgcg                                                                                                                     >  7:113989/1‑149 (MQ=255)                                    aGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCAATCG‑GTCAGgcgcgc                                                                                                                  >  2:165585/1‑149 (MQ=255)                                     ggCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCACCGCGTTCCGGCCATCG‑GTCAGGcgcccg                                                                                                                 >  8:120153/1‑149 (MQ=255)                                         aTCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCt                                                                                                             <  4:211562/149‑1 (MQ=255)                                            gAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCgg                                                                                                          >  2:70554/1‑148 (MQ=255)                                            gAGCAGGGTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGGTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAGCG‑GTCAGGCGCGCGGGCTCg                                                                                                           <  7:120153/148‑1 (MQ=255)                                                  gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCGGTCGCCGTTCCAATAGTCGATCTGCGCTGGCTGGATGTGCCAAAGCAGCATTGGCAGCAGCGCGTTCAGGCCTTCG‑GTCGGGCGCGCGGGCTCGATCatat                                                                                                    >  8:266705/1‑149 (MQ=255)                                                  gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAgg                                                                                                                                    >  8:43487/1‑118 (MQ=255)                                                  gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAt                                                                                                    >  2:81926/1‑149 (MQ=255)                                                      gcgcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGcc                                                                                              <  8:113989/149‑1 (MQ=255)                                                        gcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGtt                                                                                                                                        >  2:86048/1‑108 (MQ=255)                                                        gcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGtt                                                                                                                                        <  1:86048/108‑1 (MQ=255)                                                            aacaacTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGc                                                                                        >  6:141840/1‑149 (MQ=255)                                                                 ctgctgCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGttt                                                                                   <  1:81926/149‑1 (MQ=255)                                                                      gcAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAg                                                                              >  2:249496/1‑149 (MQ=255)                                                                       cAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGaa                                                                               <  5:141840/147‑1 (MQ=255)                                                                          ccACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGATTCGAAGctc                                                                           >  7:81562/1‑148 (MQ=255)                                                                              aGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCAAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGctcgactc                                                                      >  8:211703/1‑149 (MQ=255)                                                                                   ccGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATgg                                                                 >  1:152707/1‑149 (MQ=255)                                                                                     ggTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCACCCGTTTCGAAGCTCGACTCGATGGCg                                                               <  1:165585/149‑1 (MQ=255)                                                                                          aTCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGaaaaa                                                          <  7:158997/149‑1 (MQ=255)                                                                                             gTGCCGTCGCCGTTCCAATAGTCGCTCTGCGCTTGCTGGAAGTGCCATAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAgc                                                        <  4:28591/148‑1 (MQ=255)                                                                                              tGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAgcgc                                                      >  4:266164/1‑149 (MQ=255)                                                                                                   tCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCgggg                                                 >  6:228048/1‑149 (MQ=255)                                                                                                      ccGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                              <  3:240964/149‑1 (MQ=255)                                                                                                      ccGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                              <  2:16130/149‑1 (MQ=255)                                                                                                      ccGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                              >  4:240964/1‑149 (MQ=255)                                                                                                       cGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCggg                                             <  3:128720/149‑1 (MQ=255)                                                                                                         ttCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGtt                                           <  5:228048/149‑1 (MQ=255)                                                                                                           ccAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGtt                                           >  3:187974/1‑147 (MQ=255)                                                                                                           ccAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGtg                                           >  3:194736/1‑146 (MQ=255)                                                                                                             aaTAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAgg                                       >  6:224377/1‑149 (MQ=255)                                                                                                              aTAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGc                                      <  4:187974/149‑1 (MQ=255)                                                                                                                     aTCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGAATCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAACGCGCCGGGGCGGGTTCAGGCTAAGGc                                >  1:238670/1‑149 (MQ=255)                                                                                                                      tCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCgg                              >  7:28898/1‑149 (MQ=255)                                                                                                                          ggcTTTCTGGAAGTGGCAAATCATCATTTGGATCAAAGCGTTCATGCCATTG‑GTGATGGGCGGGGGCTCGATCAAATTGCTCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTa                            <  1:249496/148‑1 (MQ=255)                                                                                                                                 tGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCg                                                                                 >  4:224481/1‑87 (MQ=255)                                                                                                                                  ggAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTcgcccgcc                  >  4:83185/1‑149 (MQ=255)                                                                                                                                   gAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCccgccg                 >  1:123511/1‑149 (MQ=255)                                                                                                                                    aaGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATc                                                                                                        <  1:156071/63‑1 (MQ=255)                                                                                                                                    aaGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATc                                                                                                        >  2:156071/1‑63 (MQ=255)                                                                                                                                     aGTGCAAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCgg                <  4:2844/148‑1 (MQ=255)                                                                                                                                       ggCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGc             <  3:1396/148‑1 (MQ=255)                                                                                                                                         ccAACGCAGCATTCGCAGCAACGCGTTCAGGCCATCG‑TTCAGGCGCGCGGGCACGATCAAAT‑‑CGCCCTCAGCCCGTTCGAAGCTCGACTCGATGGCGAAAAAGCCCCGTGGCGGGTTCAGGCTAAGGCGGTAGACGCCCGCCGGGGc             >  3:243899/1‑147 (MQ=255)                                                                                                                                         ccAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGc             >  4:219192/1‑147 (MQ=255)                                                                                                                                          cAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAg          <  2:152707/149‑1 (MQ=255)                                                                                                                                          cAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAg          <  7:35310/149‑1 (MQ=255)                                                                                                                                          cAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAg          >  2:187738/1‑149 (MQ=255)                                                                                                                                           aaaGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAgg         >  2:47060/1‑149 (MQ=255)                                                                                                                                               cagcaTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAACT‑‑CGCCCTCAGCCGTATCGACGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCa     <  7:266705/149‑1 (MQ=255)                                                                                                                                                 gcaTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGt   >  5:136117/1‑149 (MQ=255)                                                                                                                                                 gcaTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGt   >  1:189749/1‑149 (MQ=255)                                                                                                                                                 gcaTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGt   <  5:224377/149‑1 (MQ=255)                                                                                                                                                 gcaATGGCAGCAACGAGTTCAGGCCATCG‑GTCAGGCGCGAGGGCTCGATCAAAT‑‑CGCCATCAGACGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCTGTAGTAGCACGCCGGGGCGAGGTCCAGt   >  2:84861/1‑149 (MQ=255)                                                                                                                                                  caTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGTg  <  3:112953/149‑1 (MQ=255)                                                                                                                                                  caTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGTg  <  3:83185/149‑1 (MQ=255)                                                                                                                                                  caTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGTg  >  1:203371/1‑149 (MQ=255)                                                                                                                                                  caTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGTg  >  7:2715/1‑149 (MQ=255)                                                                                                                                                  caTTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGTg  >  8:69154/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTCCAGTG  >  NC_002947/5770298‑5770591

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