Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,581,765	T→C	8.2%	D338G (GAC→GGC)	ypfG ←	DUF1176 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,581,765	0	T	C	8.2%	293.8 / 11.9	97	D338G (GAC→GGC)	ypfG	DUF1176 family protein
Reads supporting (aligned to +/- strand):  ref base T (41/48);  new base C (5/3);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATC  >  NC_000913/2581618‑2581913                                                                                                                                                    |                                                                                                                                                     gCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTc                                                                                                                                                     <  1:1097136/149‑1 (MQ=255)   tGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                   >  3:236963/1‑149 (MQ=255)    gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCAACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGTCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                                  <  3:93534/149‑1 (MQ=255)           aaaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                                                                           <  1:492095/149‑1 (MQ=255)            aaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGcc                                                                                                                                          >  4:722432/1‑149 (MQ=255)             aGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCa                                                                                                                                         <  4:548851/149‑1 (MQ=255)             aGGTGCCTGATGGTGATTGACGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCa                                                                                                                                         >  3:636322/1‑149 (MQ=255)              ggTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgttgt                                                                                                                                        <  3:800267/149‑7 (MQ=255)              ggTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACATGGGCCAGACGCCTGGCCCacgttgt                                                                                                                                        <  1:214832/149‑7 (MQ=39)               gTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGt                                                                                                                                       <  4:423852/149‑1 (MQ=255)               gTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGt                                                                                                                                       >  1:536130/1‑149 (MQ=255)                  ccTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTAt                                                                                                                                    <  4:679869/149‑1 (MQ=255)                  ccTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTAt                                                                                                                                    <  2:232392/149‑1 (MQ=255)                    tGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCg                                                                                                                                  >  3:721578/1‑149 (MQ=255)                    tGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCg                                                                                                                                  >  1:830098/1‑149 (MQ=255)                    tGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCg                                                                                                                                  >  2:1081462/1‑149 (MQ=255)                     gATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGc                                                                                                                                 >  3:215818/1‑149 (MQ=255)                       tGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                                                                               >  4:814879/1‑149 (MQ=255)                         gTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                   >  1:1076267/1‑127 (MQ=255)                         gTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                   <  2:1076267/127‑1 (MQ=255)                           gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                           <  1:782933/149‑1 (MQ=255)                               gccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTc                                                                                                                       <  2:962007/149‑1 (MQ=255)                                   gatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcg                                                                                                                   <  4:636322/149‑1 (MQ=255)                                        tgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTa                                                                                                              >  2:1015681/1‑149 (MQ=255)                                         gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                             <  4:843473/149‑1 (MQ=255)                                         gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                             >  2:234301/1‑149 (MQ=255)                                          attgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAAc                                                                                                            <  3:320461/149‑1 (MQ=255)                                           ttgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACg                                                                                                           <  4:764979/149‑1 (MQ=255)                                             gccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCa                                                                                                         <  4:394235/149‑1 (MQ=255)                                                      gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                                                                                >  3:1013919/1‑149 (MQ=255)                                                         gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAAc                                                                                             <  1:1081462/149‑1 (MQ=255)                                                           tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                           <  1:1048774/149‑1 (MQ=255)                                                           tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                           >  1:853696/1‑149 (MQ=255)                                                           tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                           <  4:236963/149‑1 (MQ=255)                                                               gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGAc                                                                                       >  2:1065741/1‑149 (MQ=255)                                                                  tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                    >  1:800371/1‑149 (MQ=255)                                                                  tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTGATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                    >  2:392514/1‑149 (MQ=255)                                                                    ttATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCa                                                                                  >  3:641639/1‑149 (MQ=255)                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                                                                     <  1:8298/127‑1 (MQ=255)                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                                                                     >  2:8298/1‑127 (MQ=255)                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                                                                               <  3:722432/149‑1 (MQ=255)                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                                                                               <  1:592227/149‑1 (MQ=255)                                                                        ccGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgc                                                                              >  3:443695/1‑149 (MQ=255)                                                                         cGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcg                                                                             >  1:505419/1‑149 (MQ=255)                                                                           gCCTACAGCGTCTTGTCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa                                                                                        >  3:84159/1‑136 (MQ=255)                                                                           gCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa                                                                                        <  4:84159/136‑1 (MQ=255)                                                                           gCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGcc                                                                           <  3:346985/149‑1 (MQ=255)                                                                                  gCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                             <  3:543748/108‑1 (MQ=255)                                                                                  gCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                             >  4:543748/1‑108 (MQ=255)                                                                                       ttATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                   >  4:294913/1‑81 (MQ=255)                                                                                       ttATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                   <  3:294913/81‑1 (MQ=255)                                                                                        tATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCGCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAAt                                                              <  2:505419/149‑1 (MQ=255)                                                                                         aTCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATg                                                             >  4:128628/1‑149 (MQ=255)                                                                                          tCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGc                                                            <  4:1015471/149‑1 (MQ=255)                                                                                           ccGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGcc                                                           <  3:278145/149‑1 (MQ=255)                                                                                           ccGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGcc                                                           <  4:640933/149‑1 (MQ=255)                                                                                                 tACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCCGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATc                                                     <  1:833197/149‑1 (MQ=255)                                                                                                  aCATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCg                                                    >  4:257092/1‑149 (MQ=255)                                                                                                   cATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                   <  3:1039588/149‑1 (MQ=255)                                                                                                    aTGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGTCGTCTGGCCCATGACAGTTAACGCAGGTGGGTTCTGCGGCGTAACGCACCATGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCt                                                  <  4:46348/149‑1 (MQ=255)                                                                                                    aTGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCt                                                  >  3:594768/1‑149 (MQ=255)                                                                                                    aTGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCt                                                  >  1:175665/1‑149 (MQ=255)                                                                                                    aTGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCt                                                  <  4:1032571/149‑1 (MQ=255)                                                                                                    aTGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCt                                                  <  4:279486/149‑1 (MQ=255)                                                                                                    aTGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCt                                                  <  1:833719/149‑1 (MQ=255)                                                                                                        cccGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAAt                                              <  3:468695/149‑1 (MQ=255)                                                                                                         ccGCCATTTTGTTTAACGGGGGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAAACGCTTAATc                                             >  2:11042/1‑149 (MQ=255)                                                                                                          cGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                            >  1:41604/1‑149 (MQ=255)                                                                                                           gCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa                                                                                        >  4:610195/1‑104 (MQ=255)                                                                                                           gCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa                                                                                        <  3:610195/104‑1 (MQ=255)                                                                                                            ccATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCcgc                                          <  4:443695/149‑1 (MQ=255)                                                                                                                    tttAACGGGTGATCCACAACGTGGGCCAGGCGTTTGGCCCATGCCAGTTATCGCAGGTGGGCTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGGCAAAGCGCCACCGCGCGTGAATGCCACAATCGCTCAATCCGCGCCCTTTc                                  >  2:369642/1‑149 (MQ=255)                                                                                                                       aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGc                                                            <  1:134439/120‑1 (MQ=255)                                                                                                                       aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGc                                                            >  2:134439/1‑120 (MQ=255)                                                                                                                       aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCt                               >  2:452806/1‑149 (MQ=255)                                                                                                                         cGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTaa                             <  4:364731/149‑1 (MQ=255)                                                                                                                          gggTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAg                            >  2:565385/1‑149 (MQ=255)                                                                                                                               aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                       <  2:800371/149‑1 (MQ=255)                                                                                                                                 ccACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCa                     >  3:298479/1‑149 (MQ=255)                                                                                                                                   acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGt                   >  2:35870/1‑149 (MQ=255)                                                                                                                                        atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTaa                             >  1:533355/2‑134 (MQ=255)                                                                                                                                        atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTaa                             <  2:533355/133‑1 (MQ=255)                                                                                                                                        atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg              >  1:956489/2‑149 (MQ=255)                                                                                                                                        atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg              >  1:28916/2‑149 (MQ=255)                                                                                                                                         tGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcgc             <  4:298479/149‑1 (MQ=255)                                                                                                                                         tGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcgc             >  3:904384/1‑149 (MQ=255)                                                                                                                                          gggCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcgcg            >  3:12979/1‑149 (MQ=255)                                                                                                                                              cAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGAttt        <  3:653106/149‑1 (MQ=255)                                                                                                                                               aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                       <  2:400382/133‑1 (MQ=255)                                                                                                                                               aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                       >  1:400382/1‑133 (MQ=255)                                                                                                                                               aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtc       <  4:736631/149‑1 (MQ=255)                                                                                                                                               aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCACCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtc       >  1:641406/1‑149 (MQ=255)                                                                                                                                                 gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtctc     <  1:108163/149‑1 (MQ=255)                                                                                                                                                 gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtctc     <  2:810315/149‑1 (MQ=255)                                                                                                                                                  cGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGCAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCtca    <  4:885361/149‑1 (MQ=255)                                                                                                                                                   gTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCtcat   >  1:14364/1‑149 (MQ=255)                                                                                                                                                   gTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCtcat   >  4:284884/1‑149 (MQ=255)                                                                                                                                                    tCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCTCCAGTTCACGCGATTTCtcatc  <  4:1013919/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCATC  >  NC_000913/2581618‑2581913

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