Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,088,626	T→C	5.3%	H168R (CAT→CGT)	yeeZ ←	putative epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,088,626	0	T	C	5.3%	250.5 / 3.4	75	H168R (CAT→CGT)	yeeZ	putative epimerase
Reads supporting (aligned to +/- strand):  ref base T (42/29);  new base C (2/2);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

TGAGCTGGCGCACATATATTATAGATGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGTGCCTTGCGC  >  NC_000913/2088479‑2088774                                                                                                                                                    |                                                                                                                                                     tGAGCTGGCGCACATATATTATAGATGTGTCCGCCTTTAGGTGCCTGTAACGACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGAt                                                                                                                                                      >  2:858364/1‑148 (MQ=255)           cacaTATATTATAGATGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGACGTCATCCCgg                                                                                                                                            <  2:85211/148‑1 (MQ=255)             caTATATTATAGATGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGACGACGTCCCGGTcc                                                                                                                                         <  4:73991/149‑1 (MQ=255)               tatatTATAGATGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCg                                                                                                                                        <  2:862166/148‑1 (MQ=255)                  attatAGATGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCa                                                                                                                                    >  4:60586/1‑149 (MQ=255)                   ttatAGATGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAg                                                                                                                                   >  3:482328/1‑149 (MQ=255)                         aTGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGc                                                                                                                             >  3:402549/1‑149 (MQ=255)                           gtgtCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCa                                                                                                                            <  1:374243/148‑1 (MQ=255)                              tCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGa                                                                                                                                       >  4:839425/1‑134 (MQ=255)                              tCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGa                                                                                                                                       <  3:839425/134‑1 (MQ=255)                               ccGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGAc                                                                                                                        >  1:136/1‑148 (MQ=255)                                 gCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCaa                                                                                                                           >  4:305234/1‑143 (MQ=255)                                 gCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCaa                                                                                                                           <  3:305234/143‑1 (MQ=255)                                   cTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTc                                                                                                                                          <  4:789526/126‑1 (MQ=255)                                   cTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTc                                                                                                                                          >  3:789526/1‑126 (MQ=255)                                        ggTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAatat                                                                                                               <  2:769187/148‑1 (MQ=255)                                              tGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGAcct                                                                                                         >  2:93432/1‑147 (MQ=255)                                               gTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGaa                                                                                                       >  4:870908/1‑149 (MQ=255)                                                tAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGaa                                                                                                       >  1:685697/1‑148 (MQ=255)                                                 aacaacAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAg                                                                                                      >  1:750952/1‑148 (MQ=255)                                                 aacaacAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGt                                                                                                     >  3:44817/1‑149 (MQ=255)                                                     acaGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTAcc                                                                                                  >  2:21901/1‑148 (MQ=255)                                                     acaGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCg                                                                                                 >  3:105734/1‑149 (MQ=255)                                                      caGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCgg                                                                                                <  4:105734/149‑1 (MQ=255)                                                      caGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCg                                                                                                 >  1:813746/1‑148 (MQ=255)                                                       agagTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCggg                                                                                               >  3:787228/1‑149 (MQ=255)                                                           tGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTaa                                                                                            <  1:200230/148‑1 (MQ=255)                                                                gcgcCAATAACATCTTCTAGATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATtgtg                                                                                      >  3:551653/1‑149 (MQ=255)                                                                 cgcCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATtgtg                                                                                      >  2:877332/1‑148 (MQ=255)                                                                          cATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCt                                                                            >  4:500643/1‑149 (MQ=255)                                                                          cATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCt                                                                            >  4:703446/1‑149 (MQ=255)                                                                            tcttctAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTc                                                                          <  4:482328/149‑1 (MQ=255)                                                                               tctAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGa                                                                        >  1:174243/1‑148 (MQ=255)                                                                                     tGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGttctt                                                                  <  1:599009/148‑1 (MQ=255)                                                                                      ggACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGttcttc                                                                 >  2:29947/1‑148 (MQ=255)                                                                                            aaTTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCaa                                                                                                                           >  4:434638/1‑84 (MQ=255)                                                                                            aaTTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCaa                                                                                                                           <  3:434638/84‑1 (MQ=255)                                                                                             aTTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACt                                                          >  2:292850/1‑148 (MQ=255)                                                                                              ttAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTc                                                         >  2:69047/1‑148 (MQ=255)                                                                                              ttAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCg                                                        >  3:416730/1‑149 (MQ=255)                                                                                                  cacCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATa                                                            >  3:765749/1‑141 (MQ=255)                                                                                                  cacCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATa                                                            <  4:765749/141‑1 (MQ=255)                                                                                                   acCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCAc                                                   >  3:184008/1‑149 (MQ=255)                                                                                                   acCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCAc                                                   >  4:35856/1‑149 (MQ=255)                                                                                                    ccATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCAc                                                   >  2:549295/1‑148 (MQ=255)                                                                                                    ccATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACt                                                  >  3:348514/1‑149 (MQ=255)                                                                                                     cATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTg                                                 >  3:692095/1‑149 (MQ=255)                                                                                                       tGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGtt                                               <  3:500643/149‑1 (MQ=255)                                                                                                       tGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAACCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGt                                                <  2:750952/148‑1 (MQ=255)                                                                                                        gTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGtt                                               >  2:196444/1‑148 (MQ=255)                                                                                                             ccATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTa                                          >  2:259119/1‑148 (MQ=255)                                                                                                             ccATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTa                                          >  1:753443/1‑148 (MQ=255)                                                                                                               aTCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACt                                        <  1:21901/148‑1 (MQ=255)                                                                                                               aTCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACt                                        <  2:489604/148‑1 (MQ=255)                                                                                                               aTCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACt                                        <  2:590230/148‑1 (MQ=255)                                                                                                               aTCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTGCGAGTTCGTCTAATACTCGTCCACTGTTGGTTACt                                        <  1:858364/148‑1 (MQ=255)                                                                                                                tCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTg                                       >  2:506585/1‑148 (MQ=255)                                                                                                                 cAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGa                                     <  4:402549/149‑1 (MQ=255)                                                                                                                         gTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCgg                             <  3:423743/149‑1 (MQ=255)                                                                                                                             ttCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGggtg                         <  3:870908/149‑1 (MQ=255)                                                                                                                             ttCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGggt                          <  1:292850/148‑1 (MQ=255)                                                                                                                                 ggCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGtt                      >  1:159356/1‑148 (MQ=255)                                                                                                                                  gCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGttt                     >  1:861210/1‑148 (MQ=255)                                                                                                                                  gCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGttt                     <  1:196444/148‑1 (MQ=255)                                                                                                                                  gCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGttt                     >  2:756977/1‑148 (MQ=255)                                                                                                                                   cAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGtttc                    <  1:877332/148‑1 (MQ=255)                                                                                                                                    aaaGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGtttct                   <  1:259119/148‑1 (MQ=255)                                                                                                                                     aagaagCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGtttctt                  >  1:707113/1‑148 (MQ=255)                                                                                                                                      agaagCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGtttctttc                <  4:348514/149‑1 (MQ=255)                                                                                                                                       gaagCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGtttctttc                <  2:861210/148‑1 (MQ=255)                                                                                                                                           ggaccGGGACGTCATCCCGGACCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGt           >  3:228/4‑149 (MQ=255)                                                                                                                                            gaccGGGACGTCATCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGt           >  1:3902/3‑148 (MQ=255)                                                                                                                                             tCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGTg          >  1:874488/1‑148 (MQ=255)                                                                                                                                             tCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGTGc         >  4:243011/1‑149 (MQ=255)                                                                                                                                                  gATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGTGCCTTgc    <  3:342806/149‑1 (MQ=255)                                                                                                                                                  gATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGTGCCTTg     <  1:805786/148‑1 (MQ=255)                                                                                                                                                    tGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGGAGATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGTGCCTTgcgc  >  4:678833/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGAGCTGGCGCACATATATTATAGATGTGTCCGCCTTTAGGTGCCTGTAACAACAGAGTGATAGCGCCAATAACATCTTCTAAATGGACTAAATTAACACCATGTTCACCATCAGGCGCGGTTTTTCCGGCAAAGAAGCGTCCGGGATGACGTCCCGGTCCGACCAGGCCCGCAAGACGCAGAATATCGACCGAAGTACCGGGTAAATTGTGCAGCCAGTCTTCGAGTTCTTCTAATACTCGTCCACTGTTGGTTACTGGATTACGCGGGGTGGTTTCTTTCACCGTGCCTTGCGC  >  NC_000913/2088479‑2088774

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