Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,027,278	Δ1 bp	7.3%	coding (588/993 nt)	ocd →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,027,278	0	G	.	7.3%	205.6 / 10.6	69	coding (588/993 nt)	ocd	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base . (2/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

TCAAAGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGAT  >  NC_002947/5027133‑5027425                                                                                                                                                  |                                                                                                                                                     tCAAAGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAgc                                                                                                                                                   >  7:261205/1‑149 (MQ=255)   aaaGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAgcgc                                                                                                                                                 >  2:351645/1‑149 (MQ=255)   aaaGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTAGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAgcgc                                                                                                                                                 >  3:192392/1‑149 (MQ=255)    aaGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCt                                                                                                                                                <  8:261205/149‑1 (MQ=255)          tCCTTCTTGTGCGGGATATCATGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACAtt                                                                                                                                         >  3:219179/1‑149 (MQ=255)            cttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTg                                                                                                                                        <  4:257026/148‑1 (MQ=255)            cttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTg                                                                                                                                        <  1:351645/148‑1 (MQ=255)            cttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTg                                                                                                                                        <  2:113580/148‑1 (MQ=255)             ttcttGTGCGGGATATCAAGCGAGTTCGTGTATTGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTATACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGct                                                                                                                                      <  2:114302/149‑2 (MQ=255)             ttcttGTGCGGGATATCAAGCGAGTTAGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGACTTGCAGCTATCGAGACTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGct                                                                                                                                      >  1:114302/1‑148 (MQ=255)               cttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCa                                                                                                                                      <  1:141303/147‑1 (MQ=255)               cttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGct                                                                                                                                      <  6:370978/147‑2 (MQ=255)                ttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGt                                                                                                                                    >  2:182525/1‑148 (MQ=255)                  gtgCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGcttc                                                                                                                                    <  1:258500/146‑4 (MQ=255)                   tgCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCAcc                                                                                                                                >  8:106708/1‑149 (MQ=255)                   tgCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCAc                                                                                                                                 <  6:124606/148‑1 (MQ=255)                    gCGGTATATCAAGCGAGTTCGGGTATGGGCTAGTAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCt                                                                                                                               <  6:282277/149‑1 (MQ=255)                              atgCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAACGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAg                                                                                                                     <  4:219179/147‑1 (MQ=255)                              aaGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCaa                                                                                                                      <  7:106708/148‑1 (MQ=255)                               aGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGc                                                                                                                    <  1:182525/149‑1 (MQ=255)                               aGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTAACCTGCACGCCAAg                                                                                                                     >  8:83122/1‑148 (MQ=255)                                   aGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATc                                                                                                                <  8:113158/149‑1 (MQ=255)                                     ttCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCaa                                                                                                              >  4:160773/1‑149 (MQ=255)                                       cGGGTATGGTCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCGTTGCAGCTATCTAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAAc                                                                                                             <  3:233111/148‑1 (MQ=255)                                        gggTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAAcc                                                                                                            >  3:35416/1‑148 (MQ=255)                                                    ggAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCAACGGCCTGGAAGCCAATGACGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAACCAATCAACAAATTTTACATCGc                                                                                               >  5:214618/1‑149 (MQ=255)                                                     gAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCa                                                                                              <  2:53440/149‑1 (MQ=255)                                                        ttCCGCAAAAGCGCAGCAACTGGCCTTGCATCTAGCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACTCCAAGCAATCAACCAATTTTACATCGCAgg                                                                                            <  6:63724/148‑1 (MQ=255)                                                        ttCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATc                                                                                                                <  6:220771/128‑1 (MQ=255)                                                        ttCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATc                                                                                                                >  5:220771/1‑128 (MQ=255)                                                           cGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCaa                                                                                                                  <  8:107974/123‑1 (MQ=255)                                                           cGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCaa                                                                                                                  >  7:107974/1‑123 (MQ=255)                                                           cGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGACACTGCCGCAGAGTCGGTAGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACt                                                                                        >  3:262989/1‑149 (MQ=255)                                                             cAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGc                                                                                           >  3:255299/1‑144 (MQ=255)                                                             cAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGc                                                                                           <  4:255299/144‑1 (MQ=255)                                                             cAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTg                                                                                      >  6:142975/1‑149 (MQ=255)                                                                  gcgcAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGcc                                                                                 <  4:106254/149‑1 (MQ=255)                                                                    gcagcaACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCg                                                                               <  8:220501/149‑1 (MQ=255)                                                                    gcagcaACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCg                                                                               <  8:112882/149‑1 (MQ=255)                                                                    gcagcaACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCg                                                                               >  7:272191/1‑149 (MQ=255)                                                                      agcaACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTAGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCgg                                                                              >  7:333795/1‑148 (MQ=255)                                                                           cTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCa                                                                        >  1:88591/1‑149 (MQ=255)                                                                              gCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCatat                                                                     <  4:35416/149‑1 (MQ=255)                                                                                   gCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGTAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCTTATCACgg                                                                >  2:300448/1‑149 (MQ=255)                                                                                    cAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCAAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATTTCACGGc                                                               >  1:110659/1‑149 (MQ=255)                                                                                      gCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTa                                                             >  4:170820/1‑149 (MQ=255)                                                                                       cTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTAt                                                            >  8:357364/1‑149 (MQ=255)                                                                                          tCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTAgggg                                                         >  4:197942/1‑149 (MQ=255)                                                                                            gAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCa                                                                                              <  2:164644/110‑1 (MQ=255)                                                                                            gAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCa                                                                                              >  1:164644/1‑110 (MQ=255)                                                                                             aGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGa                                                                                                                                                     <  8:361313/55‑1 (MQ=255)                                                                                                     cGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GCCATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACt                                               >  6:273119/1‑148 (MQ=255)                                                                                                        ccTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACGAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                            >  4:61842/1‑148 (MQ=255)                                                                                                          tGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCCAGCGCT‑GACATTGAAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAaca                                          >  4:15405/1‑148 (MQ=255)                                                                                                              aGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACaaa                                      >  4:317173/1‑148 (MQ=255)                                                                                                                   cTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAg                                <  3:170820/149‑1 (MQ=255)                                                                                                                          agagTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTc                         <  3:160773/149‑1 (MQ=255)                                                                                                                            agTCGGTCGATGAAGCAAAGTCGAGCTCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTAGGGGGGCAGACACTGAACACAAAAACGAGGttcaggct                       >  8:149722/1‑143 (MQ=255)                                                                                                                                  tcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGa                     >  1:69885/1‑145 (MQ=255)                                                                                                                                  tcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGa                     <  2:69885/145‑1 (MQ=255)                                                                                                                                    gATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTTCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTATCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAAt               <  4:262989/149‑1 (MQ=255)                                                                                                                                      tGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGa                                             <  5:150146/117‑1 (MQ=255)                                                                                                                                      tGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGa                                             >  6:150146/1‑117 (MQ=255)                                                                                                                                      tGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATag             >  4:374034/1‑149 (MQ=255)                                                                                                                                         tgCAATGTC‑AGCGCTCGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcag          >  7:300210/2‑149 (MQ=255)                                                                                                                                          gCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcag          >  4:96608/1‑148 (MQ=255)                                                                                                                                           cAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcag          >  7:64101/1‑147 (MQ=255)                                                                                                                                             aTGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACAGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGt      >  5:327803/1‑149 (MQ=255)                                                                                                                                                 cGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGAt  >  4:279199/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TCAAAGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGAT  >  NC_002947/5027133‑5027425

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