Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	236,395	C→G	7.5%	G147A (GGC→GCC)	argH ←	argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	236,395	0	C	G	7.5%	191.5 / 2.4	62	G147A (GGC→GCC)	argH	argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base G (3/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCGTGCTT  >  NC_002947/236251‑236542                                                                                                                                                 |                                                                                                                                                       cAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTGGAGCAGGcagag                                                                                                                                                    <  6:280271/149‑5 (MQ=255) cAGTCGACCAGGCGCTCATAGTAGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGATTCGGCTTAACGCTCTGCCTGCTCCAGCATGCCct                                                                                                                                                     >  6:216268/1‑148 (MQ=255)  aGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCctgc                                                                                                                                                   >  3:204953/1‑149 (MQ=255)       aCCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAg                                                                                                                                           <  4:4593/149‑1 (MQ=255)          aGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCg                                                                                                                                        >  4:58114/1‑149 (MQ=255)                  aTAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGAttt                                                                                                                                >  7:60857/1‑149 (MQ=255)                    aGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCg                                                                                                                              >  1:67894/1‑149 (MQ=255)                       cgcgGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGTCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGcc                                                                                                                            <  5:216268/148‑1 (MQ=255)                        gcgGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGccc                                                                                                                           <  8:128418/148‑1 (MQ=255)                              gAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGAt                                                                                                                      >  1:34131/1‑147 (MQ=255)                              gAGCATTTCGAACCACGCCAGCAGGTGGTGGACAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGCTGAAACCAGGCATGATGGTTTAGGCTTCACTCTCTGCCTGCTCCAGCAGCCACTGCT‑‑‑GCAGGCGGGTGATTTCGCCAAGGAt                                                                                                                      >  8:156075/1‑147 (MQ=255)                               aGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGGCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATc                                                                                                                     <  3:76599/147‑1 (MQ=255)                                      cGAACCACGCCAACAGGTGGTGGCCTAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCTGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGGCCTGCT‑‑‑GCAGGCGGGTGATTTGGCCCAGGATCAGGGCGAt                                                                                                             <  5:108661/148‑1 (MQ=255)                                        aaCCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTtc                                                                                                          >  2:289636/1‑149 (MQ=255)                                        aaCCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTtc                                                                                                          <  4:204953/149‑1 (MQ=255)                                         aCCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGTAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTtcg                                                                                                         >  6:27452/1‑149 (MQ=255)                                           cACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGccc                                                                                                                           >  4:318225/1‑129 (MQ=255)                                           cACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGccc                                                                                                                           <  3:318225/129‑1 (MQ=255)                                           cACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTtcgtc                                                                                                       <  2:34131/149‑1 (MQ=255)                                              gCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTcgcg                                                                                                    <  5:60222/149‑1 (MQ=255)                                                cagcagGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTcgcg                                                                                                    >  4:112144/1‑147 (MQ=255)                                                      gtggtgGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATc                                                                                                                     <  3:174165/124‑1 (MQ=255)                                                      gtggtgGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATc                                                                                                                     >  4:174165/1‑124 (MQ=255)                                                         gtggCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGc                                                                                         <  6:207530/149‑1 (MQ=255)                                                          tggCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCg                                                                                        <  4:35603/149‑1 (MQ=255)                                                            gCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGa                                                                                      >  2:71000/1‑149 (MQ=255)                                                            gCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTAACGCTCTGCCTGCTCCAGAAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTAGCCCAGGAAAAGGACGATATCGTCCCCAAGCAAAAGGCga                                                                                       >  2:288847/1‑147 (MQ=255)                                                               aaaGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGt                                                                                   <  3:58114/149‑1 (MQ=255)                                                                 aGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGATACCAGGCATGATTGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCTTCGCGCAGCCAAAGGCGGATGt                                                                                   <  7:156075/147‑1 (MQ=255)                                                                   gTGACCGGCTGCGCCGTCTGCAGGTGGGTGGAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCgg                                                                                       <  7:293738/141‑1 (MQ=255)                                                                   gTGACCGGCTGCGCCGTCTGCAGGTGGGTGGAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCgg                                                                                       >  8:293738/1‑141 (MQ=255)                                                                   gTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGt                                                                               >  5:233801/1‑149 (MQ=255)                                                                         ggCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCAc                                                                         <  3:200939/149‑1 (MQ=255)                                                                          gCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCAcc                                                                        <  1:288481/149‑1 (MQ=255)                                                                           cTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCt                                                                       <  3:92878/149‑1 (MQ=255)                                                                              cgcCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTgg                                                                     <  6:233801/148‑1 (MQ=255)                                                                                  gTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGtt                                                                <  1:71000/149‑1 (MQ=255)                                                                                       taGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGc                                                              <  5:27452/145‑1 (MQ=255)                                                                                        aGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGc                                                           >  2:65903/1‑148 (MQ=255)                                                                                                  aaCCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCATGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCTCTCAGCCAAAGGCGGATGTCGTTGGCCACCTGGTCGTTGCGGCTAGGCCCAGt                                                 <  1:262685/148‑1 (MQ=255)                                                                                                     cAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGc                                             <  7:118535/149‑1 (MQ=255)                                                                                                          atgatgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGctt                                        >  1:195733/1‑149 (MQ=255)                                                                                                           tgatgGTTTCGGCTTCACGCTCTGCGTGCTCCAGCGGGGCGGGCG‑‑‑GCGGGCGGGTGATTTCGCCCAGGATCGGGTCGAGTTGGTCGCGCGGCGCAAGGGGGATGTCGGTGGCCGCCGGGTGGTGGCGGCGGCGGGCAGTGTGCAGcttg                                       >  4:220597/1‑148 (MQ=255)                                                                                                           tgatgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGcttc                                       <  3:220597/149‑1 (MQ=255)                                                                                                             atgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGcttctt                                     >  1:115924/1‑149 (MQ=255)                                                                                                             atgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGcttctt                                     <  1:288847/149‑1 (MQ=255)                                                                                                                 tttCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCg                                 <  1:65903/149‑1 (MQ=255)                                                                                                                       cTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCgg                                <  8:125725/144‑1 (MQ=255)                                                                                                                       cTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGtagat                           >  7:125725/1‑145 (MQ=255)                                                                                                                           aCGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGa                       >  5:25232/1‑149 (MQ=255)                                                                                                                              ctctGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGc                    <  8:60857/149‑1 (MQ=255)                                                                                                                                gggcctgcTGGAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGCTGCCGATGcgg                  >  5:186216/3‑149 (MQ=255)                                                                                                                                 tgcctgcTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGTGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcg                   <  1:289636/147‑1 (MQ=255)                                                                                                                                  gcctgcTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcggtc                >  5:156133/1‑149 (MQ=255)                                                                                                                                   cctgcTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcggtcg               >  6:213809/1‑149 (MQ=255)                                                                                                                                   ccgcctgcagcagGGCCTGCTGGAGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTATGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcgg                  >  3:210376/8‑149 (MQ=255)                                                                                                                                     tgcTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCGGCCAATGGCGGATGTCGGTGGCAACCTGGTAGTTACGGCTACGCCCAGTATGCAGCTTCTTGCGGGTCATGCGGATGcggtcgt              >  1:102648/1‑147 (MQ=255)                                                                                                                                      cctgcagcagGGCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGctt                                        <  4:250051/117‑1 (MQ=255)                                                                                                                                      cctgcagcagGGCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGctt                                        >  3:250051/5‑121 (MQ=255)                                                                                                                                         ccAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAgg         >  7:194353/1‑149 (MQ=255)                                                                                                                                            gcagGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCGt      <  6:92075/149‑1 (MQ=255)                                                                                                                                            gcagGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCGt      <  6:156133/149‑1 (MQ=255)                                                                                                                                                gCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCGTGCtt  >  3:101476/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                       CAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCT‑‑‑GCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCGTGCTT  >  NC_002947/236251‑236542

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