Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	555,233	A→T	6.3%	E208V (GAA→GTA)	cysS →	cysteinyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	555,233	0	A	T	6.3%	208.5 / 3.8	63	E208V (GAA→GTA)	cysS	cysteinyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base T (2/2);  total (29/34)
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.81e-01

TCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATGGTCAGTATGTGAACTACTGGATGCAC  >  NC_000913/555087‑555378                                                                                                                                                   |                                                                                                                                                  tctcGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGaa                                                                                                                                                  >  1:21128/1‑148 (MQ=255)     gACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGaa                                                                                                                                                  <  2:59107/144‑1 (MQ=255)     gACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGGCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGa                                                                                                                                                   <  7:99417/143‑1 (MQ=255)              cagccagcCGCTCGCGTTCACGGGGTCGACGACAACCGCAACCCCATCGACTTCGTACTGTGGAAGAGGTCCAAAGAGGGCCAACCCAGCTGGCCTTCTCCGTGGGGCGCGGGTCTTCCTGGCTGGCACATTGAATGTTCGGCAAt                                                                                                                                       <  6:27020/146‑1 (MQ=255)                    ggcGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCa                                                                                                                                         >  4:25020/1‑138 (MQ=255)                    ggcGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCa                                                                                                                                         <  3:69082/138‑1 (MQ=255)                    ggcGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCa                                                                                                                                         >  4:69082/1‑138 (MQ=255)                    ggcGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCa                                                                                                                                         <  3:25020/138‑1 (MQ=255)                                 aCGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTaa                                                                                                                  <  8:28746/148‑1 (MQ=255)                                  cGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAAc                                                                                                                 >  3:123047/1‑148 (MQ=255)                                       tCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGaa                                                                                                                                    >  4:231190/1‑124 (MQ=255)                                       tCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGGCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGaa                                                                                                                                    <  3:231190/124‑1 (MQ=255)                                        cgacgaCAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGTATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACttt                                                                                                           >  7:111202/1‑148 (MQ=255)                                        cgacgaCAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCa                                                                                                                                         >  3:356078/1‑118 (MQ=255)                                        cgacgaCAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCa                                                                                                                                         <  4:356078/118‑1 (MQ=255)                                        cgacgaCAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGc                                                                                                                                <  8:295019/127‑1 (MQ=255)                                        cgacgaCAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGc                                                                                                                                >  7:295019/1‑127 (MQ=255)                                         gacgacAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACttt                                                                                                           >  2:233980/1‑147 (MQ=255)                                         gacgacAAACGCAACCCAATGGACTTCGTGCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTg                                                                                                          <  8:227239/148‑1 (MQ=255)                                           cgacAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGa                                                                                                         >  7:339982/1‑147 (MQ=255)                                             acAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGa                                                                                                                                                   <  8:312864/103‑1 (MQ=255)                                             acAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGa                                                                                                                                                   >  7:312864/1‑103 (MQ=255)                                               aaaCGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCGGGGTAACCACTTTGATATcc                                                                                                    <  5:93927/148‑1 (MQ=255)                                                aaCGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCa                                                                                                   >  2:128736/1‑148 (MQ=255)                                                       cccAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATc                                                                                                     <  2:342424/139‑1 (MQ=255)                                                       cccAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATc                                                                                                     >  1:342424/1‑139 (MQ=255)                                                            tGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGtt                                                                                        >  3:139533/1‑147 (MQ=255)                                                            tGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGtt                                                                                        >  1:169824/1‑147 (MQ=255)                                                             ggACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGtt                                                                                        >  2:313721/1‑146 (MQ=255)                                                             ggACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGtt                                                                                        >  6:372993/1‑146 (MQ=255)                                                               aCTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGa                                                                                    <  1:168927/148‑1 (MQ=255)                                                                   cGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGata                                                                                                       >  2:170606/1‑125 (MQ=255)                                                                   cGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGata                                                                                                       <  1:170606/125‑1 (MQ=255)                                                                    gTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTg                                                                                <  1:128736/147‑1 (MQ=255)                                                                       cTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATg                                                                             <  1:313721/147‑1 (MQ=255)                                                                             aaGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGc                                                                      <  1:233980/148‑1 (MQ=255)                                                                                     gAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCacga                                                               <  2:21128/147‑1 (MQ=255)                                                                                      aaaGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCacgaa                                                              >  7:108401/1‑147 (MQ=255)                                                                                           gggCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGTATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCacgaaaacgaa                                                        <  8:111202/148‑1 (MQ=255)                                                                                           gggCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCacgaaaacgaa                                                        <  5:372993/148‑1 (MQ=255)                                                                                            ggCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCCCACCacgaaaacgaaa                                                       <  4:123047/148‑1 (MQ=255)                                                                                             gCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCAcggcggc                                                                                            <  1:150257/110‑1 (MQ=255)                                                                                             gCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCAcggcggc                                                                                            >  2:150257/1‑110 (MQ=255)                                                                                               gAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGTATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATcg                                                    >  1:323327/1‑148 (MQ=255)                                                                                                aaCCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATcgc                                                   <  2:169824/148‑1 (MQ=255)                                                                                                 aCCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATcgcg                                                  <  8:339982/148‑1 (MQ=255)                                                                                                  ccGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCgcgc                                                 <  4:139533/148‑1 (MQ=255)                                                                                                    gAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGa                                                                                    <  1:371279/111‑1 (MQ=255)                                                                                                    gAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGa                                                                                    >  2:371279/1‑111 (MQ=255)                                                                                                    gAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGc                                                                      <  8:71849/125‑1 (MQ=255)                                                                                                    gAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGc                                                                      >  7:71849/1‑125 (MQ=255)                                                                                                               ctcCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCt                                        <  6:196745/144‑1 (MQ=255)                                                                                                               ctcCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCt                                        >  5:196745/1‑144 (MQ=255)                                                                                                                 ccGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGcc                                   >  7:24087/1‑147 (MQ=255)                                                                                                                    tGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCatg                               <  8:108401/148‑1 (MQ=255)                                                                                                                        gcgcGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACTAAATCGCGCAGTCCACCTGTGCCCATGATgg                           <  5:105533/148‑1 (MQ=255)                                                                                                                          gcgGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATGGt                          >  2:73620/1‑147 (MQ=255)                                                                                                                             ggTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATGGTCAGt                      <  8:24087/148‑1 (MQ=255)                                                                                                                              gtcgtcCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGc                                                                      <  5:249747/99‑1 (MQ=255)                                                                                                                              gtcgtcCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGc                                                                      >  6:249747/1‑99 (MQ=255)                                                                                                                              gtcgtcCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATGGTCAGTa                     <  1:73620/148‑1 (MQ=255)                                                                                                                                               aTTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCAc                                          >  2:1508/1‑110 (MQ=255)                                                                                                                                               aTTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCAc                                          <  1:1508/110‑1 (MQ=255)                                                                                                                                                ttGTATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATCGTCAGTATGTGAACTACTGGATGCa   <  2:323327/148‑1 (MQ=255)                                                                                                                                                 tGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTcc                                                                         <  6:29844/77‑1 (MQ=255)                                                                                                                                                 tGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTcc                                                                         >  5:29844/1‑77 (MQ=255)                                                                                                                                                 tGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATGGTCAGTATGTGAACTACTGGATGCAc  >  6:345592/1‑148 (MQ=255)                                                                                                                                                 tGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATGGTCAGTATGTGAACTACTGGATGCAc  >  6:22008/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGGCGGCGGTTCAGACCTGATGTTCCCGCACCACGAAAACGAAATCGCGCAGTCCACCTGTGCCCATGATGGTCAGTATGTGAACTACTGGATGCAC  >  NC_000913/555087‑555378

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