Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,646,266	A→T	7.7%	D540V (GAC→GTC)	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,646,266	0	A	T	7.7%	241.1 / 7.6	86	D540V (GAC→GTC)	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base T (5/2);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

CACCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAA  >  NC_002947/5646120‑5646407                                                                                                                                                   |                                                                                                                                                  cacCCTGGTCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACt                                                                                                                                                 <  7:155999/149‑1 (MQ=255) cacCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACt                                                                                                                                                 <  2:341577/149‑1 (MQ=255)  acCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTg                                                                                                                                                >  2:307595/1‑149 (MQ=255)   cccTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTgg                                                                                                                                               <  7:288786/149‑1 (MQ=255)   cccTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTgg                                                                                                                                               >  6:381609/1‑149 (MQ=255)        gCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAggg                                                                                                                                          >  4:362882/1‑149 (MQ=255)           aCTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGcc                                                                                                                                        >  3:243081/1‑148 (MQ=255)               ccGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACgg                                                                                                                                   <  1:238367/149‑1 (MQ=255)               ccGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACgg                                                                                                                                   >  3:309585/1‑149 (MQ=255)                 gcgATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCt                                                                                                                                 <  4:243081/149‑1 (MQ=255)                  cgATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTg                                                                                                                                <  7:100532/149‑1 (MQ=255)                         ccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACt                                                                                                                                                 <  1:228763/125‑1 (MQ=255)                         ccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACt                                                                                                                                                 >  2:228763/1‑125 (MQ=255)                          cTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGgc                                                                                                                        <  5:381609/149‑1 (MQ=255)                             ttGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCt                                                                                                                     <  4:137724/148‑1 (MQ=255)                             aTGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGgcg                                                                                                                       <  5:340295/147‑1 (MQ=255)                               ggTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTgg                                                                                                                   >  6:133508/1‑149 (MQ=255)                                gTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGCTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGc                                                                                                                  >  6:234384/1‑149 (MQ=255)                                   cTGTGCTGCTGCAACCAGATCGGCGATGGCTACGGCCTGGGGCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCaa                                                                                                               <  1:216610/149‑1 (MQ=255)                                           cTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATAccc                                                                                                          <  1:196850/146‑1 (MQ=255)                                           cTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTgg                                                                                                       >  3:165471/1‑149 (MQ=255)                                             tCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGc‑‑‑‑c                                                                                                 >  3:218629/1‑149 (MQ=255)                                             tCAACCAGATCGGCGATGGCTACGGCCTGTACCTGCCGCGCCTGTTAGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCATGCCTGACTGGAAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGc‑‑‑‑c                                                                                                 >  2:90283/1‑149 (MQ=255)                                               aaCCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTg                                                                                               >  4:282249/1‑149 (MQ=255)                                                 ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTgcg                                                                                             <  3:362882/149‑1 (MQ=255)                                                  cataTCTGCGATGGCTACGGCCTGTTCCTGCAGCGCCTGTTCTATACCCTGGTCGGCAGCGTTATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGCCGCTGGCCAATACCCTGGC‑‑‑‑CTgcgc                                                                                            <  2:226624/146‑1 (MQ=255)                                                    gATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCa                                                                                          <  4:218629/149‑1 (MQ=255)                                                      tCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGc                                                                                        <  2:312067/149‑1 (MQ=255)                                                          cGAGGGCGACGGCCTGGTCCTGCCGCGCGTGTTCGATACGCTGGTCGGCAGCCTGCTTGCCATTCTTGCAGTGTGCCTGTTCCTGCCTGACTGGCCGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCgtgt                                                                                    <  8:296085/149‑1 (MQ=255)                                                            aTGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTAc                                                                                  <  2:40523/149‑1 (MQ=255)                                                               gCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTg                                                                               <  6:334857/149‑1 (MQ=255)                                                                   cGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTgcgc                                                                            <  3:351492/148‑1 (MQ=255)                                                                   cGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTgcgc                                                                            <  2:102899/148‑1 (MQ=255)                                                                      ccTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTgcgc                                                                            >  6:149319/1‑145 (MQ=255)                                                                       cTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTgcgc                                                                            >  7:279540/1‑144 (MQ=255)                                                                           tCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                   >  5:180645/1‑149 (MQ=255)                                                                            ccTGCCGCGGCTGTTCTTTACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCGTGCCTGACTGGCAGGGCCGAGGGCTGAACAAGCCGCTGGCCAATAGCCTGGC‑‑‑‑GTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                   <  4:280366/148‑1 (MQ=255)                                                                            ccTGCCGCGCCTGCTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                   <  5:234384/148‑1 (MQ=255)                                                                              tGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATgcag                                                                >  7:284101/1‑149 (MQ=255)                                                                                ccGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGAc                                                                                                                                                  >  6:345034/1‑69 (MQ=255)                                                                                ccGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGAc                                                                                                                                                  <  5:345034/69‑1 (MQ=255)                                                                                ccGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATgcagca                                                              <  1:337699/149‑1 (MQ=255)                                                                                ccGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATgcagca                                                              <  8:112927/149‑1 (MQ=255)                                                                                ccGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATgcagca                                                              <  8:331681/149‑1 (MQ=255)                                                                                  gcgcCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCg                                                                                       <  1:164947/122‑1 (MQ=255)                                                                                  gcgcCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCg                                                                                       >  2:164947/1‑122 (MQ=255)                                                                                     ccTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATg                                                         <  4:142445/149‑1 (MQ=255)                                                                                      cTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGc                                                        <  8:284101/149‑1 (MQ=255)                                                                                       tGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGcc                                                       >  7:246465/1‑149 (MQ=255)                                                                                           cGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCAc                                                    <  6:180645/148‑1 (MQ=255)                                                                                           cGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACg                                                   >  6:135802/1‑149 (MQ=255)                                                                                            gATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCTTGCGTGACTGGCAGGGCCGACGGCTGAACAAGGGGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTTTGCCCACg                                                   <  6:174024/148‑1 (MQ=255)                                                                                                 ccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgc                                             >  8:351353/1‑149 (MQ=255)                                                                                                 ccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgc                                             >  6:122194/1‑149 (MQ=255)                                                                                                   tGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgcgc                                           >  7:54274/1‑149 (MQ=255)                                                                                                     gTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGcga                                         >  1:338924/1‑149 (MQ=255)                                                                                                             ccTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTgg                                 <  4:165471/149‑1 (MQ=255)                                                                                                              cTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGc                                >  7:75460/1‑149 (MQ=255)                                                                                                               tGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGcc                               >  5:197835/1‑149 (MQ=255)                                                                                                                  ttGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTAc                            >  2:114166/1‑149 (MQ=255)                                                                                                                     ccATTCTTGCATTGTTCCTGTTCCAGCCTGTCTGGCCGGGCCGACGGCGGTACAACTCGCTGGGCAAAACCCTGCCTCAGCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCt                              <  1:307595/148‑1 (MQ=255)                                                                                                                     ccATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGc                         >  4:113305/1‑149 (MQ=255)                                                                                                                     ccATTCTCGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCgacgac                                     <  7:329539/137‑1 (MQ=255)                                                                                                                     ccATTCTCGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCgacgac                                     >  8:329539/1‑137 (MQ=255)                                                                                                                      cATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGcc                        >  6:369055/1‑149 (MQ=255)                                                                                                                      cATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGcc                        <  8:279540/149‑1 (MQ=255)                                                                                                                       aTTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCt                       >  3:339352/1‑149 (MQ=255)                                                                                                                        ttCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTg                      >  2:47937/1‑149 (MQ=255)                                                                                                                             gCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCg                 >  7:152064/1‑149 (MQ=255)                                                                                                                              cAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCgg                >  5:269529/1‑149 (MQ=255)                                                                                                                               aGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgc               >  4:366535/1‑149 (MQ=255)                                                                                                                                gtgtTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgcg              >  6:287359/1‑149 (MQ=255)                                                                                                                                     gtcggccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATg         >  5:103970/7‑149 (MQ=255)                                                                                                                                     cctgttcctgCGTTACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCTTGCAGCAGTATGCCCACGTCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAAt          <  1:90283/148‑1 (MQ=255)                                                                                                                                       tgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTCCCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGcc       >  3:33109/1‑149 (MQ=255)                                                                                                                                       tgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGcc       >  4:371935/1‑149 (MQ=255)                                                                                                                                        gttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGcc       >  1:361779/1‑148 (MQ=255)                                                                                                                                         ttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGcc       <  5:149319/147‑1 (MQ=255)                                                                                                                                         gccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgc               <  8:138102/137‑1 (MQ=255)                                                                                                                                         gccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgc               >  7:138102/3‑139 (MQ=255)                                                                                                                                         gccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca     >  2:328355/3‑149 (MQ=255)                                                                                                                                          ccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca     >  8:351662/2‑148 (MQ=255)                                                                                                                                           cctgcctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTg                                  <  6:44247/118‑1 (MQ=255)                                                                                                                                           cctgcctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTg                                  >  5:44247/1‑118 (MQ=255)                                                                                                                                            ctgcctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACaa  <  2:116167/149‑1 (MQ=255)                                                                                                                                            ctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACaa  >  8:105576/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CACCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGC‑‑‑‑CTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAA  >  NC_002947/5646120‑5646407

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