Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,144	G→T	11.5%	W37L (TGG→TTG) ‡	PP_3533 →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,004,144	0	G	T	11.5%	212.5 / 15.3	76	W37L (TGG→TTG) ‡	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base G (29/38);  new base T (4/5);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

TGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTG  >  NC_002947/4004000‑4004289                                                                                                                                                 |                                                                                                                                                   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                                                                                                                                     <  8:130677/142‑1 (MQ=255)                   tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                     >  7:130677/1‑142 (MQ=255)                   tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcg                                                                                                                              >  7:190487/1‑149 (MQ=255)                   tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATGTCAAACGCCTGCCAGCGTTTGaaatcg                                                                                                                                   >  5:364499/1‑139 (MQ=255)                   tAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATGTCAAACGCCTGCCAGCGTTTGaaatcg                                                                                                                                   <  6:364499/144‑6 (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                                                                                <  8:190487/149‑1 (MQ=255)                                                                tGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCg                                                                                <  3:124675/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:352458/1‑149 (MQ=255)                                                                                                       ttGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAAc                                         >  3:358000/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:352458/148‑1 (MQ=255)                                                                                                                                         cAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCAt        <  4:358000/148‑1 (MQ=255)                                                                                                                                           aaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATa       >  3:54523/1‑147 (MQ=255)                                                                                                                                           aaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATa       <  4:54523/147‑1 (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)                                                                                                                                                 |                                                                                                                                                   TGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTG  >  NC_002947/4004000‑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: ‑