Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,671,075	A→G	6.1%	K495R (AAG→AGG)	treF →	cytoplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,671,075	0	A	G	6.1%	204.7 / 3.7	65	K495R (AAG→AGG)	treF	cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (20/41);  new base G (2/2);  total (22/43)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAATTGGT  >  NC_000913/3670928‑3671223                                                                                                                                                    |                                                                                                                                                     cTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAg                                                                                                                                                     >  3:15010/1‑149 (MQ=255)  tGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGa                                                                                                                                                    <  3:504534/149‑1 (MQ=255)      gggATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGt                                                                                                                                                <  2:189561/149‑1 (MQ=255)          ttCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAAt                                                                                                                                            <  3:287631/149‑1 (MQ=255)            cTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTAAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCa                                                                                                                                          <  1:408839/149‑1 (MQ=255)              ggCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGt                                                                                                                                        <  1:456362/149‑1 (MQ=255)               gCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGGCGGTGAATCAGtt                                                                                                                                       >  2:24163/1‑149 (MQ=255)               gCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGtt                                                                                                                                       >  4:55195/1‑149 (MQ=255)               gCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGtt                                                                                                                                       <  1:327900/149‑1 (MQ=255)               gCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGtt                                                                                                                                       <  4:242232/149‑1 (MQ=255)                cAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTc                                                                                                                                      <  1:197258/149‑1 (MQ=255)                cAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTc                                                                                                                                      <  3:489786/149‑1 (MQ=255)                 aaGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCt                                                                                                                                     <  4:298786/149‑1 (MQ=255)                 aaGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCt                                                                                                                                     <  4:298526/149‑1 (MQ=255)                     gAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGAGCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCt                                                                                                                                 <  3:112350/149‑1 (MQ=255)                      aGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTg                                                                                                                                <  1:232177/149‑1 (MQ=255)                       gTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTgg                                                                                                                               >  4:541730/1‑149 (MQ=255)                                 ggTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACaa                                                                                                                     >  4:10920/1‑149 (MQ=255)                                  gTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACaaa                                                                                                                    >  1:77346/1‑149 (MQ=255)                                   tGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAAc                                                                                                                   <  4:15010/149‑1 (MQ=255)                                   tGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAAc                                                                                                                   <  2:77346/149‑1 (MQ=255)                                                  aCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAcc                                                                                                    >  1:474119/1‑149 (MQ=255)                                                  aCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAcc                                                                                                    <  4:398407/149‑1 (MQ=255)                                                    ccAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAGTGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCat                                                                                                  >  3:387661/1‑149 (MQ=255)                                                    ccAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCat                                                                                                  >  1:293244/1‑149 (MQ=255)                                                     cAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCata                                                                                                 <  2:236916/149‑1 (MQ=255)                                                     cAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCata                                                                                                 <  2:62065/149‑1 (MQ=255)                                                      aaCGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAGGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCatat                                                                                                >  1:511102/1‑149 (MQ=255)                                                      aaCGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAGGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCatat                                                                                                >  1:116275/1‑149 (MQ=255)                                                             gggCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGAt                                                                                         <  4:542385/149‑1 (MQ=255)                                                              ggCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATg                                                                                        <  1:278532/149‑1 (MQ=255)                                                                 aCCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGgtg                                                                                     >  3:165321/1‑149 (MQ=255)                                                                            ggatggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGaa                                                                          >  1:199658/1‑149 (MQ=255)                                                                              atggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAgg                                                                        <  2:120604/149‑1 (MQ=255)                                                                                  cGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggcggc                                                                    <  3:10920/149‑1 (MQ=255)                                                                                     ttCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGt                                                                 <  2:293244/149‑1 (MQ=255)                                                                                      tCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTg                                                                <  3:33569/149‑1 (MQ=255)                                                                                       cAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGATGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTgg                                                               <  2:488910/149‑1 (MQ=255)                                                                                        aGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGc                                                              <  4:251895/149‑1 (MQ=255)                                                                                        aGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGc                                                              >  2:426359/1‑149 (MQ=255)                                                                                        aGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGc                                                              <  2:183436/149‑1 (MQ=255)                                                                                        aGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGc                                                              >  1:183436/1‑149 (MQ=255)                                                                                          ggATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACCGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGa                                                            <  4:387661/149‑1 (MQ=255)                                                                                          ggATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGa                                                            >  3:96707/1‑149 (MQ=255)                                                                                                      tACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCa                                                <  2:103837/149‑1 (MQ=255)                                                                                                         ggCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGa                                             >  2:406101/1‑149 (MQ=255)                                                                                                         ggCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGa                                             >  1:540133/1‑149 (MQ=255)                                                                                                          gCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGAt                                            <  2:313047/149‑1 (MQ=255)                                                                                                            gATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATgg                                          >  3:392855/1‑149 (MQ=255)                                                                                                               gACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGtt                                       <  3:527055/149‑1 (MQ=255)                                                                                                                aCCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGttt                                      >  1:371210/1‑149 (MQ=255)                                                                                                                   ttCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGc                                   <  2:199658/149‑1 (MQ=255)                                                                                                                    tCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCt                                  <  3:541730/149‑1 (MQ=255)                                                                                                                     cTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTg                                 <  3:55195/149‑1 (MQ=255)                                                                                                                     cTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTg                                 <  1:24163/149‑1 (MQ=255)                                                                                                                        ggTGATGAAATCGCGCGAAGCTGGCTGAGGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGAc                              <  2:511102/149‑1 (MQ=255)                                                                                                                        ggTGATGAAATCGCGCGAAGCTGGCTGAGGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGAc                              <  2:116275/149‑1 (MQ=255)                                                                                                                             tgaAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAATATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACg                         <  4:336338/149‑1 (MQ=255)                                                                                                                             tgaAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACg                         <  4:61895/149‑1 (MQ=255)                                                                                                                             tgaAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACg                         <  3:566191/149‑1 (MQ=255)                                                                                                                                      gcgAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTAc                >  2:12286/1‑149 (MQ=255)                                                                                                                                           gctggctgAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGt           >  4:501963/1‑149 (MQ=255)                                                                                                                                                ctgAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAAt      <  4:392855/149‑1 (MQ=255)                                                                                                                                                  gAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAATTg    <  4:165321/149‑1 (MQ=255)                                                                                                                                                  gAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAATTg    <  3:101693/149‑1 (MQ=255)                                                                                                                                                  caaGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCCAAAAATACCGTATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAATTg    <  3:50976/148‑1 (MQ=255)                                                                                                                                                   aaGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAATTgg   <  2:375656/149‑1 (MQ=255)                                                                                                                                                    aGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAATTGGt  <  4:137630/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGGTGGCGAGTATCCGTTGCAGGATGGGTTTGGCTGGACTAACGGTGTGGTACGCCGTTTAATTGGT  >  NC_000913/3670928‑3671223

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