Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,757,455	A→T	5.7%	S93C (AGC→TGC) ‡	PP_3321 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,757,455	0	A	T	5.7%	215.0 / 2.5	70	S93C (AGC→TGC) ‡	PP_3321	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base T (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TGGCCTGCCTGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCGA  >  NC_002947/3757309‑3757597                                                                                                                                                    |                                                                                                                                                tGGCCTGCCTGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGc                                                                                                                                               >  8:128501/1‑149 (MQ=255)   gcctgcctGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCCCCTGATGCGA‑TCCTGGTTCCGACCCGGTCCCACCAAGGGCGTGACTGCGACGACAGCCg                                                                                                                                             >  8:164143/1‑149 (MQ=255)   gcctgcctGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGcc                                                                                                                                              <  7:256470/148‑1 (MQ=255)     ctgcctgGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTc                                                                                                                                           >  2:233858/1‑149 (MQ=255)      tgcctgGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCg                                                                                                                                          >  1:23584/1‑149 (MQ=255)         ctgGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCaa                                                                                                                                       <  4:176534/149‑1 (MQ=255)            gCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAAc                                                                                                                                      <  8:89242/147‑1 (MQ=255)            gCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGGGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTa                                                                                                                                    <  7:164143/149‑1 (MQ=255)                 tGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑g                                                                                                                              <  7:128501/149‑1 (MQ=255)                   ccAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTc                                                                                                                                >  4:107420/1‑146 (MQ=255)                   ccAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTc                                                                                                                                >  5:226450/1‑146 (MQ=255)                    cAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑gaga                                                                                                                           <  1:106181/149‑1 (MQ=255)                         cggcgCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTc                                                                                                                      <  5:155042/149‑1 (MQ=255)                                   ccGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTa                                                                                                               <  2:113687/146‑1 (MQ=255)                                   ccGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAac                                                                                                             >  3:250589/1‑148 (MQ=255)                                     gTTGGGCCGCGACCGTGGCCACACGTGATGAGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAacgac                                                                                                          >  1:29584/1‑149 (MQ=255)                                            cgcgACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCt                                                                                                                    <  8:93788/132‑1 (MQ=255)                                            cgcgACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCt                                                                                                                    >  7:93788/1‑132 (MQ=255)                                            cgcgACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTgg                                                                                                   >  6:270316/1‑149 (MQ=255)                                                aCCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGTCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCtg                                                                                               <  1:233858/149‑1 (MQ=255)                                                 ccGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCtgt                                                                                              <  1:296257/149‑1 (MQ=255)                                                   gTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGa                                                                                            <  4:250589/149‑1 (MQ=255)                                                    cggCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGAc                                                                                           <  3:63567/148‑1 (MQ=255)                                                       ccACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATc                                                                                        >  6:296704/1‑149 (MQ=255)                                                             gTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACAt                                                                                         <  8:94909/142‑1 (MQ=255)                                                             gTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACAt                                                                                         >  7:94909/1‑142 (MQ=255)                                                                   tgGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCag                                                                             <  1:267264/148‑2 (MQ=255)                                                                     gATCAGGTGCTGGCGTTGTACGCACCTGATGCGATTCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCaag                                                                            >  1:48056/1‑148 (MQ=255)                                                                           gTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCAGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGt                                                                    >  2:255858/1‑149 (MQ=255)                                                                              cTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAg                                                                 >  6:52326/1‑149 (MQ=255)                                                                                 gCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGcc                                                               <  5:296704/148‑1 (MQ=255)                                                                                        aCGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTg                                                       >  5:274402/1‑149 (MQ=255)                                                                                            aCCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGa                                                   >  6:90003/1‑149 (MQ=255)                                                                                             ccTGATGCGA‑TCCTGGTTCCGAACCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGAc                                                  <  6:226450/149‑1 (MQ=255)                                                                                               tGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCg                                                >  1:51441/1‑149 (MQ=255)                                                                                                 aTGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTg                                              <  3:107420/149‑1 (MQ=255)                                                                                                   gCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGt                                            >  7:40601/1‑149 (MQ=255)                                                                                                    cGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTc                                           <  2:29584/149‑1 (MQ=255)                                                                                                      a‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGt                                         >  7:29903/1‑149 (MQ=255)                                                                                                        tCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGtt                                        >  4:73648/1‑149 (MQ=255)                                                                                                          cTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTg                                       >  5:97147/1‑148 (MQ=255)                                                                                                          cTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTg                                       <  6:119585/148‑1 (MQ=255)                                                                                                             gTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAg                                    <  8:29903/148‑1 (MQ=255)                                                                                                              ttCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGc                                  <  5:52326/149‑1 (MQ=255)                                                                                                                 cGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCa                                <  2:51441/148‑1 (MQ=255)                                                                                                                 cGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCa                                <  3:16460/148‑1 (MQ=255)                                                                                                                     ccTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTaca                           >  2:216061/1‑149 (MQ=255)                                                                                                                      cTGTCCAACCAAGTCCGTGACTGAGACGACAGCCGTCGCAACTATATC‑GAGAATTTCCTCGCTAAAGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTaaac                          >  4:301724/1‑149 (MQ=255)                                                                                                                            aaccaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGGAGAATTTCCACGCAAACGACGGGCTGGTCTGTGACATCCAGTTGTTCaagaag                                                                         >  1:62846/1‑97 (MQ=255)                                                                                                                            aaccaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCaagaag                                                                         <  2:62846/96‑1 (MQ=255)                                                                                                                            aaccaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGa                                                   >  5:10065/1‑118 (MQ=255)                                                                                                                            aaccaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGa                                                   <  6:10065/118‑1 (MQ=255)                                                                                                                             accaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTc                   <  2:23584/149‑1 (MQ=255)                                                                                                                              ccaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTg                                       >  7:242873/1‑128 (MQ=255)                                                                                                                              ccaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTg                                       <  8:242873/128‑1 (MQ=255)                                                                                                                                     cGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTga           <  5:270316/149‑1 (MQ=255)                                                                                                                                        gACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGc                                  <  8:45499/123‑1 (MQ=255)                                                                                                                                        gACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGc                                  >  7:45499/1‑123 (MQ=255)                                                                                                                                          ttGCGACGGCTGTCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTg                                       >  1:125639/2‑116 (MQ=255)                                                                                                                                          ttGCGACGGCTGTCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTg                                       <  2:125639/115‑1 (MQ=255)                                                                                                                                          ttGCGACGGCTGTCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGacgacg      >  5:178872/2‑149 (MQ=255)                                                                                                                                           tGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACgg     >  2:5038/1‑149 (MQ=255)                                                                                                                                             cgacggCTGTCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGc    <  3:88203/148‑1 (MQ=255)                                                                                                                                             cgacgaCAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTTGGGACCGTGGTCGTTGGAGTCCAGTACACCTTCGTCTTTAGTGAAGACGGc    >  4:189841/1‑148 (MQ=255)                                                                                                                                             cgacgaCAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTc               <  2:233632/137‑1 (MQ=255)                                                                                                                                             cgacgaCAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCg   >  8:38569/1‑149 (MQ=255)                                                                                                                                             cgacgaCAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACAGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTc               >  1:233632/1‑137 (MQ=255)                                                                                                                                              gacgacAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCGa  <  2:48056/149‑1 (MQ=255)                                                                                                                                              gacgacAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCGa  <  5:90003/149‑1 (MQ=255)                                                                                                                                                   cAGCCGTCGCAACTATTTC‑TAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTc               <  4:306087/131‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TGGCCTGCCTGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGA‑TCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTC‑GAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCGA  >  NC_002947/3757309‑3757597

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