Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,717,410	A→C	6.1%	Q34P (CAA→CCA)	ghrB →	glyoxylate/hydroxypyruvate reductase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,717,410	0	A	C	6.1%	200.3 / 4.8	66	Q34P (CAA→CCA)	ghrB	glyoxylate/hydroxypyruvate reductase B
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base C (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGAT  >  NC_000913/3717269‑3717558                                                                                                                                              |                                                                                                                                                     tCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGt                                                                                                                                                <  5:231295/148‑1 (MQ=255) tCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTc                                                                                                                                               >  5:16967/1‑149 (MQ=255)             gATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATg                                                                                                                                     >  7:160365/1‑147 (MQ=255)              aTCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATgca                                                                                                                                   <  8:163921/148‑1 (MQ=255)                 gCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCaa                                                                                                                               >  5:116479/1‑149 (MQ=255)                          gTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGcc                                                                                                                       >  8:229858/1‑148 (MQ=255)                          gTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGcc                                                                                                                       >  8:172304/1‑148 (MQ=255)                          gTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGcc                                                                                                                       >  3:57427/1‑148 (MQ=255)                             aaTGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAg                                                                                                                   <  5:210370/149‑1 (MQ=255)                                  agagaAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCCCCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTg                                                                                                                >  5:123307/1‑147 (MQ=255)                                   gagaAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAg                                                                                                             <  3:37363/149‑1 (MQ=255)                                    agaAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAgg                                                                                                            <  4:197578/149‑1 (MQ=255)                                       aGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGttt                                                                                                         <  5:34710/149‑1 (MQ=255)                                           tGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTg                                                                                                     >  4:94328/1‑149 (MQ=255)                                           tGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTg                                                                                                     >  4:90931/1‑149 (MQ=255)                                             aaGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTgg                                                                                                    <  6:16967/148‑1 (MQ=255)                                               ggcgtccgtTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCTCAAACCGTCGAACATAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGtt                                                                                                 <  6:4952/147‑1 (MQ=255)                                               gCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGtt                                                                                                 <  7:229858/149‑1 (MQ=255)                                                     cgtTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAAc                                                                                            <  7:17849/148‑1 (MQ=255)                                                     cgtTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAAc                                                                                            <  8:106622/148‑1 (MQ=255)                                                          tCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAAt                                                                                      >  1:92212/1‑149 (MQ=255)                                                            ctctACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATg                                                                                     >  5:237657/1‑148 (MQ=255)                                                             tctACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGATTGt                                                                                    >  7:198022/1‑148 (MQ=255)                                                                aCAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAAt                                                                                >  7:248401/1‑149 (MQ=255)                                                                 cAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATg                                                                               >  1:92693/1‑149 (MQ=255)                                                                  aaaGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGc                                                                              >  1:5962/1‑149 (MQ=255)                                                                        ttACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCAGCAtt                                                                        >  1:202654/1‑149 (MQ=255)                                                                            cTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTg                                                                    <  1:136826/149‑1 (MQ=255)                                                                              gatgatTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTgg                                                                   >  4:55340/1‑148 (MQ=255)                                                                                   tttACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTg                                                                    <  1:136458/142‑1 (MQ=255)                                                                                   tttACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTg                                                                    >  2:136458/1‑142 (MQ=255)                                                                                   tttACTGCAACGCCTGCAAGAGCATGTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAAGGGGAATGTAACTGCCGCATTGCTGGaaaaaa                                                             >  6:186678/1‑149 (MQ=255)                                                                                          cAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGa                                                       >  2:7181/1‑148 (MQ=255)                                                                                               ccTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGc                                                 <  3:94328/149‑1 (MQ=255)                                                                                               ccTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGc                                                 <  3:90931/149‑1 (MQ=255)                                                                                                cTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           >  5:151015/1‑138 (MQ=255)                                                                                                cTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           <  6:151015/138‑1 (MQ=255)                                                                                                 tGCAAGAGCATTTCCCCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAA‑TTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTg                                              >  2:160359/1‑149 (MQ=255)                                                                                                  gCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTg                                              >  6:36983/1‑149 (MQ=255)                                                                                                   cAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGATTGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGt                                               <  8:198022/147‑1 (MQ=255)                                                                                                    aaGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCCTTGCTGGAAAACATGCCGAAACTGCGTGc                                             >  5:62701/1‑148 (MQ=255)                                                                                                    aaGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGc                                             >  3:60764/1‑148 (MQ=255)                                                                                                    aaGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGAAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGca                                            >  6:63583/1‑148 (MQ=255)                                                                                                      gagCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCca                                           >  4:82193/1‑148 (MQ=255)                                                                                                      gagCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCca                                           >  5:209773/1‑148 (MQ=255)                                                                                                      gagCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCca                                           >  4:85548/1‑148 (MQ=255)                                                                                                         cATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATc                                       <  6:209773/149‑1 (MQ=255)                                                                                                          aTTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           <  7:85938/128‑1 (MQ=255)                                                                                                          aTTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           >  8:85938/1‑128 (MQ=255)                                                                                                              cACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGa                                  >  8:154184/1‑149 (MQ=255)                                                                                                                ccGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGAtc                                <  8:248401/149‑1 (MQ=255)                                                                                                                   ttCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTcc                             >  7:107268/1‑149 (MQ=255)                                                                                                                      aCCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTc                          >  1:67128/1‑149 (MQ=255)                                                                                                                          ggTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCt                      <  5:36983/149‑1 (MQ=255)                                                                                                                           gTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTa                     <  5:63583/149‑1 (MQ=255)                                                                                                                            tGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTAt                    <  7:197002/149‑1 (MQ=255)                                                                                                                              gCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTAACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGa                  <  6:153288/149‑1 (MQ=255)                                                                                                                                 aaCCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACaa               <  3:109990/149‑1 (MQ=255)                                                                                                                                  aCCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAAt              <  7:172304/149‑1 (MQ=255)                                                                                                                                   ccTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCTACATCAACGATCTCCGTCGGCTATGACAAtt             >  5:88755/1‑149 (MQ=255)                                                                                                                                    cTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAAttt            <  8:107268/149‑1 (MQ=255)                                                                                                                                    atCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAAttt            >  7:175678/2‑149 (MQ=255)                                                                                                                                        gCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGAt        <  8:160365/149‑1 (MQ=255)                                                                                                                                        gCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGAt        <  2:92693/149‑1 (MQ=255)                                                                                                                                          ccACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTTATGt      <  3:85548/149‑1 (MQ=255)                                                                                                                                          ccACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGt      <  3:82193/149‑1 (MQ=255)                                                                                                                                            acaAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGTTTTAGTGGTTTCAAACGAGAATGTAAATGCCGCATTCCTGAAAAAAATGCCGAAACGGCGTCCCACAACAACGATGTCCTGCGGCTATGACAATTTGGATGTc     >  8:52478/1‑148 (MQ=255)                                                                                                                                            acaAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTCCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGACGGCTATGACAAATTTGATGTc     >  3:73244/1‑148 (MQ=255)                                                                                                                                              aaaCCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCCTCGGCTATGACAATTTTGATGTCGAt  <  6:123307/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     TCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGAT  >  NC_000913/3717269‑3717558

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