Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,599,308	G→A	9.5%	G167G (GGC→GGT)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,599,308	0	G	A	9.5%	227.5 / 9.8	84	G167G (GGC→GGT)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base A (3/5);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

TACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAGG  >  NC_002947/4599164‑4599455                                                                                                                                                  |                                                                                                                                                    tACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGcccc                                                                                                                                                 >  2:358462/1‑149 (MQ=255)   ccGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGt                                                                                                                                               >  6:404695/1‑149 (MQ=255)    cGAAGACTCAACGTTAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                                              >  8:83633/1‑149 (MQ=255)    cGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                                              >  1:411301/1‑149 (MQ=255)     gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGCTAAAGGCACCGGGCGGGGCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgg                                                                                                                                             <  4:282044/149‑2 (MQ=37)     gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGc                                                                                                                                             <  4:391401/149‑1 (MQ=255)      aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGcc                                                                                                                                            >  6:374378/1‑149 (MQ=255)      aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTggg                                                                                                                                            <  4:77446/149‑3 (MQ=38)       aGACTCAACGCCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgggg                                                                                                                                           <  7:35316/149‑4 (MQ=37)         aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                         >  8:8978/1‑149 (MQ=255)         aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                         >  8:2671/1‑149 (MQ=255)         aCTCAACGTCAGGGGAGGCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                         <  3:151650/149‑1 (MQ=255)          cTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcg                                                                                                                                        <  2:35839/149‑1 (MQ=255)           tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                       >  5:201125/1‑149 (MQ=255)           tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                       >  7:20697/1‑149 (MQ=255)           tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                                 >  7:191986/1‑139 (MQ=38)           tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                                 <  8:191986/139‑1 (MQ=38)            cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                                      <  7:2671/149‑1 (MQ=255)            cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                                      <  7:8978/149‑1 (MQ=255)             aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                       <  6:378962/147‑1 (MQ=255)             aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                                     >  5:332642/1‑149 (MQ=255)             aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                                     >  5:52232/1‑149 (MQ=255)             aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                                     >  7:105538/1‑149 (MQ=255)                       gAGTCACTGCCCAGCGAATCTGTGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCg                                                                                                                           >  6:12035/1‑149 (MQ=255)                       gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGc                                                                                                                            <  6:81427/148‑1 (MQ=255)                          tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                                        >  8:51009/1‑149 (MQ=255)                          tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                                        <  5:54930/149‑1 (MQ=255)                              tGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                                    >  7:214131/1‑149 (MQ=255)                                cccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                                    >  7:377068/1‑147 (MQ=255)                                  cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                                                 >  7:132207/1‑148 (MQ=255)                                  cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACAAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                                                 >  7:339325/1‑148 (MQ=255)                                   aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                                >  5:122338/1‑148 (MQ=255)                                    gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                                              <  6:122338/149‑1 (MQ=255)                                        aTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACttt                                                                                                          <  2:411301/149‑1 (MQ=255)                                             cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTg                                                                                                         <  5:441286/145‑1 (MQ=255)                                             cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATc                                                                                                      <  6:201125/148‑1 (MQ=255)                                                ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  <  6:52232/149‑1 (MQ=255)                                                ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  <  5:12035/149‑1 (MQ=255)                                                ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  >  2:69170/1‑149 (MQ=255)                                                 gCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  >  8:31085/1‑148 (MQ=255)                                                  ccGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCa                                                                                                <  6:332642/149‑1 (MQ=255)                                                   cGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAg                                                                                               <  7:51009/149‑1 (MQ=255)                                                     ttCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTc                                                                                             <  7:83633/149‑1 (MQ=255)                                                         tgcagATAAACGCACCGGGGGGGTCGATGATCCCCGCATCATTGGTGATACACCCCTCCCTCGCGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCaa                                                                                           <  1:358462/144‑1 (MQ=255)                                                              aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                                       <  6:234841/130‑1 (MQ=255)                                                              aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                                       >  5:234841/1‑130 (MQ=255)                                                                 aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgctgct                                                                                 <  3:212510/149‑1 (MQ=255)                                                                  aGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGctc                                                                                <  1:231808/149‑1 (MQ=255)                                                                        cgggcgggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCg                                                                          <  7:31085/149‑1 (MQ=255)                                                                              ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                                    <  5:374378/149‑1 (MQ=255)                                                                              ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                                    >  4:425101/1‑149 (MQ=255)                                                                              ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                                    <  6:325981/149‑1 (MQ=255)                                                                                    tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                                              <  2:6542/149‑1 (MQ=255)                                                                                    tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                                              <  8:20697/149‑1 (MQ=255)                                                                                     gatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcgc                                                             >  6:313327/1‑149 (MQ=255)                                                                                       tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                                 >  1:401858/1‑143 (MQ=255)                                                                                       tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                                 <  2:401858/143‑1 (MQ=255)                                                                                        cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                                >  2:44751/1‑95 (MQ=255)                                                                                        cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                                <  1:44751/95‑1 (MQ=255)                                                                                         aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTATGCGTGCGCCTTCGCTGGGCGGGCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCa                                                          <  2:299253/148‑1 (MQ=255)                                                                                         aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAg                                                         >  4:247087/1‑149 (MQ=255)                                                                                               tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATGAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                                                   <  8:48266/149‑1 (MQ=255)                                                                                               tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                                                   >  7:257281/1‑149 (MQ=255)                                                                                                catTGCTGATCCACCAC‑CCCAGGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGGCAACGAGAACCTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGCgg                                                  >  8:386002/1‑149 (MQ=255)                                                                                                 atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAgg                                                  <  3:421283/148‑1 (MQ=255)                                                                                                 atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGt                                                 <  8:257281/149‑1 (MQ=255)                                                                                                     cTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCa                                             >  7:300901/1‑149 (MQ=255)                                                                                                             ccac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATc                                     <  5:270398/149‑1 (MQ=255)                                                                                                              cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                                    >  7:352522/1‑149 (MQ=255)                                                                                                              cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                                    >  4:430497/1‑149 (MQ=255)                                                                                                               ac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCgg                                   >  6:135943/1‑149 (MQ=255)                                                                                                                c‑cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATATAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGt                                  <  4:107898/149‑1 (MQ=255)                                                                                                                  cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa                                 <  8:132207/149‑1 (MQ=255)                                                                                                                   gcACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa                                 >  1:99909/2‑148 (MQ=38)                                                                                                                       ggggccggCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcg                            >  4:214933/7‑149 (MQ=38)                                                                                                                       ggggTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcg                            <  8:339325/149‑1 (MQ=255)                                                                                                                             ggCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGAc                      <  3:247087/149‑1 (MQ=255)                                                                                                                                     ccTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg              >  5:286600/1‑149 (MQ=255)                                                                                                                                     ccTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg              >  2:91759/1‑149 (MQ=255)                                                                                                                                     caggcgcTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg              <  2:99909/146‑1 (MQ=255)                                                                                                                                      cTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCgg             >  5:169797/1‑149 (MQ=255)                                                                                                                                          gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTTCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt         <  8:214131/149‑1 (MQ=255)                                                                                                                                          gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGGCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt         >  7:274543/1‑149 (MQ=255)                                                                                                                                          gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt         >  5:167623/1‑149 (MQ=255)                                                                                                                                                cggccCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAg   <  8:377068/149‑1 (MQ=255)                                                                                                                                                 ggccCCGTGCCTGAGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAgg  <  8:105538/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAGG  >  NC_002947/4599164‑4599455

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