Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,066,335	T→C	7.2%	S89G (AGC→GGC)	PP_0924 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,066,335	0	T	C	7.2%	193.2 / 4.8	69	S89G (AGC→GGC)	PP_0924	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GGCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCG  >  NC_002947/1066189‑1066483                                                                                                                                                   |                                                                                                                                                     ggCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGgccgc                                                                                                                                                    >  7:322234/1‑149 (MQ=255)      gTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTggcg                                                                                                                                               <  5:103048/149‑1 (MQ=255)            cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCATCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTgg                                                                                                                                         >  6:128201/1‑149 (MQ=255)              cAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgggc                                                                                                                                       <  8:322234/149‑3 (MQ=255)                aaGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGAt                                                                                                                                     <  1:219914/149‑1 (MQ=255)                aaGCACGGCGATGGAACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGAt                                                                                                                                     <  8:321813/149‑1 (MQ=255)                 aGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                    <  2:7947/149‑1 (MQ=255)                 aGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                    <  3:392610/149‑1 (MQ=255)                  gCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                    >  7:49730/1‑148 (MQ=255)                   cACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGa                                                                                                                                  <  8:294929/149‑1 (MQ=255)                        cGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAgg                                                                                                                             <  3:138665/149‑1 (MQ=255)                               aCCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtgg                                                                                                                      >  7:1815/1‑149 (MQ=255)                                 cgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                    >  1:158458/1‑149 (MQ=255)                                 cgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                    <  6:177254/149‑1 (MQ=255)                                 cgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                    >  8:184409/1‑149 (MQ=255)                                              gTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGc                                                                                                       <  8:49730/149‑1 (MQ=255)                                               tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGctgctg                                                                                                                                                  >  4:386297/1‑105 (MQ=255)                                               tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGctgctg                                                                                                                                                  <  3:386297/105‑1 (MQ=255)                                               tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgg                                                                                                                                         <  2:355962/114‑1 (MQ=255)                                               tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgg                                                                                                                                         >  1:355962/1‑114 (MQ=255)                                                 gggggCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGAAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCg                                                                                                     <  7:184409/148‑1 (MQ=255)                                                          gCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGt                                                                                           >  2:8555/1‑149 (MQ=255)                                                           cAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTc                                                                                          <  2:49584/149‑1 (MQ=255)                                                           cAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTc                                                                                          <  4:9717/149‑1 (MQ=255)                                                            aGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTcg                                                                                         >  2:118197/1‑149 (MQ=255)                                                               tACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTACAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTcgcg                                                                                       >  2:67572/1‑148 (MQ=255)                                                                 cGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTa                                                                                     <  3:282230/148‑1 (MQ=255)                                                                 cGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTAGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTa                                                                                     <  1:310655/148‑1 (MQ=255)                                                                   gCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCATGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTc                                                                                                         >  5:110111/1‑126 (MQ=255)                                                                   gCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCATGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTc                                                                                                         <  6:110111/126‑1 (MQ=255)                                                                    cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                  >  5:152013/1‑148 (MQ=255)                                                                    cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                  >  7:68427/1‑148 (MQ=255)                                                                    cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                  <  8:68427/148‑1 (MQ=255)                                                                            gAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggc                                                                         >  2:8186/1‑149 (MQ=255)                                                                            gAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggc                                                                         >  2:298524/1‑149 (MQ=255)                                                                                   gggCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGt                                                                   <  3:365536/148‑1 (MQ=255)                                                                                    ggCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTg                                                                  >  2:113083/1‑148 (MQ=255)                                                                                      ccTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAc                                                               >  3:165354/1‑149 (MQ=255)                                                                                       cTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAcg                                                              <  3:196632/149‑1 (MQ=255)                                                                                       cTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAcg                                                              <  1:45022/149‑1 (MQ=255)                                                                                        ttGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGc                                                                 >  3:281084/1‑145 (MQ=255)                                                                                        ttGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGc                                                                 <  4:281084/145‑1 (MQ=255)                                                                                            aGCTACAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTa                                                         <  1:67572/149‑1 (MQ=255)                                                                                             gCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAg                                                        >  3:355911/1‑149 (MQ=255)                                                                                                gCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAgtcg                                                     >  4:45704/1‑149 (MQ=255)                                                                                                 cAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAgtcg                                                     <  1:298524/148‑1 (MQ=255)                                                                                                  aaCAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAgtcgtc                                                   >  4:292168/1‑149 (MQ=255)                                                                                                                gCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCt                                     >  1:194144/1‑149 (MQ=255)                                                                                                                 cTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTg                                    <  8:1815/149‑1 (MQ=255)                                                                                                                    aaTTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTc                                 <  3:341549/149‑1 (MQ=255)                                                                                                                      ttcttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCgg                               >  3:313837/1‑149 (MQ=255)                                                                                                                      ttcttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCgg                               >  7:187229/1‑149 (MQ=255)                                                                                                                         ttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAc                            >  6:52245/1‑149 (MQ=255)                                                                                                                          tcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAcc                           >  4:71535/1‑149 (MQ=255)                                                                                                                               gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTg                       <  3:45704/148‑1 (MQ=255)                                                                                                                               gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTg                       <  7:245982/148‑1 (MQ=255)                                                                                                                                cTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTTGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAg                     <  5:128201/149‑1 (MQ=255)                                                                                                                                   tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTGGTAGTCGTCCATGTCCAGCCCCTGCTCGGCACCtttta                      >  6:265007/1‑143 (MQ=255)                                                                                                                                   tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCtt                  >  7:37856/1‑149 (MQ=255)                                                                                                                                   tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTGTGAGCtt                  >  4:231607/1‑149 (MQ=255)                                                                                                                                   tAGAAACCGCCCAGCCTGCTGGCGGCCTGGTGATCGAGCAGGCGGCTGGTGGCAACGAGTTCGCCGGTGCTCAGGACGCGTACACCAATGTGCCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGa                      <  5:265007/145‑1 (MQ=255)                                                                                                                                     gAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTg                 <  3:400981/148‑1 (MQ=255)                                                                                                                                      tcgatcaccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGc               >  2:203550/8‑149 (MQ=255)                                                                                                                                      aaaCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGc               >  8:212443/1‑149 (MQ=255)                                                                                                                                        aCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCt             <  2:194144/149‑1 (MQ=255)                                                                                                                                         cGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt            <  1:8186/149‑1 (MQ=255)                                                                                                                                         atcaccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt            >  7:55355/5‑149 (MQ=255)                                                                                                                                           gCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGa          <  4:165354/149‑1 (MQ=255)                                                                                                                                            cccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGaa         <  6:152013/149‑1 (MQ=255)                                                                                                                                             ccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAt        >  6:20476/1‑149 (MQ=255)                                                                                                                                             ccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGCTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAt        >  1:266481/1‑149 (MQ=255)                                                                                                                                                ggCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCt     >  5:1112/1‑149 (MQ=255)                                                                                                                                                   tgctgGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTgcg  >  3:328541/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCG  >  NC_002947/1066189‑1066483

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