Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,927,030	G→T	5.1%	Q604K (CAG→AAG)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,927,030	0	G	T	5.1%	276.3 / 3.9	79	Q604K (CAG→AAG)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base G (33/41);  new base T (2/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACC  >  NC_000913/1926883‑1927175                                                                                                                                                    |                                                                                                                                                  cATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtg                                                                                                                                                   >  2:248727/1‑148 (MQ=255)  aTGGCTTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtgt                                                                                                                                                  <  1:144417/148‑1 (MQ=255)  aTGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtgt                                                                                                                                                  <  2:330476/148‑1 (MQ=255)    ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGTCTGTGc                                                                                                                                                >  4:86525/1‑148 (MQ=255)    ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGc                                                                                                                                                >  2:229111/1‑148 (MQ=255)          aGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg                                                                                                                                            <  8:129245/146‑1 (MQ=255)          aGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg                                                                                                                                            <  1:248727/146‑1 (MQ=255)          aGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGa                                                                                                                                           <  1:229111/147‑1 (MQ=255)           gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACg                                                                                                                                         >  8:401848/1‑148 (MQ=255)           gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACg                                                                                                                                         <  6:29972/148‑1 (MQ=255)           gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATAGGTTTTCAGCTAGCGCAATTTAGCGACCAAATTTGCCGGTTCCCAATATTGCACCTGAGAATAGTTCAAACCGGTCGTTACCAGTAAAAGCGGAAAAGACTGTGCGGTGACg                                                                                                                                         >  5:118274/1‑148 (MQ=255)               aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGt                                                                                                                                     <  8:165062/148‑1 (MQ=255)                tGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTc                                                                                                                                    <  7:126334/148‑1 (MQ=255)                 gTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCa                                                                                                                                   >  4:313315/1‑148 (MQ=255)                 gTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAATTAGCCTGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTCTGCGGTGACGTTGTc                                                                                                                                    >  1:339871/1‑147 (MQ=255)                  tCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCAt                                                                                                                                  <  8:76733/148‑1 (MQ=255)                     gTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACg                                                                                                                               >  1:349885/1‑148 (MQ=255)                           agCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCAt                                                                                                                                  <  8:304382/139‑1 (MQ=255)                           agCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCAt                                                                                                                                  >  7:304382/1‑139 (MQ=255)                           agCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGt                                                                                                                         <  1:15806/148‑1 (MQ=255)                              aaTAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACggg                                                                                                                             <  5:188760/141‑1 (MQ=255)                              aaTAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACggg                                                                                                                             >  6:188760/1‑141 (MQ=255)                                 aaaaGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCttt                                                                                                                   >  6:279355/1‑148 (MQ=255)                                 aaaaGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCttt                                                                                                                   >  6:208602/1‑148 (MQ=255)                                  aaaGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCtttt                                                                                                                  >  8:346993/1‑148 (MQ=255)                                   aaGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTc                                                                                                                 <  7:346993/148‑1 (MQ=255)                                   aaGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTc                                                                                                                 <  5:347782/148‑1 (MQ=255)                                    aGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCa                                                                                                                <  5:208602/148‑1 (MQ=255)                                        ggTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTa                                                                                                            >  1:90994/1‑148 (MQ=255)                                         gTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTAc                                                                                                           <  7:401848/148‑1 (MQ=255)                                               ggTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTa                                                                                                      <  5:117331/147‑1 (MQ=255)                                                gTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGTCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAat                                                                                                    <  3:86525/148‑1 (MQ=255)                                                gTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAat                                                                                                    <  3:313315/148‑1 (MQ=255)                                                      aGCTCGCGCAATTTAGGGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGATTCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGc                                                                                               <  2:280931/147‑1 (MQ=255)                                                           tcgaattttAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCTGTTAAGCCTGTGCGGTGACGTTTTCATACTGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGAt                                                                                           <  6:118274/140‑1 (MQ=255)                                                           gcgcAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCTGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATTTTGCGGAt                                                                                           <  2:349885/146‑1 (MQ=255)                                                           gcgcAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGAt                                                                                           <  2:90994/146‑1 (MQ=255)                                                              cAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGc                                                                                      <  5:279355/148‑1 (MQ=255)                                                                           cATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCAct                                                                         >  2:82560/1‑148 (MQ=255)                                                                              tttGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTa                                                                                                      >  5:373647/1‑116 (MQ=255)                                                                              tttGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTa                                                                                                      <  6:373647/116‑1 (MQ=255)                                                                              tttGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCtt                                                                      >  3:298617/1‑148 (MQ=255)                                                                              tttGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCtt                                                                      >  3:263441/1‑148 (MQ=255)                                                                              tttGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCtt                                                                      >  3:109376/1‑148 (MQ=255)                                                                               ttGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGa                                                                                            <  2:128596/125‑1 (MQ=255)                                                                               ttGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGa                                                                                            >  1:128596/1‑125 (MQ=255)                                                                                  ccGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaa                                                                   <  4:109376/147‑1 (MQ=255)                                                                                      ttCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCa                                                                           >  6:316936/1‑135 (MQ=255)                                                                                      ttCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCa                                                                           <  5:316936/135‑1 (MQ=255)                                                                                        cccAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAc                                                             <  1:82560/147‑1 (MQ=255)                                                                                         ccAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAcc                                                            <  7:16808/147‑1 (MQ=255)                                                                                         ccAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAcc                                                            <  4:263441/147‑1 (MQ=255)                                                                                           aaTATTGCACCTGTGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATAtt                                                                                                 <  3:150732/108‑1 (MQ=255)                                                                                           aaTATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATAtt                                                                                                 >  4:150732/1‑108 (MQ=255)                                                                                           aaTATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaaa                                                                  >  3:56877/1‑139 (MQ=255)                                                                                           aaTATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaaa                                                                  <  4:56877/139‑1 (MQ=255)                                                                                              atTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTgg                                                       >  4:110230/1‑147 (MQ=255)                                                                                              atTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTgg                                                       >  6:194060/1‑147 (MQ=255)                                                                                                       gagaATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAAt                                             >  8:218482/1‑148 (MQ=255)                                                                                                            tCGTGCAAACCGGTCGTTACCAGTAAATGCGGATCAGCCTGCGCGGGGACGCTGCCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGAtt                                        <  4:298617/148‑1 (MQ=255)                                                                                                             cGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTTTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTAcccc                                                                            <  2:244186/111‑1 (MQ=255)                                                                                                             cGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTTTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTAcccc                                                                            >  1:244186/1‑111 (MQ=255)                                                                                                             cGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTTTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTAcccc                                                                            >  1:238324/1‑111 (MQ=255)                                                                                                             cGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTTTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTAcccc                                                                            <  2:238324/111‑1 (MQ=255)                                                                                                             cGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAcc                                                            <  2:23087/127‑1 (MQ=255)                                                                                                             cGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAcc                                                            >  1:23087/1‑127 (MQ=255)                                                                                                               tGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATtcat                                     <  7:218482/148‑1 (MQ=255)                                                                                                                               ccAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCa                                                           >  1:268405/1‑110 (MQ=255)                                                                                                                               ccAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCa                                                           <  2:268405/110‑1 (MQ=255)                                                                                                                                 aGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTacaa                   >  6:338383/1‑148 (MQ=255)                                                                                                                                   tAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTacaaca                 >  1:88126/1‑148 (MQ=255)                                                                                                                                    aaaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTacaacaa                >  2:316853/1‑147 (MQ=255)                                                                                                                                    aaaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACAt                >  2:193739/1‑148 (MQ=255)                                                                                                                                     aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCTCTCTTCAAACTCACCAGTGGTAAGAGTAATTGATTCATCAAGCATCGTTGTTACAACATc               <  2:339871/148‑1 (MQ=255)                                                                                                                                     aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATc               <  3:110230/148‑1 (MQ=255)                                                                                                                                     aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATc               <  5:194060/148‑1 (MQ=255)                                                                                                                                        gctgATAAGCCTGGGCGGTGACGTTGTCATTCGTGATGTAGAATTTCATGTACTAGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATcaa             >  8:332838/4‑147 (MQ=255)                                                                                                                                             tAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACg       >  7:394697/1‑148 (MQ=255)                                                                                                                                               aGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACgg      >  7:160274/1‑147 (MQ=255)                                                                                                                                                gCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTa    <  8:160274/148‑1 (MQ=255)                                                                                                                                                  cTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAACCTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTAcc  <  5:338383/148‑1 (MQ=255)                                                                                                                                                  cTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTAcc  <  1:316853/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACC  >  NC_000913/1926883‑1927175

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