Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,044,485	G→A	6.6%	D519N (GAC→AAC)	PP_4448 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,044,485	0	G	A	6.6%	272.8 / 5.5	91	D519N (GAC→AAC)	PP_4448	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (45/40);  new base A (3/3);  total (48/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 = 2.48e-01

GCGTAAATACGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGACATTCCATACTGG  >  NC_002947/5044343‑5044628                                                                                                                                               |                                                                                                                                                 gagTAAATACGTTGTGACCTTGTTGCTGTCCAGTCAACGACCATCCGAGATATCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGt                                                                                                                                            <  8:321609/147‑1 (MQ=255) gCGTAAATACGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACGATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTTCGCAGGCTGCTGCCCGACACGt                                                                                                                                            >  7:321609/1‑149 (MQ=255)  cGTAAATACGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGtt                                                                                                                                           >  8:68117/1‑149 (MQ=255)      aaTACGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTggg                                                                                                                                        >  3:169777/1‑148 (MQ=255)      aaTACGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTggg                                                                                                                                        >  5:46727/1‑148 (MQ=255)          cGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAg                                                                                                                                    <  8:292616/148‑1 (MQ=255)          cGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAg                                                                                                                                    <  7:68117/148‑1 (MQ=255)          cGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAg                                                                                                                                    <  2:134576/148‑1 (MQ=255)          cGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAg                                                                                                                                    <  7:302474/148‑1 (MQ=255)                 aCCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTga                                                                                                                            <  4:169777/149‑1 (MQ=255)                   cTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTgaga                                                                                                                          <  1:71669/149‑1 (MQ=255)                     tgttgCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGaa                                                                                                                         >  7:148380/1‑148 (MQ=255)                       ttgCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAAt                                                                                                                      >  8:6582/1‑149 (MQ=255)                           tGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGtt                                                                                                                                           >  1:48599/1‑124 (MQ=255)                           tGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGtt                                                                                                                                           <  2:48599/124‑1 (MQ=255)                           tGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGcc                                                                                                                   >  5:136787/1‑148 (MQ=255)                             tCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGaaaa                                                                                                                       >  2:267976/1‑142 (MQ=255)                             tCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGaaaa                                                                                                                       <  1:267976/142‑1 (MQ=255)                             tCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCa                                                                                                                 >  8:227159/1‑148 (MQ=255)                             tCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCAACGTGTCGGGCagcag                                                                                                                                  >  2:226936/1‑127 (MQ=255)                             tCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCAACGTGTCGGGCagcag                                                                                                                                  <  1:226936/131‑5 (MQ=255)                                  gcAACCCCCCTCCCCGATCTCAGAGACCATCTTCACTCAACGCAACAACTATGCCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTCCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑gg                                                                                                          <  2:291950/148‑1 (MQ=255)                                    aaCGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGgc                                                                                                        >  4:206639/1‑149 (MQ=255)                                     aCGCCCATCCGAGATCTCAGAGACCATCTTCTGTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGAt                                                                                                                               <  5:70674/126‑1 (MQ=255)                                     aCGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTAAGTGCGCAGGCTGCTGACCGACACGTTGGGCAAGCTGag                                                                                                                               >  6:70674/1‑125 (MQ=255)                                       gCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGcc                                                                                                     >  1:176068/1‑149 (MQ=255)                                        cccATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCg                                                                                                    <  8:93583/149‑1 (MQ=255)                                           aTCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGt                                                                                                  >  3:298940/1‑148 (MQ=255)                                            tCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCAACGTGTCGGGCag                                                                                                                                     <  3:313998/113‑2 (MQ=255)                                            tCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCAACGTGTCGGGCag                                                                                                                                     >  4:313998/1‑112 (MQ=255)                                            tCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTAAGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTg                                                                                                 >  5:211906/1‑148 (MQ=255)                                                   tctcAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGGCCGACACCTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGctc                                                                                         >  2:105002/1‑149 (MQ=255)                                                     tcAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCtt                                                                                       >  5:65757/1‑149 (MQ=255)                                                        gagaCCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGAt                                                                                    <  6:46727/149‑1 (MQ=255)                                                           aCCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCgg                                                                                 <  2:320470/149‑1 (MQ=255)                                                           aCCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCgg                                                                                 >  1:320470/1‑149 (MQ=255)                                                            ccATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCggg                                                                                <  7:286444/149‑1 (MQ=255)                                                            ccATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCggg                                                                                <  6:136787/149‑1 (MQ=255)                                                            ccATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCggg                                                                                <  3:206639/149‑1 (MQ=255)                                                             cATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTgg                                                                                                                                         <  1:142808/92‑1 (MQ=255)                                                             cATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTgg                                                                                                                                         >  2:142808/1‑92 (MQ=255)                                                             cATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTAAGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGAATGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCggg                                                                                >  2:96234/1‑148 (MQ=255)                                                              aTCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGa                                                                              >  1:192336/1‑149 (MQ=255)                                                                cttctCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGAtt                                                                            >  2:269417/1‑149 (MQ=255)                                                                   ctctcAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGAtt                                                                            <  2:176068/146‑1 (MQ=255)                                                                   ctctcAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGAtt                                                                            <  1:269417/146‑1 (MQ=255)                                                                      tcAATGCAACAACTTTGTCGCAATGCGCCTTACGACCCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGCCACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTgg                                                                      >  7:54905/1‑149 (MQ=255)                                                                      tcAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCAc                                                                         <  8:148380/146‑1 (MQ=255)                                                                      tcAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTgg                                                                      >  6:250088/1‑149 (MQ=255)                                                                        aaTGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTgag                                                                                                                           >  4:180381/1‑95 (MQ=255)                                                                        aaTGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTgag                                                                                                                           >  6:339392/1‑95 (MQ=255)                                                                        aaTGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTgag                                                                                                                           <  3:180381/95‑1 (MQ=255)                                                                        aaTGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTgag                                                                                                                           <  5:339392/95‑1 (MQ=255)                                                                        aaTGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTg                                                                       <  7:6582/146‑1 (MQ=255)                                                                          tGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGAt                                                                  <  7:227159/149‑1 (MQ=255)                                                                          tGCAACAACTTTGTCACAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGATAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGAt                                                                  <  6:290563/149‑1 (MQ=255)                                                                             aacaacTTTGTCGCAAAGCCCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGATGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCCGGGATTCACTGGTGATGc                                                                >  6:322368/1‑148 (MQ=255)                                                                                     tGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTg                                                        <  2:192336/148‑1 (MQ=255)                                                                                     tGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTg                                                        <  6:65757/148‑1 (MQ=255)                                                                                        cGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTg                                                     >  2:43981/1‑148 (MQ=255)                                                                                         gCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGaa                                                   <  5:250088/149‑1 (MQ=255)                                                                                          cAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTTGGGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGaa                                                   <  1:96234/149‑1 (MQ=255)                                                                                          cAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCAACTCTGGTg                                                     >  7:116269/1‑146 (MQ=255)                                                                                           aaTGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGctc                                                                                         <  7:154483/109‑1 (MQ=255)                                                                                           aaTGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGctc                                                                                         >  8:154483/1‑109 (MQ=255)                                                                                           aaTGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTAGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTg                                                     >  2:8400/1‑145 (MQ=255)                                                                                                     aCGAACCCATGCTTTCAGTATTATGTGCGCAGGCTGCTGTCGGACACGTTGGGCAAGCTGATTGAGAATATGCCCAACTT‑GGGCGCCTGTGAAGCGCTGTTGTTCGGGGATTCACTGGTTATGCCATCTCTGGTGAAGATCGACGTc                                         <  5:322368/147‑1 (MQ=255)                                                                                                     aCGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCg                                              >  3:154792/1‑142 (MQ=255)                                                                                                     aCGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCg                                              <  4:154792/142‑1 (MQ=255)                                                                                                       gAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTg                                                                       >  3:122321/1‑115 (MQ=255)                                                                                                       gAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTg                                                                       <  4:122321/115‑1 (MQ=255)                                                                                                          cccAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGAAAGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGAt                                   >  4:29897/1‑148 (MQ=255)                                                                                                          cccAAGCGATCAGAATTATGTGCGCAGGCTCCTGCCCGACACCTTCCGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCCGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATc                                  >  4:324903/1‑149 (MQ=255)                                                                                                             aaGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAg                               >  1:191661/1‑149 (MQ=255)                                                                                                                           aTGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGACTGc                  <  6:91413/148‑1 (MQ=255)                                                                                                                           aTGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGACTGc                  >  5:91413/1‑148 (MQ=255)                                                                                                                              tgCTCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGAc              <  8:116269/149‑1 (MQ=255)                                                                                                                              tgCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGa               <  4:298940/148‑1 (MQ=255)                                                                                                                               gcgcAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGACa             >  4:70399/1‑149 (MQ=255)                                                                                                                                 gcAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATc                        >  3:292202/1‑136 (MQ=255)                                                                                                                                 gcAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATc                        <  4:292202/136‑1 (MQ=255)                                                                                                                                   aGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGACATTc          >  6:337939/1‑148 (MQ=255)                                                                                                                                     gctgctGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACt                                                                        <  4:80955/84‑1 (MQ=255)                                                                                                                                     gctgctGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACt                                                                        >  3:80955/1‑84 (MQ=255)                                                                                                                                      ctgctgCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GTGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCa                                <  5:118590/123‑1 (MQ=255)                                                                                                                                      ctgctgCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCa                                >  6:118590/1‑123 (MQ=255)                                                                                                                                        gctgcCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGACATTCCATACt    >  6:358713/1‑149 (MQ=255)                                                                                                                                         ttgcCCAACGTGTCGGGCAAGCTGATTGAGAAAATGCCCAACTT‑gg                                                                                                          <  4:323431/45‑1 (MQ=21)                                                                                                                                         ttgcCCAACGTGTCGGGCAAGCTGATTGAGAAAATGCCCAACTT‑gg                                                                                                          >  3:323431/2‑46 (MQ=21)                                                                                                                                           gcCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGACATTCCATACTgg  >  3:106201/1‑148 (MQ=255)                                                                                                                                            cccGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGcc                 >  8:95692/1‑132 (MQ=255)                                                                                                                                               |                                                                                                                                                 GCGTAAATACGGTGTGACCTTGTTGCTGTCCAGTCAACGCCCATCCGAGATCTCAGAGACCATCTTCTCTCAATGCAACAACTTTGTCGCAATGCGCCTTACGAACCCAAGCGATCAGAATTATGTGCGCAGGCTGCTGCCCGACACGTTGGGCAAGCTGATTGAGAAAATGCCCAACTT‑GGGCGCCGGTGAAGCGCTCTTGATCGGGGATTCACTGGTGATGCCATCTCTGGTGAAGATCGACGTCTGCGATCCAGCACCATCGTCTGCCGACATTCCATACTGG  >  NC_002947/5044343‑5044628

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