Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,148	G→A	10.8%	Q38Q (CAG→CAA) ‡	PP_3533 →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,004,148	0	G	A	10.8%	196.3 / 11.6	74	Q38Q (CAG→CAA) ‡	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base G (27/39);  new base A (4/4);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGG  >  NC_002947/4004000‑4004293                                                                                                                                                     |                                                                                                                                                   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                                                                                                                                 >  2:272696/1‑148 (MQ=255)                     aCAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgc                                                                                                                                 <  1:272696/148‑1 (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: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)                                                                                                                                                  gCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATgg  <  6:89910/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCA‑GCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGG  >  NC_002947/4004000‑4004293

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