Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,148	G→A	6.1%	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	6.1%	194.3 / 2.1	65	Q38Q (CAG→CAA)	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGG  >  NC_002947/4004002‑4004293                                                                                                                                                    |                                                                                                                                                  ctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAg                                                                                                                                                   <  7:154477/147‑1 (MQ=255)  tgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCg                                                                                                                                                >  4:377374/1‑149 (MQ=255)    ggCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGt                                                                                                                                               <  6:169705/148‑1 (MQ=255)     gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                             <  5:204134/149‑1 (MQ=255)      cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                             >  1:133027/1‑148 (MQ=255)          cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCCTGCCAGCGTTTGaaat                                                                                                                                        <  1:152136/149‑4 (MQ=255)                aTAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgc                                                                                                                                  <  7:205934/149‑1 (MQ=255)                    cAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgcgc                                                                                                                              >  6:130304/1‑149 (MQ=255)                    cAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgcgc                                                                                                                              <  2:89033/149‑1 (MQ=255)                              gagCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagc                                                                                                                     >  1:296185/1‑148 (MQ=255)                                  aTGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTc                                                                                                                <  2:133027/149‑1 (MQ=255)                                     aCGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGa                                                                                                             <  2:296185/149‑1 (MQ=255)                                          tttCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTgg                                                                                                        >  6:347367/1‑149 (MQ=255)                                          tttCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTgg                                                                                                        <  5:347367/149‑1 (MQ=255)                                              aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCCTGCCAGCGttt                                                                                                                                             <  6:220827/108‑1 (MQ=255)                                              aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCCTGCCAGCGttt                                                                                                                                             >  5:220827/1‑108 (MQ=255)                                                  aTGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATc                                                                                                <  5:130304/149‑1 (MQ=255)                                                       ccAACAATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCAGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCCTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACt                                                                                           >  4:24801/1‑149 (MQ=255)                                                        cAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACt                                                                                           <  3:377374/148‑1 (MQ=255)                                                         aaCCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGCGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGGTGCCCGCCACTCGGAAGTGGGCGGGATCGAACTGa                                                                                         >  8:339960/1‑149 (MQ=255)                                                          aCCATGTCGGACTTGGTGCATGACATTGGCGTAGCGACGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGAt                                                                                        >  1:235015/1‑149 (MQ=255)                                                           ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATg                                                                                       >  3:79952/1‑149 (MQ=255)                                                              tGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCg                                                                                    <  2:236792/149‑1 (MQ=255)                                                                  ggACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                             >  7:74973/1‑88 (MQ=255)                                                                  ggACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                             <  8:74973/88‑1 (MQ=255)                                                                   gACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGc                                                                               <  8:167121/149‑1 (MQ=255)                                                                     cTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGTGCGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGcc                                                                              >  1:298795/1‑149 (MQ=255)                                                                     cTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGcc                                                                              >  5:374498/1‑148 (MQ=255)                                                                        ggTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGa                                                                            <  4:79952/147‑1 (MQ=255)                                                                           gCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaa                                                                         >  1:302059/1‑147 (MQ=255)                                                                            cATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAAc                                                                                            <  3:214413/127‑1 (MQ=255)                                                                            cATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGGCGCCGACaaaa                                                                       <  6:383486/148‑1 (MQ=255)                                                                            cATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGACACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGGCGACAAAAg                                                                      >  1:306309/1‑149 (MQ=255)                                                                            cATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCAGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAAc                                                                                            >  4:214413/1‑127 (MQ=255)                                                                                 cATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCg                                                                                               <  1:101575/119‑1 (MQ=255)                                                                                 cATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCg                                                                                               >  2:101575/1‑119 (MQ=255)                                                                                 cATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGc                                                                  >  1:34949/1‑148 (MQ=255)                                                                                  aTTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGt                                                                                                          <  6:24817/107‑1 (MQ=255)                                                                                   ttGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTAt                                                               <  2:302059/149‑1 (MQ=255)                                                                                     ggCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATg                                                              <  2:298795/148‑1 (MQ=255)                                                                                              ccGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATa                                                    >  8:180523/1‑149 (MQ=255)                                                                                               cGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATAc                                                   <  4:280503/149‑1 (MQ=255)                                                                                               cGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATAc                                                   <  6:374498/149‑1 (MQ=255)                                                                                                gTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACg                                                  >  8:167488/1‑149 (MQ=255)                                                                                                gTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACg                                                  >  3:24478/1‑149 (MQ=255)                                                                                                  ttATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGt                                                 <  4:325112/148‑1 (MQ=255)                                                                                                     ttGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                      >  5:357085/1‑108 (MQ=255)                                                                                                     ttGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                      <  6:357085/108‑1 (MQ=255)                                                                                                     ttGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAAc                                             >  1:277912/1‑149 (MQ=255)                                                                                                     ttGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAAc                                             >  1:272735/1‑149 (MQ=255)                                                                                                                  cGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGcc                                 >  1:107515/1‑148 (MQ=255)                                                                                                                  cGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGcc                                 >  5:278100/1‑148 (MQ=255)                                                                                                                       cccTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACt                           >  7:59140/1‑149 (MQ=255)                                                                                                                            cGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCa                        <  2:235015/147‑1 (MQ=255)                                                                                                                             ggg‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGc                     <  7:339960/149‑1 (MQ=255)                                                                                                                                   gATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCgg                                   >  6:377575/1‑130 (MQ=255)                                                                                                                                   gATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCgg                                   <  5:377575/130‑1 (MQ=255)                                                                                                                                   gATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCt                >  8:118722/1‑149 (MQ=255)                                                                                                                                    ttgtCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg               >  2:244229/4‑149 (MQ=255)                                                                                                                                     tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTAt                                                               >  5:219204/1‑100 (MQ=255)                                                                                                                                     tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTAt                                                               <  6:219204/100‑1 (MQ=255)                                                                                                                                              cTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCCCCCGGCCAATACTGCACGCACCCTGCATACGGTGa     >  7:260621/1‑149 (MQ=255)                                                                                                                                              cTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACgg        >  1:311203/1‑146 (MQ=255)                                                                                                                                              cTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACgg        <  2:311203/146‑1 (MQ=255)                                                                                                                                               tGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGAt    >  6:347160/1‑149 (MQ=255)                                                                                                                                                 gCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATgg  >  2:156883/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGG‑ACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGG  >  NC_002947/4004002‑4004293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑