Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,147	A→C	5.8%	Q38P (CAG→CCG) ‡	PP_3533 →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,004,147	0	A	C	5.8%	278.3 / 4.7	85	Q38P (CAG→CCG) ‡	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base C (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATG  >  NC_002947/4004002‑4004292                                                                                                                                                  |                                                                                                                                                    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)                                                                                                                                                gCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATg  <  7:43512/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTC‑‑GGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATG  >  NC_002947/4004002‑4004292

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