Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,092,672	A→G	6.0%	K265K (AAA→AAG)	gshB →	glutathione synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,092,672	0	A	G	6.0%	206.0 / 2.9	66	K265K (AAA→AAG)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand):  ref base A (38/24);  new base G (2/2);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTT  >  NC_000913/3092524‑3092814                                                                                                                                                     |                                                                                                                                               gCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGaaaaa                                                                                                                                                <  4:205181/149‑1 (MQ=255)  cGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAg                                                                                                                                               <  6:142031/149‑1 (MQ=255)    aGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAggg                                                                                                                                             <  4:113409/149‑1 (MQ=255)    aGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAggg                                                                                                                                             <  3:284035/149‑1 (MQ=255)    aGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAAGCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCTCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAggg                                                                                                                                             <  6:45775/149‑1 (MQ=255)        ggTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCGGa                                                                                                                                         >  2:295268/1‑149 (MQ=255)         gTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGAt                                                                                                                                        >  2:108888/1‑149 (MQ=255)           aCCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAGCGGCTGAttt                                                                                                                                      >  8:196258/1‑149 (MQ=255)           aCCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGAttt                                                                                                                                      >  4:146839/1‑149 (MQ=255)             cGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAAGCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGAttttt                                                                                                                                    <  6:220132/149‑1 (MQ=255)             cGTACTGCCTGGCGCCTATTCCGCAGGGGGGCGAACCCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTACGCTGACGGAAAGTGAATGGAAAATCGCCCGTCAGATCGGGCCGCCGCTGAAAGAAAAAGGGCTGAttttt                                                                                                                                    >  3:238397/1‑149 (MQ=255)               tACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAAGCTGGCTGCCGGTGGTCGCGGTGAAGCTCGTCCGCTGTCGGAAAGTGTCTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAACAAAAAGGGCTGATTTttgt                                                                                                                                  <  1:213735/149‑1 (MQ=255)                 cTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTttgttg                                                                                                                                <  2:130530/149‑1 (MQ=255)                 cTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTttgtt                                                                                                                                 <  1:125969/148‑1 (MQ=255)                  tGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTgg                                                                                                                               >  4:250253/1‑149 (MQ=255)                      tGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTg                                                                                                                           >  6:247827/1‑149 (MQ=255)                          gcgTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGata                                                                                                                       >  7:324593/1‑149 (MQ=255)                          gcgTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGata                                                                                                                       <  4:275814/149‑1 (MQ=255)                           cgTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGatat                                                                                                                      >  6:218632/1‑149 (MQ=255)                             tATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGa                                                                                                                         <  7:248482/144‑1 (MQ=255)                             tATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGa                                                                                                                         >  8:248482/1‑144 (MQ=255)                             tATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATatca                                                                                                                    <  1:108888/149‑1 (MQ=255)                               ttCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATatcatc                                                                                                                  >  6:115853/1‑149 (MQ=255)                                         gggCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAGGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTc                                                                                                        >  2:295863/1‑149 (MQ=255)                                         gggCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTAAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTc                                                                                                        >  5:211404/1‑149 (MQ=255)                                          ggCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAGGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTct                                                                                                       <  1:295863/149‑1 (MQ=255)                                           gCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGAttttt                                                                                                                                    >  2:14968/1‑119 (MQ=255)                                            cGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTctga                                                                                                     <  6:12291/149‑1 (MQ=255)                                             gAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTctgac                                                                                                    <  3:250253/149‑1 (MQ=255)                                                 cccGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTctgactga                                                                                                 >  5:16134/1‑148 (MQ=255)                                                  ccGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGaaa                                                                                               >  2:45804/1‑149 (MQ=255)                                                          aTCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCGGAAAGAAAAAGGGCTGATTTTTGTTGGTCGGGATATCATCGGCGCCCGTCTGACTGAAATTACCgg                                                                                        >  4:14161/1‑147 (MQ=255)                                                                    cGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCaa                                                                             >  1:316530/1‑149 (MQ=255)                                                                    cGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCaa                                                                             >  5:38900/1‑149 (MQ=255)                                                                             cgGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAGGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTc                                                                    <  7:196258/149‑1 (MQ=255)                                                                              ggTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCg                                                                   >  7:15897/1‑149 (MQ=255)                                                                                tGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGt                                                                  >  8:219205/1‑148 (MQ=255)                                                                                  aaCCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTga                                                                <  6:16134/148‑1 (MQ=255)                                                                                    ccTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGAt                                                             >  7:275348/1‑149 (MQ=255)                                                                                     cTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGAtt                                                            >  3:172434/1‑149 (MQ=255)                                                                                        gTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGaa                                                         <  3:222272/149‑1 (MQ=255)                                                                                        gTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGa                                                          <  3:146839/148‑1 (MQ=255)                                                                                           cGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAg                                                        >  3:286816/1‑147 (MQ=255)                                                                                            gCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAg                                                        >  1:9995/1‑146 (MQ=255)                                                                                             cTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCaga                                                    <  5:247827/149‑1 (MQ=255)                                                                                                  ggAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTc                                               >  1:275850/1‑149 (MQ=255)                                                                                                            tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGAt                                     >  2:227893/1‑149 (MQ=255)                                                                                                              gAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGAttttt                                                                                                                                    >  6:211351/1‑52 (MQ=255)                                                                                                                     gCCCGTCAGATCGGGCCGACGCTGACAGAAAAAGGGCGGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCGAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGaa                             <  2:215301/148‑1 (MQ=255)                                                                                                                      cccGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATg                           <  4:172434/149‑1 (MQ=255)                                                                                                                          tCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGCTCACCGGAATGTTaa                       >  6:127782/1‑149 (MQ=255)                                                                                                                          tCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTaa                       >  7:111049/1‑149 (MQ=255)                                                                                                                            aGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAAt                      <  2:316530/148‑1 (MQ=255)                                                                                                                             gATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATgg                    >  2:61748/1‑149 (MQ=255)                                                                                                                             gATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATg                     >  3:105111/1‑148 (MQ=255)                                                                                                                              aTCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGa                   <  6:51909/149‑1 (MQ=255)                                                                                                                              aTCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGa                   <  5:218632/149‑1 (MQ=255)                                                                                                                                     cGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATc            >  4:194957/1‑149 (MQ=255)                                                                                                                                        cGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGGCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGaa         <  8:324593/149‑1 (MQ=255)                                                                                                                                            gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc     >  6:1694/1‑149 (MQ=255)                                                                                                                                            gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc     >  1:168884/1‑149 (MQ=255)                                                                                                                                             aaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACg    >  7:233842/1‑149 (MQ=255)                                                                                                                                             aaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACg    >  3:199190/1‑149 (MQ=255)                                                                                                                                               agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAAt                      >  7:298259/1‑129 (MQ=255)                                                                                                                                               agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAAt                      <  8:298259/129‑1 (MQ=255)                                                                                                                                               agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt  >  5:65581/1‑149 (MQ=255)                                                                                                                                               agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt  >  4:9797/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               GCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGGGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTT  >  NC_000913/3092524‑3092814

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