Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	21,472	C→T	7.0%	D425D (GAC→GAT)	PP_0016 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	21,472	0	C	T	7.0%	204.0 / 4.2	71	D425D (GAC→GAT)	PP_0016	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

TAAGCCTCTGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCGGTAGGGGCGTG  >  NC_002947/21324‑21615                                                                                                                                                      |                                                                                                                                                 tAAGCCTCTGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGAc                                                                                                                                                  >  5:53725/1‑149 (MQ=255) tAAGCCTCTGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGAc                                                                                                                                                  >  3:143459/1‑149 (MQ=255) tAAGCCTCTGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGAc                                                                                                                                                  >  3:137516/1‑149 (MQ=255) tAAGCCTCTGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATAGATTGATGAGAAGCTGCCTACCA‑GCTCTGAc                                                                                                                                                  >  4:356119/1‑149 (MQ=255)   aGCCTCTGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGAc                                                                                                                                                  <  8:591/147‑1 (MQ=255)         tGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCt                                                                                                                                          >  6:357342/1‑149 (MQ=255)           tCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTgg                                                                                                                                        >  3:143584/1‑149 (MQ=255)           tCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTg                                                                                                                                         <  5:357342/148‑1 (MQ=255)           tCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTg                                                                                                                                         <  3:346307/148‑1 (MQ=255)             gCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGgtg                                                                                                                                      >  1:31033/1‑149 (MQ=255)             gCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGgt                                                                                                                                       <  1:181746/148‑1 (MQ=255)                    ggTTAGGCAGCTGCTAGCTGGTGTTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCa                                                                                                                               <  4:143459/149‑1 (MQ=255)                    ggTTAGGCAGCTGCTAGCTGGGGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCa                                                                                                                               <  4:137516/149‑1 (MQ=255)                     gTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCaa                                                                                                                              <  3:401152/149‑1 (MQ=255)                      ttAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAac                                                                                                                             <  6:396868/149‑1 (MQ=255)                          gCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAactact                                                                                                                         >  2:55008/1‑149 (MQ=255)                              ctgctAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGc                                                                                                                     <  2:31033/149‑1 (MQ=255)                                 ctAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCAGGCTCTGGTCAGAGctg                                                                                                                                          <  2:185539/126‑3 (MQ=255)                                 ctAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATAGATGAGAAGCTGCCTACCAGGCTCTGGTCAGAGctg                                                                                                                                          >  1:185539/1‑124 (MQ=255)                                    gCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCa                                                                                                               >  8:343588/1‑149 (MQ=255)                                      tggtggTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCAtt                                                                                                             >  6:297634/1‑149 (MQ=255)                                              cccTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCCAGCAATCTCTG‑cc                                                                                                     >  8:347962/1‑148 (MQ=255)                                                 tCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCaa                                                                                                                  <  8:180022/133‑1 (MQ=255)                                                 tCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCaa                                                                                                                  >  7:180022/1‑133 (MQ=255)                                                   cAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAgg                                                                                                <  8:265224/148‑1 (MQ=255)                                                         tttGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAac                                                                                                                             >  8:256083/1‑114 (MQ=255)                                                         tttGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAac                                                                                                                             <  7:256083/114‑1 (MQ=255)                                                            ggtggAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGgggtgggt                                                                                       <  6:53725/148‑1 (MQ=255)                                                            ggtggAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGtt                                                                                      >  2:377563/1‑149 (MQ=255)                                                             gtggAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGATCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTg                                                                                     >  5:222147/1‑149 (MQ=255)                                                             gtggAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTg                                                                                     <  4:143584/149‑1 (MQ=255)                                                             gtggAGAACCGCGCCCAGGAGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTg                                                                                     <  3:12637/149‑1 (MQ=255)                                                                gagaACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGATCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTc                                                                                  <  6:222147/149‑1 (MQ=255)                                                                  gaACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGa                                                                                <  5:150092/149‑1 (MQ=255)                                                                    aCCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGt                                                                              <  7:343588/149‑1 (MQ=255)                                                                     ccGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGa                                                                                >  2:312644/1‑146 (MQ=255)                                                                     ccGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGa                                                                                <  1:312644/146‑1 (MQ=255)                                                                       gcgcCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTgg                                                                           <  3:356119/149‑1 (MQ=255)                                                                        cgcCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGt                                                                          <  1:377563/149‑1 (MQ=255)                                                                             aGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAgg                                                                     >  3:216320/1‑149 (MQ=255)                                                                              ggCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTGCCCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAgg                                                                     <  1:55008/149‑1 (MQ=255)                                                                                       tGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAgggggg                                                                                            >  1:397110/1‑116 (MQ=255)                                                                                       tGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAgggggg                                                                                            <  2:397110/116‑1 (MQ=255)                                                                                       tGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGc                                                           <  5:297634/149‑1 (MQ=255)                                                                                           tGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGc                                                           <  4:253281/145‑1 (MQ=255)                                                                                           tGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGc                                                           >  3:253281/1‑145 (MQ=255)                                                                                           tGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGt                                                       >  8:332791/1‑149 (MQ=255)                                                                                            gATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAgggg                                                                                              <  8:251699/109‑1 (MQ=255)                                                                                                   tGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTgg                                                                           <  5:162567/121‑1 (MQ=255)                                                                                                   tGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTgg                                                                           >  6:162567/1‑121 (MQ=255)                                                                                                            aTCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTa                                         >  6:45795/1‑146 (MQ=255)                                                                                                            aTCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTa                                         <  5:45795/146‑1 (MQ=255)                                                                                                            aTCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCT‑GTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGc                                     >  2:256847/1‑149 (MQ=255)                                                                                                             tCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGc                                     >  2:11326/1‑149 (MQ=255)                                                                                                             tCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTAGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGc                                     >  4:364984/1‑149 (MQ=255)                                                                                                               aaCGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGcc                                    >  3:78410/1‑148 (MQ=255)                                                                                                                        attgatGAGAAGCGGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCACCTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGaa                          <  6:348473/149‑1 (MQ=255)                                                                                                                            atgaGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGc                      >  5:364243/1‑149 (MQ=255)                                                                                                                            atgaGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTAGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATgg                       <  3:364984/148‑1 (MQ=255)                                                                                                                             tgaGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGGTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCt                     >  3:224943/1‑149 (MQ=255)                                                                                                                              gagaAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTc                    >  4:36786/1‑149 (MQ=255)                                                                                                                               agaATCTGCCTACCA‑GCTCTGACCTGAGCCTGGTGTCGATCAACTACTTGGCGAAGCATTCTCTG‑CCGAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATAGGTGCAGGTGAACTTCGCGCTTACGGCCCCGCTGAGATTGGCTc                    <  7:347962/148‑1 (MQ=255)                                                                                                                                gaAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCgg                  <  4:216320/149‑1 (MQ=255)                                                                                                                                    cTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGaa                          >  3:32134/1‑137 (MQ=255)                                                                                                                                    cTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGaa                          <  4:32134/137‑1 (MQ=255)                                                                                                                                    cTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGc              >  2:189572/1‑149 (MQ=255)                                                                                                                                     tGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAAAATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCg             >  1:100350/1‑149 (MQ=255)                                                                                                                                       ccTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACTTGCAGGTGAACTTCGCGATTACGGCCCAGCTGAGAAATGCTCGGTTGCgg            >  2:17183/1‑148 (MQ=255)                                                                                                                                       ccTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCg             <  2:100350/147‑1 (MQ=255)                                                                                                                                          aCCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGGCAAGCATTCTCTG‑CCCAAGGGGGGAGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGGCCAGCTGAGAATGGCTCGGTTGCGGTAgg        >  5:338264/1‑149 (MQ=255)                                                                                                                                           ccA‑GCTCTGACCAGAGCCT‑GTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCGGTAggg       <  1:256847/148‑1 (MQ=255)                                                                                                                                            caggCTCTGGTCAGAG‑CTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCGGTAgggg      >  2:50932/4‑149 (MQ=255)                                                                                                                                            cA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCGGTAgggg      >  6:214790/1‑149 (MQ=255)                                                                                                                                             a‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCGGTAGGGGc     >  7:12525/1‑149 (MQ=255)                                                                                                                                                ttctGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGATGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGATGCGGTAGGGGCGt   >  2:183312/2‑149 (MQ=255)                                                                                                                                                   tGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCGGTAGGGGCGTg  <  6:364243/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 TAAGCCTCTGTCGCAATCGGGTTAGGCAGCTGCTAGCTGGTGGTACCCTCCAGTTCTTTGGTGGAGAACCGCGCCCAGGCGAAAGGTGGCTGATTGACTGCCGCTCGATCAACGAATTGATTGATGAGAAGCTGCCTACCA‑GCTCTGACCAGAGCCTGGTGTCGATCAACTACTTGGCCAAGCATTCTCTG‑CCCAAGGGGGGTGGGTTGGTCGAGTTGGTTCAGGCAATACGTGCAGGTGAACTTCGCGCTTACGGCCCAGCTGAGAATGGCTCGGTTGCGGTAGGGGCGTG  >  NC_002947/21324‑21615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑