Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,808,058	C→T	9.9%	V547V (GTG→GTA)	PP_3365 ←	acetolactate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,808,058	0	C	T	9.9%	163.1 / 6.1	60	V547V (GTG→GTA)	PP_3365	acetolactate synthase
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base T (3/3);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGAGCTGGGGTGAATGGGATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCCAGCAGCGGCAGCAGGCCTTCGGCGCTTTCCACGCGGTGG  >  NC_002947/3807914‑3808205                                                                                                                                                   |                                                                                                                                                    cTTGAGCTGGGGTGAATGGGATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTacac                                                                                                                                                   <  5:151363/147‑1 (MQ=255)    gAGCTGGGGTGAATGGGATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAACAGACGGCGCCGGCCTGTCAGGGCTTATACCGc                                                                                                                                                >  6:239659/1‑147 (MQ=255)    gAGCTGGGGTGAATGGGATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGc                                                                                                                                                >  6:82043/1‑147 (MQ=255)    gAGCTGGGGTGAATGGGATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGc                                                                                                                                                >  2:311350/1‑147 (MQ=255)               aaTGGGATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTc                                                                                                                                    <  1:37117/148‑1 (MQ=255)                  gggATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACg                                                                                                                                 >  5:93105/1‑148 (MQ=255)                    gATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCt                                                                                                                                     >  3:243175/1‑142 (MQ=255)                    gATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCt                                                                                                                                     <  4:243175/142‑1 (MQ=255)                       aCTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTATACCGCCAGGGCGCGCTCACGCAGCTc                                                                                                                           >  5:208659/1‑149 (MQ=255)                         tGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATACGGGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCCCCAGGGCGCGCTCACGCAGCTCg                                                                                                                          >  1:102895/1‑149 (MQ=255)                            agcagcCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTATACCGCCAGGGCGCGCTCACGCAGCTCGCTGt                                                                                                                      <  6:208659/149‑1 (MQ=255)                               agcCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACGGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTATACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGa                                                                                                                   <  5:239659/149‑1 (MQ=255)                               agcCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGa                                                                                                                   <  2:224470/149‑1 (MQ=255)                                gcCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAg                                                                                                                  <  2:495/149‑1 (MQ=255)                                 ccTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAg                                                                                                                  >  8:52178/1‑148 (MQ=255)                                 ccTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAg                                                                                                                  >  3:240199/1‑148 (MQ=255)                                        ttGCGAAAAGGCCCGTAAAGCCGACAGCATTACCGGCCGCTTTAC‑GGGCCGCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGCCGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCgg                                                                                                           >  1:4458/1‑149 (MQ=38)                                        ttGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCgg                                                                                                           <  6:93105/148‑1 (MQ=255)                                        ttGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGt                                                                                                          >  2:106583/1‑149 (MQ=255)                                          gCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGctg                                                                                                           <  8:279810/146‑3 (MQ=255)                                          gCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCgg                                                                                                           >  7:279810/1‑146 (MQ=255)                                           cGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCg                                                                                                        <  6:261725/148‑1 (MQ=255)                                            gAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCTTGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGtt                                                                                                      <  7:265843/149‑1 (MQ=255)                                            gAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGtt                                                                                                      >  8:310695/1‑149 (MQ=255)                                                aGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCg                                                                                                  <  1:311350/149‑1 (MQ=255)                                                  gCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCt                                                                                                <  5:82043/149‑1 (MQ=255)                                                    ccGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCt                                                                                                >  5:68823/1‑147 (MQ=255)                                                    ccGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCt                                                                                                <  6:68823/147‑1 (MQ=255)                                                    ccGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTg                                                                                               >  2:55634/1‑148 (MQ=255)                                                      gTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTaa                                                                                            <  4:241978/149‑1 (MQ=255)                                                      gTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTaa                                                                                            <  2:102895/149‑1 (MQ=255)                                                              gACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCg                                                                                                        <  5:108316/129‑1 (MQ=255)                                                              gACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCg                                                                                                        >  6:108316/1‑129 (MQ=255)                                                                   caTTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCTCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGa                                                                                        <  1:269918/141‑1 (MQ=255)                                                                   caTTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGa                                                                                        >  2:269918/1‑141 (MQ=255)                                                                    aTTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGt                                                                               >  5:174178/1‑149 (MQ=255)                                                                    aTTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGt                                                                               >  1:215294/1‑149 (MQ=255)                                                                        ccGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTTGCTGTAATCGACCGGGCAGTCGAt                                                                           <  7:52178/149‑1 (MQ=255)                                                                        ccGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGAt                                                                           <  7:310695/149‑1 (MQ=255)                                                                           tCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTATACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACTGCCCGGTCGAtta                                                                         >  3:119207/1‑146 (MQ=37)                                                                           tCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCAc                                                                        <  7:78560/149‑1 (MQ=255)                                                                             aGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTATACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACTGCCCGGTCGATTAca                                                                       <  4:119207/148‑2 (MQ=25)                                                                             aGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACTGCCCGGTCGATTAca                                                                       <  1:106583/148‑2 (MQ=25)                                                                              gCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTg                                                                     <  6:101201/149‑1 (MQ=255)                                                                              gCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTg                                                                     <  2:332985/149‑1 (MQ=255)                                                                                      gggCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAgg                                                                                                                                            <  7:210649/71‑1 (MQ=255)                                                                                      gggCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAgg                                                                                                                                            >  8:210649/1‑71 (MQ=255)                                                                                       ggCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTg                                                             <  4:16189/149‑1 (MQ=255)                                                                                          ctctTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTgggg                                                          >  5:2098/1‑149 (MQ=255)                                                                                           tctTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGt                                                         >  3:54971/1‑149 (MQ=255)                                                                                                 cagcaCCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACAGCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGAt                                                   <  4:240199/149‑1 (MQ=255)                                                                                                 cagcaCCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGAt                                                   >  1:48509/1‑149 (MQ=255)                                                                                                        aGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGc                                            <  2:4458/149‑1 (MQ=255)                                                                                                                ccGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCcagcag                                    <  1:55634/149‑1 (MQ=255)                                                                                                                         aGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCCAGCAGCGGcagcag                           <  2:48509/149‑1 (MQ=255)                                                                                                                           aCGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCGCGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCCAGCAGCGGCAGCAGGc                         >  6:54239/1‑149 (MQ=255)                                                                                                                                 ccGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCCAGCAGCGGCAGCAGGCCTTCgg                   <  6:2098/149‑1 (MQ=255)                                                                                                                                  cGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCg                                                                             <  5:104630/90‑1 (MQ=255)                                                                                                                                  cGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCg                                                                             >  6:104630/1‑90 (MQ=255)                                                                                                                                             ggCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTcc                                         >  4:121589/1‑115 (MQ=34)                                                                                                                                             ggCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTcc                                         <  3:121589/115‑1 (MQ=34)                                                                                                                                                  acacCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCCAGCAGCGGCAGCAGGCCTTCGGCGCTTTCCACGCGGTgg  >  1:125579/1‑149 (MQ=32)                                                                                                                                                  acacCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCCAGCAGCGGCAGCAGGCCTTCGGCGCTTTCCACGCGGTg   <  6:232580/148‑1 (MQ=32)                                                                                                                                                   |                                                                                                                                                    CTTGAGCTGGGGTGAATGGGATACTGAAGCAGCCTTGTATTGCGAAAAGGCCCGTAAAGCCGACA‑CATTACCGTCAGCTATAC‑GGGCCTCTTTGCAGCACCAGGCCGCGCCGGCAAAAAGACGGCGCCGGCCTGTCAGGGCTTACACCGCCAGGGCGCGCTCACGCAGCTCGCTGTTGAGGATGCGGTCGTTCTCGCTGTAATCGACCGGGCAGTCGATCACGTGCACACCTGGGGTCTTGATGCAGTGCTCCAGCAGCGGCAGCAGGCCTTCGGCGCTTTCCACGCGGTGG  >  NC_002947/3807914‑3808205

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