Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,143	T→C	10.6%	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	10.6%	201.2 / 11.0	75	W37R (TGG→CGG) ‡	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base T (29/38);  new base C (4/4);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACTGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTG  >  NC_002947/4003995‑4004289                                                                                                                                                     |                                                                                                                                                    cGCACTGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCt                                                                                                                                                     <  1:126854/149‑1 (MQ=255)      tgctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGAAg                                                                                                                                                >  3:111942/1‑149 (MQ=255)       gctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGGCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAgg                                                                                                                                               >  8:69005/1‑149 (MQ=255)        ctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGc                                                                                                                                              >  3:363656/1‑149 (MQ=255)            gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGttt                                                                                                                                          <  5:5560/149‑1 (MQ=255)             cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                          >  8:54443/1‑148 (MQ=255)             cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATGTCAAACGCCTGCCAGCGTTTg                                                                                                                                         <  7:215948/149‑1 (MQ=255)              aaGCGCACAATAACAAGGTGCTAGAGCATGACGTTTTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTTTTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGa                                                                                                                                        >  6:101792/1‑149 (MQ=255)              aaGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGa                                                                                                                                        <  2:94627/149‑1 (MQ=255)                gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGaaa                                                                                                                                      <  4:76878/149‑3 (MQ=255)                        tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                     >  7:130677/1‑142 (MQ=255)                        tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                     <  8:130677/142‑1 (MQ=255)                        tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcg                                                                                                                              >  7:190487/1‑149 (MQ=255)                        tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATGTCAAACGCCTGCCAGCGTTTGaaatcg                                                                                                                                   <  6:364499/144‑6 (MQ=255)                        tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATGTCAAACGCCTGCCAGCGTTTGaaatcg                                                                                                                                   >  5:364499/1‑139 (MQ=255)                         aaCAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgc                                                                                                                             >  4:251135/1‑149 (MQ=255)                          aCAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgc                                                                                                                             <  1:272696/148‑1 (MQ=255)                          aCAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgc                                                                                                                             >  2:272696/1‑148 (MQ=255)                           cAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgcgc                                                                                                                           >  1:164424/1‑149 (MQ=255)                               gTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCGTGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTg                                                                                                                       <  7:180935/149‑1 (MQ=255)                                  cAAGGGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcc                                                                                                                     <  4:35631/148‑1 (MQ=255)                                  cAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca                                                                                                                    <  8:31375/149‑1 (MQ=255)                                   aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca‑g                                                                                                                  >  1:319130/1‑149 (MQ=255)                                     gagCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca‑gc                                                                                                                 >  7:265182/1‑148 (MQ=255)                                        cATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTGaaatc                                                                                                                                    <  8:309657/127‑5 (MQ=255)                                           gACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCg                                                                                                           >  7:264433/1‑148 (MQ=255)                                            aCGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCgg                                                                                                          <  3:31007/148‑1 (MQ=255)                                               tATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGt                                                                                                      <  5:296310/149‑1 (MQ=255)                                                 tttCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTgg                                                                                                    >  5:288747/1‑149 (MQ=255)                                                   tCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTggg                                                                                                   <  3:260836/148‑1 (MQ=255)                                                    cATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGc                                                                                                  <  3:53461/148‑1 (MQ=255)                                                    cATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCg                                                                                                 >  3:157974/1‑149 (MQ=255)                                                    cATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGAGCCGTTTATTGGCGAGCTTGCCGATGACCTGCGGGACGATTTCAAACTCTGGCAGGCGTTTGACAATTACGCGCGCGTTTAAA‑GCCACTAGTAAGTGTGCg                                                                                                 >  7:298814/1‑149 (MQ=255)                                                        gATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGa                                                                                              <  7:69005/148‑1 (MQ=255)                                                          tGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATc                                                                                            <  4:363656/148‑1 (MQ=255)                                                             tCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAAc                                                                                        >  2:92533/1‑149 (MQ=255)                                                              ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCACGCCACTCGGAAGTGGGCGTGATCGaa                                                                                         <  1:4158/148‑1 (MQ=255)                                                                 aCCATGTCCGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGAt                                                                                    >  4:124675/1‑149 (MQ=255)                                                                  ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATg                                                                                   >  1:315681/1‑149 (MQ=255)                                                                     tGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCg                                                                                <  3:124675/149‑1 (MQ=255)                                                                     tGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCg                                                                                <  8:190487/149‑1 (MQ=255)                                                                         ggACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTc                                                                             <  7:114222/148‑1 (MQ=255)                                                                           aCTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGcc                                                                          <  5:347885/149‑1 (MQ=255)                                                                               ggTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGa                                                                        <  2:164424/147‑1 (MQ=255)                                                                                  gCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaa                                                                     >  3:235526/1‑147 (MQ=255)                                                                                   cATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAg                                                                  >  8:163379/1‑149 (MQ=255)                                                                                       aCATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGc                                                              <  8:265182/149‑1 (MQ=255)                                                                                               gTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCgg                                                                               <  7:28234/124‑1 (MQ=255)                                                                                               gTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCgg                                                                               >  8:28234/1‑124 (MQ=255)                                                                                               gTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCtt                                                      <  8:264433/149‑1 (MQ=255)                                                                                                tAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTc                                                     <  7:264968/149‑1 (MQ=255)                                                                                                          tATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTc                                            >  5:89910/1‑148 (MQ=255)                                                                                                            ttGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAAc                                         >  3:358000/1‑149 (MQ=255)                                                                                                            ttGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAAc                                         >  3:352458/1‑149 (MQ=255)                                                                                                                  aGCTTGCCGCTGCCCTGCGGGACGAATTCAAACGCTGGCAGGCGTTTGACAAGCCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCAc                                   <  5:351720/149‑1 (MQ=255)                                                                                                                        ccGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGc                              <  6:288747/148‑1 (MQ=255)                                                                                                                              cccTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATAc                        <  4:157974/148‑1 (MQ=255)                                                                                                                               ccTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATCCCTTCAAATACGTCAACGGCCACCCGGCCAATACTg                      <  5:101792/149‑1 (MQ=255)                                                                                                                                   cGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCa                    <  2:303446/147‑1 (MQ=255)                                                                                                                                   cGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCAc                   <  7:163379/148‑1 (MQ=255)                                                                                                                                    gggACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACg                  >  8:22841/1‑148 (MQ=255)                                                                                                                                      gACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCaa               <  2:319130/149‑1 (MQ=255)                                                                                                                                         gATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCt            >  1:237378/1‑149 (MQ=255)                                                                                                                                           tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGc          <  7:54443/149‑1 (MQ=255)                                                                                                                                            gtCAAACGNNTGCCNGCGTNTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAcc             <  6:355118/144‑1 (MQ=255)                                                                                                                                            gtCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCttga         >  5:355118/2‑146 (MQ=255)                                                                                                                                            gtCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa         >  3:348393/2‑149 (MQ=255)                                                                                                                                            gtCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa         >  3:235932/2‑149 (MQ=255)                                                                                                                                              cAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCAt        <  4:358000/148‑1 (MQ=255)                                                                                                                                              cAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCAt        <  4:352458/148‑1 (MQ=255)                                                                                                                                                aaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATa       <  4:54523/147‑1 (MQ=255)                                                                                                                                                aaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATa       >  3:54523/1‑147 (MQ=255)                                                                                                                                                aaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACg     >  3:190629/1‑149 (MQ=255)                                                                                                                                                aaCGCTGGCAGGCGTTTGACAAGTCCGCGAGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACg     <  5:366447/149‑1 (MQ=255)                                                                                                                                                 aCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACgg    <  2:356836/149‑1 (MQ=255)                                                                                                                                                  cGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGt   <  3:45325/149‑1 (MQ=255)                                                                                                                                                   gCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTg  <  2:237378/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGCACTGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTG  >  NC_002947/4003995‑4004289

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