Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,143	T→C	6.0%	W37R (TGG→CGG) ‡	PP_3533 →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,004,143	0	T	C	6.0%	249.0 / 3.7	83	W37R (TGG→CGG) ‡	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base T (40/38);  new base C (3/2);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGAT  >  NC_002947/4004002‑4004291                                                                                                                                              |                                                                                                                                                       ctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGc                                                                                                                                                 <  7:342839/149‑1 (MQ=255)    ggCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGt                                                                                                                                               <  3:125573/148‑1 (MQ=255)     gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                             <  7:263590/149‑1 (MQ=255)      cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCg                                                                                                                                                <  7:179695/145‑1 (MQ=255)      cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCg                                                                                                                                                >  8:179695/1‑145 (MQ=255)         gcgcACTATAACAAGGTGCAAGAGCATGACGTATATCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTGaaa                                                                                                                                         <  7:252187/149‑3 (MQ=255)         gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGAc                                                                                                                                          >  4:161913/1‑148 (MQ=255)          cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                        <  6:246738/149‑1 (MQ=255)          cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                        <  8:287226/149‑1 (MQ=255)          cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                        <  4:359504/149‑1 (MQ=255)            cacaATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGt                                                                                                                                      >  2:194659/1‑149 (MQ=255)             acaATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTc                                                                                                                                     >  1:89481/1‑149 (MQ=255)             acaATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTc                                                                                                                                     >  1:94454/1‑149 (MQ=255)              caATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGATGGCAGGCGTTTGACAAGTcc                                                                                                                                    >  7:310353/1‑149 (MQ=255)              caATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTACCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTcc                                                                                                                                    <  3:286365/149‑1 (MQ=255)                  aaCAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgc                                                                                                                                <  3:254772/149‑1 (MQ=255)                   aCAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgc                                                                                                                                <  6:294509/148‑1 (MQ=255)                    cAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgcgc                                                                                                                              >  7:340235/1‑149 (MQ=255)                    cAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgcgc                                                                                                                              >  4:250585/1‑149 (MQ=255)                       ggTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGtt                                                                                                                           <  2:89481/149‑1 (MQ=255)                       ggTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGtt                                                                                                                           <  2:94454/149‑1 (MQ=255)                           cAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca                                                                                                                       <  3:250585/149‑1 (MQ=255)                            aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccag                                                                                                                      >  5:12620/1‑149 (MQ=255)                             agagCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagc                                                                                                                     >  2:61853/1‑149 (MQ=255)                                   tGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑g                                                                                                             <  6:140166/149‑1 (MQ=255)                                    gACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCAcgc‑‑g                                                                                                             >  4:128996/1‑148 (MQ=255)                                    gACGTATTTCACTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTg                                                                                                                                            <  6:34580/119‑1 (MQ=255)                                    gACGTATTTCACTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTg                                                                                                                                            >  5:34580/1‑119 (MQ=255)                                    gACGTATTTAATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGACGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑g                                                                                                             >  4:106603/1‑148 (MQ=255)                                     aCGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑gg                                                                                                            <  3:141665/148‑1 (MQ=255)                                        tATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGt                                                                                                        >  3:286580/1‑149 (MQ=255)                                        tATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGt                                                                                                        <  5:270054/149‑1 (MQ=255)                                          tttCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTg                                                                                                       <  5:276204/148‑1 (MQ=255)                                              aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagc                                                                                                                     <  1:197470/132‑1 (MQ=255)                                              aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagc                                                                                                                     >  2:197470/1‑132 (MQ=255)                                                 gATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGa                                                                                                <  3:161913/148‑1 (MQ=255)                                                   tGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAACTGGGCGTGATc                                                                                              >  5:274193/1‑148 (MQ=255)                                                     ttCCAACCATGGCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGaa                                                                                           <  6:104479/149‑1 (MQ=255)                                                       ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAAc                                                                                          <  2:298245/148‑1 (MQ=255)                                                       ccAACAATGTCGGACTTGGTGCATGACATTGGCTTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCAGCGCGCGTTGCCAGCCACAC‑‑GGAAGTGGGCGTGATCGAACt                                                                                         >  8:194856/1‑149 (MQ=255)                                                           ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATg                                                                                     >  6:313226/1‑149 (MQ=255)                                                           ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATg                                                                                     <  2:366687/149‑1 (MQ=255)                                                            cATGGCGGACTTGGTGCATTTCATGTGCGTAGCGGCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGAGGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGc                                                                                    <  2:351144/149‑1 (MQ=255)                                                                tCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGt                                                                                >  6:336529/1‑149 (MQ=255)                                                                  ggACTTGGTGCACGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTc                                                                               <  4:122996/148‑1 (MQ=255)                                                                    aCTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑TGAAGTGGGCGTGATCGAACTGATGCAGGTCGc                                                                             >  3:65240/1‑148 (MQ=255)                                                                    aCTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGcc                                                                            >  7:148939/1‑149 (MQ=255)                                                                    aCTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGAAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGc                                                                             >  3:60080/1‑148 (MQ=255)                                                                         gTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACa                                                                        >  1:266545/1‑148 (MQ=255)                                                                           gCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATGTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                     <  4:306473/149‑1 (MQ=255)                                                                                   ttGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACTAAAGCCGCTAt                                                             <  6:274193/149‑1 (MQ=255)                                                                                   ttGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTAt                                                             <  6:12620/149‑1 (MQ=255)                                                                                   ttGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTAt                                                             >  8:111582/1‑149 (MQ=255)                                                                                     ggCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGc                                                           >  3:49467/1‑149 (MQ=255)                                                                                         tAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTc                                                       <  2:154297/149‑1 (MQ=255)                                                                                           gcgcCGTTTATTGGAGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCa                                                      >  6:81853/1‑148 (MQ=255)                                                                                               cGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATAc                                                 <  5:313226/149‑1 (MQ=255)                                                                                                       ggCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACg                                          <  5:81853/148‑1 (MQ=255)                                                                                                       ggCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACg                                          <  8:340235/148‑1 (MQ=255)                                                                                                        gCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGc                                        >  1:155033/1‑149 (MQ=255)                                                                                                         cGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGcc                                       <  1:194659/149‑1 (MQ=255)                                                                                                         cGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGc                                        >  3:319010/1‑148 (MQ=255)                                                                                                             cTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGCTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCAccc                                   >  5:141643/1‑149 (MQ=255)                                                                                                                gCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCgg                                 <  8:148939/148‑1 (MQ=255)                                                                                                                  cGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTc                                                       <  4:162505/124‑1 (MQ=255)                                                                                                                  cGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTc                                                       >  3:162505/1‑124 (MQ=255)                                                                                                                       cccTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGAGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATa                           >  1:111559/1‑149 (MQ=255)                                                                                                                           gCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCAc                     >  8:43512/1‑149 (MQ=255)                                                                                                                             gggACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACg                    >  6:42622/1‑148 (MQ=255)                                                                                                                             gggACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGc                   <  3:106603/149‑1 (MQ=255)                                                                                                                               gACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCaa                 >  3:118252/1‑149 (MQ=255)                                                                                                                                  gATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCt              >  6:152209/1‑149 (MQ=255)                                                                                                                                  gATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCAACCCGGCCAATACTGCACGCAACCt              >  2:225442/1‑149 (MQ=255)                                                                                                                                   aTTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCAGTC‑‑GGAAGTGGGCGGGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAACTACGTCAACGGGCCCCCGGCCAATACTGCACGGAACCTg             >  8:93960/1‑149 (MQ=255)                                                                                                                                   aTTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg             >  5:47768/1‑149 (MQ=255)                                                                                                                                   aTTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg             <  4:118252/149‑1 (MQ=255)                                                                                                                                   aTTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg             <  1:61853/149‑1 (MQ=255)                                                                                                                                    tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg             >  1:15941/1‑148 (MQ=255)                                                                                                                                      tCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa           >  4:338866/1‑148 (MQ=255)                                                                                                                                      tCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa           >  4:320430/1‑148 (MQ=255)                                                                                                                                        aaaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATAc        >  4:328662/1‑149 (MQ=255)                                                                                                                                           cGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGAGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGt     >  8:290630/1‑149 (MQ=255)                                                                                                                                             cTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTTCACGCAACCTGCATACGGTGa   <  8:310353/149‑1 (MQ=255)                                                                                                                                              tGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGAt  <  3:128996/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                       CTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGAT  >  NC_002947/4004002‑4004291

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