Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,609,522	C→T	5.1%	Q346* (CAA→TAA)	hyfG →	hydrogenase 4, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,609,522	0	C	T	5.1%	254.4 / 2.5	78	Q346* (CAA→TAA)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base T (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTGCGCTTTGATCACCCC  >  NC_000913/2609378‑2609668                                                                                                                                                 |                                                                                                                                                   ggTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCaa                                                                                                                                                  <  7:159191/147‑1 (MQ=255) ggTTTTGATACCGACTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCaa                                                                                                                                                  <  5:202691/147‑1 (MQ=255)  gTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAAc                                                                                                                                                >  8:111689/1‑148 (MQ=255)   ttttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACg                                                                                                                                               >  8:86750/1‑148 (MQ=255)   ttttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACg                                                                                                                                               >  5:330400/1‑148 (MQ=255)   ttttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACg                                                                                                                                               >  2:314108/1‑148 (MQ=255)        aTACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGaa                                                                                                                                          <  1:135199/148‑1 (MQ=255)          aCCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAAc                                                                                                                                        >  7:330420/1‑148 (MQ=255)            ctgCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGGCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTg                                                                                                                                      <  4:45742/146‑1 (MQ=255)            cGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGAGATTCTCAAAGAGCAACGTCTGCAACCGCTGCAACTg                                                                                                                                      >  6:316656/1‑148 (MQ=255)            cGGCTTTATGCAATTTTTCCGCGTGCGGG‑AAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTgg                                                                                                                                     >  1:412661/1‑148 (MQ=255)             ggCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTgg                                                                                                                                     >  2:17464/1‑148 (MQ=255)              gCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGt                                                                                                                                    <  7:14262/148‑1 (MQ=255)                tttATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTgc                                                                                                                                  >  2:181510/1‑148 (MQ=255)                  tATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTgcg                                                                                                                                 <  3:156751/147‑1 (MQ=255)                        aTTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATg                                                                                                                          >  1:371120/1‑148 (MQ=255)                          ttttCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcg                                                                                                                        >  5:317602/1‑148 (MQ=255)                                 cgTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAACCCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGAc                                                                                                                 <  6:297024/148‑1 (MQ=255)                                    gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGTaa                                                                                                                                                  <  7:234945/112‑1 (MQ=255)                                    gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGTaa                                                                                                                                                  <  7:241111/112‑1 (MQ=255)                                    gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGTaa                                                                                                                                                  >  8:234945/1‑112 (MQ=255)                                    gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGTaa                                                                                                                                                  >  8:241111/1‑112 (MQ=255)                                     cGGGATAAGTCGATGACGATGGCTGATTTGCTGATCGGGTCGCGTAAAACCTAGGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACgtgt                                                                                                             <  4:389264/148‑1 (MQ=255)                                        gAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCgg                                                                                                          <  2:31761/148‑1 (MQ=255)                                         aaaaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGa                                                                                                         <  6:415804/148‑1 (MQ=255)                                          aaaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAg                                                                                                        >  8:286568/1‑148 (MQ=255)                                                 aTGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTag                                                                                                 <  1:17464/148‑1 (MQ=255)                                                  tGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGaa                                                                                                                                          <  4:326070/106‑1 (MQ=255)                                                  tGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGaa                                                                                                                                          >  3:326070/1‑106 (MQ=255)                                                  tGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTa                                                                                                  >  7:441514/1‑146 (MQ=255)                                                    aCGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTagag                                                                                               <  8:330420/147‑1 (MQ=255)                                                    aCGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTagag                                                                                               <  6:317602/147‑1 (MQ=255)                                                    aCGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTagag                                                                                               <  2:412661/147‑1 (MQ=255)                                                       aTGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGAtgc                                                                                           >  1:244847/1‑148 (MQ=255)                                                               aTTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCtgcttgc                                                                                    >  5:225618/1‑147 (MQ=255)                                                               aTTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGGCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCtgcttgct                                                                                   <  1:307552/148‑1 (MQ=255)                                                               aTTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGGCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCtgcttgct                                                                                   <  8:121165/148‑1 (MQ=255)                                                                 tGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCACCGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGCCGTGTCTTAGCTGTTAGAGAGTCTGCTTGATAc                                                                                 >  2:149705/1‑148 (MQ=255)                                                                 tGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTAc                                                                                 >  3:127766/1‑148 (MQ=255)                                                                      aTCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACg                                                                                                                >  6:323959/1‑112 (MQ=255)                                                                      aTCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACg                                                                                                                <  5:323959/112‑1 (MQ=255)                                                                         gggTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAat                                                                          <  7:286568/147‑1 (MQ=255)                                                                            tCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTg                                                                                                     <  8:130782/117‑1 (MQ=255)                                                                            tCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTg                                                                                                     >  7:130782/1‑117 (MQ=255)                                                                             cgcgTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGa                                                                     <  7:111689/148‑1 (MQ=255)                                                                                gTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaac                                                                   >  3:126077/1‑147 (MQ=255)                                                                                gTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaac                                                                   >  1:324913/1‑147 (MQ=255)                                                                                  aaaaCCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaacaac                                                                >  1:440894/1‑148 (MQ=255)                                                                                   aaaCCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaacaac                                                                <  2:137124/147‑1 (MQ=255)                                                                                          cGGTCTGACTCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAg                                                        <  1:314108/148‑1 (MQ=255)                                                                                             tCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggc                                                     <  6:225618/148‑1 (MQ=255)                                                                                                gAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggcatt                                                  >  4:90494/1‑148 (MQ=255)                                                                                                 aaTCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCGGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggca                                                    <  6:330400/145‑1 (MQ=255)                                                                                                    cTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggcattggca                                              >  3:10203/1‑148 (MQ=255)                                                                                                    cTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggcattggca                                              >  3:205306/1‑148 (MQ=255)                                                                                                    cTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggca                                                    <  6:229764/142‑1 (MQ=255)                                                                                                    cTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggca                                                    >  5:229764/1‑142 (MQ=255)                                                                                                       aTTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTg                                                                                                     >  7:352035/1‑90 (MQ=255)                                                                                                       aTTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggcattggcatt                                            <  2:324913/147‑1 (MQ=255)                                                                                                       aTTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGggcattggc                                               <  7:86750/144‑1 (MQ=255)                                                                                                        ttGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTTAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATtc                                           <  2:133290/147‑1 (MQ=255)                                                                                                        ttGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATtc                                           >  2:18235/1‑147 (MQ=255)                                                                                                                  tCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGTACACCGCTCTCAGGGCATTGGCATTCTCGACCGACa                                 >  6:28928/1‑147 (MQ=255)                                                                                                                     ccgcGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATc                             >  6:156712/1‑148 (MQ=255)                                                                                                                        cgATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGcc                          <  4:205306/148‑1 (MQ=255)                                                                                                                           tatTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGt                       <  2:440894/148‑1 (MQ=255)                                                                                                                              tctcAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGAt                    >  7:218275/1‑148 (MQ=255)                                                                                                                                tctaaGAGCTACGTCTGAAAACGGTGAATCTGGTGCGCGAGTTGCGCGCCGAGGTGTCGGAGATGGTAGAGATGCTGCTTGCTACGCCTAATTTGGAAGAACGCACTCAGGGCATTGGCATTCTCGAGCGACaa                                <  6:316594/131‑1 (MQ=255)                                                                                                                                 cAAAGAGCAACGTCTGCACACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCGACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGAGTTg                 >  4:314725/1‑148 (MQ=255)                                                                                                                                 cAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACaaa                               >  5:316594/1‑134 (MQ=255)                                                                                                                                  aaaGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaaca                                                                  >  4:350770/1‑98 (MQ=255)                                                                                                                                  aaaGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaaca                                                                  <  3:350770/98‑1 (MQ=255)                                                                                                                                  aaaGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTg                 <  8:441514/147‑1 (MQ=255)                                                                                                                                   aaGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTgcg               >  4:134707/1‑148 (MQ=255)                                                                                                                                           cGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCTGAGCTGGTAGAGATGCTGCTTGCTACGCCTAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTGAGCTTTGATc       <  4:127766/148‑1 (MQ=255)                                                                                                                                            gTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTGCGCTTTGATCa      <  8:218275/148‑1 (MQ=255)                                                                                                                                                gCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTGCGCTTTGATCAcccc  >  8:196365/1‑148 (MQ=255)                                                                                                                                                gCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTGCGCTTTGATCAcccc  >  2:341393/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACAACGCACTCAGGGCATTGGCATTCTCGACCGACAAATCGCCCGTGATTTGCGCTTTGATCACCCC  >  NC_000913/2609378‑2609668

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