Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,601,811:1	+C	6.2%	coding (1671/5418 nt)	PP_4924 ←	subtilase family serine protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,601,811	1	.	C	6.2%	331.9 / 13.5	84	coding (1671/5418 nt)	PP_4924	subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base . (38/40);  new base C (2/3);  total (41/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 = 6.60e-01

CCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGC  >  NC_002947/5601664‑5601957                                                                                                                                                        |                                                                                                                                                     ccAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCg                                                                                                                                                       >  7:334363/1‑148 (MQ=255)  cAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑g                                                                                                                                                  >  4:277764/1‑149 (MQ=255)   aaTCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑gc                                                                                                                                                 <  6:138567/149‑1 (MQ=255)     tcgtcgCTGGTACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑gcgc                                                                                                                                               >  5:347969/1‑149 (MQ=255)           cTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCGC‑A‑GCGCtcggc                                                                                                                                          <  8:334363/149‑5 (MQ=255)              gCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGGTGCTGCCACCGCCACCAACCAGC‑CACCTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAACGCGTCT‑‑GGACGCCGC‑A‑GCGCT‑CGGCGta                                                                                                                                      <  7:174063/149‑2 (MQ=255)                aCTGACCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCa                                                                                                                                    <  7:272651/149‑1 (MQ=255)                    aTCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGt                                                                                                                                >  8:44935/1‑149 (MQ=255)                     tCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGtg                                                                                                                               >  5:239407/1‑149 (MQ=255)                        cGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGtgctg                                                                                                                            >  3:367238/1‑149 (MQ=255)                         gAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGtgctgc                                                                                                                           >  7:402448/1‑149 (MQ=255)                          aaTACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGcc                                                                                                                          >  5:392126/1‑149 (MQ=255)                          aaTACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGcc                                                                                                                          <  7:354060/149‑1 (MQ=255)                              cATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCCCCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGa                                                                                                                      <  1:223049/149‑1 (MQ=255)                              cATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGa                                                                                                                      <  2:177052/149‑1 (MQ=255)                              cATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGa                                                                                                                      <  7:44935/149‑1 (MQ=255)                               aTCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCGC‑A‑GCGCT‑CGGCGTACAga                                                                                                                                  <  1:190802/136‑2 (MQ=255)                               aTCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCGC‑A‑GCGCT‑CGGCGTACAga                                                                                                                                  >  2:190802/1‑135 (MQ=255)                               aTCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTc                                                                                                                                      >  4:7176/1‑132 (MQ=255)                               aTCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTc                                                                                                                                      <  3:7176/132‑1 (MQ=255)                                    cgcccgccCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑GAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCGGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCActc                                                                                                                <  3:213970/149‑1 (MQ=255)                                    cgcccgccCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCActc                                                                                                                >  1:281881/1‑149 (MQ=255)                                     gcccgcccCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCActct                                                                                                               <  4:367238/149‑1 (MQ=255)                                        cgccccgccGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑TTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGc                                                                                                            <  1:95414/149‑1 (MQ=255)                                        cgccccgccGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGc                                                                                                            <  7:295165/149‑1 (MQ=255)                                          ccccgccGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCt                                                                                                           >  3:60935/1‑148 (MQ=255)                                           cccgccGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGa                                                                                                         >  7:280434/1‑149 (MQ=255)                                                cgTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAg                                                                                                    <  3:277764/149‑1 (MQ=255)                                                     tAGGTGTTGCTGCCACCGCCACCAACCAGCTCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGccc                                                                                                >  1:40842/1‑149 (MQ=255)                                                       ggTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTc                                                                                             <  2:188847/149‑1 (MQ=255)                                                        gtgtTGCTGCCACCGCCCCCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCg                                                                                            >  7:412038/1‑149 (MQ=255)                                                        gtgtTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCg                                                                                            >  4:201285/1‑149 (MQ=255)                                                          gtTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCggg                                                                                          <  3:38427/149‑1 (MQ=255)                                                          gtTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCACCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCgg                                                                                           >  8:81648/1‑148 (MQ=255)                                                              ctgcCACCGCCCCCACCCACC‑CAGCTGCCCTCTGGGGCCGCTGTCCAGAGATCCTCCGCCCCGCTGCCGTAAAGCGTCA‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGAGGTAGGCCCTTGGGGTAgg                                                                                      <  4:312986/149‑2 (MQ=255)                                                              ctgcCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAg                                                                                       <  5:85453/148‑1 (MQ=255)                                                              ctgcCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGt                                                                                      <  7:245300/149‑1 (MQ=255)                                                                gccaccgccaccAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTTGGGGGGTGG‑‑T‑CCGAGGGGGCGACCCCCGTGATGCCTTGAGCACTCTGGCTGATTTATGCCCTACGGGaa                                                                                        >  1:359597/1‑147 (MQ=255)                                                                 ccaccgccaccAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTg                                                                                     <  6:347969/147‑1 (MQ=255)                                                                  caccgccaccAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGc                                                                                  <  8:280434/149‑1 (MQ=255)                                                                   accgccaccAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGcc                                                                                 >  3:151382/1‑149 (MQ=255)                                                                    ccgccaccAACCAGC‑CAACTGCCCTCTGGGGCCGCGCTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGccc                                                                                >  5:129870/1‑149 (MQ=255)                                                                     cgccaccAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCg                                                                               >  4:111469/1‑149 (MQ=255)                                                                        caccaACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCtt                                                                            <  6:239407/149‑1 (MQ=255)                                                                         accaaccaGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAg                                                                                       <  6:208130/137‑1 (MQ=255)                                                                         accaaccaGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAg                                                                                       >  5:208130/1‑137 (MQ=255)                                                                         accaaccaGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTg                                                                           <  1:17536/149‑1 (MQ=255)                                                                          ccaaccaGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGa                                                                          <  8:402448/149‑1 (MQ=255)                                                                                aGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCGTAAGGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACg                                                                        >  1:353366/1‑148 (MQ=255)                                                                                     aaCTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCtt                                                                                              >  5:265197/1‑119 (MQ=255)                                                                                     aaCTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCtt                                                                                              <  6:265197/119‑1 (MQ=255)                                                                                     aaCTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAg                                                                 >  8:210847/1‑148 (MQ=255)                                                                                        tGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTa                                                               >  5:445254/1‑147 (MQ=255)                                                                                          cccTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGt                                                           <  3:304030/149‑1 (MQ=255)                                                                                            ctctGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGc                                                         >  6:158147/1‑149 (MQ=255)                                                                                             tctGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGt                                                                    <  4:119509/137‑1 (MQ=255)                                                                                             tctGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGt                                                                    >  3:119509/1‑137 (MQ=255)                                                                                                  ggCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGa                                                                                                         <  4:342156/95‑1 (MQ=255)                                                                                                  ggCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGa                                                                                                         >  3:342156/1‑95 (MQ=255)                                                                                                   gCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACcctc                                                  <  5:158147/149‑1 (MQ=255)                                                                                                     cgcgGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACcctcct                                                <  7:210847/149‑1 (MQ=255)                                                                                                         gTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGc                                            <  2:40842/149‑1 (MQ=255)                                                                                                            aaGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTc                                         >  8:130782/1‑149 (MQ=255)                                                                                                             aGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGTAGTCAGTAGAGTTCCACCCTCGTGGGCCTCg                                        >  8:195509/1‑149 (MQ=255)                                                                                                              gACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGa                                       <  4:60935/149‑1 (MQ=255)                                                                                                               aCATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTgg                                              >  1:209673/1‑141 (MQ=255)                                                                                                               aCATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTgg                                              <  2:209673/141‑1 (MQ=255)                                                                                                               aCATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGTCCTCGAt                                      >  4:383466/1‑149 (MQ=255)                                                                                                                catcatCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATc                                     <  6:392126/149‑1 (MQ=255)                                                                                                                catcatCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATc                                     <  6:129870/149‑1 (MQ=255)                                                                                                                 atcatcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATcc                                    <  4:151382/149‑1 (MQ=255)                                                                                                                  tcatcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCg                                   <  3:111469/149‑1 (MQ=255)                                                                                                                   catcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGaa                                                                                        >  2:65357/1‑95 (MQ=255)                                                                                                                   catcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGaa                                                                                        <  1:65357/95‑1 (MQ=255)                                                                                                                    atcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACcctc                                                  <  5:172611/132‑1 (MQ=255)                                                                                                                    atcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACcctc                                                  >  6:172611/1‑132 (MQ=255)                                                                                                                    atcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGAt                                 >  4:108879/1‑149 (MQ=255)                                                                                                                     tcGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGAt                                 >  5:263149/1‑148 (MQ=255)                                                                                                                           gcgcTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGAc                          >  5:320122/1‑149 (MQ=255)                                                                                                                            cgcTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACg                         >  6:163926/1‑149 (MQ=255)                                                                                                                            cgcTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACg                         <  6:445254/149‑1 (MQ=255)                                                                                                                                    tAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGc                 <  2:359597/149‑1 (MQ=255)                                                                                                                                        cacgcT‑‑GGACGCCGC‑A‑GCGCT‑CGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCgag              >  8:330379/5‑148 (MQ=255)                                                                                                                                                    gCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCg   >  3:159980/1‑149 (MQ=255)                                                                                                                                                      cG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGc  >  5:27923/1‑148 (MQ=255)                                                                                                                                                        |                                                                                                                                                     CCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGC‑CAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCT‑‑GTACGCCG‑‑A‑GCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGC  >  NC_002947/5601664‑5601957

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