Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,411,047	C→T	S165S (AGC→AGT)	yjfJ →	PspA/IM30 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,411,047	0	C	T	98.8%	276.6 / ‑4.7	83	S165S (AGC→AGT)	yjfJ	PspA/IM30 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (46/36);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGC  >  NC_000913/4410951‑4411142                                                                                                 |                                                                                                tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGtt                                                                                              <  1:2556068/100‑1 (MQ=255)  cAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGtt                                                                                              <  2:495999/99‑1 (MQ=255)    gTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCg                                                                                           <  2:990922/100‑1 (MQ=255)    gTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCg                                                                                           <  1:2750221/100‑1 (MQ=255)    gTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCg                                                                                           <  1:1639346/100‑1 (MQ=255)     tttGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGa                                                                                          <  1:419573/100‑1 (MQ=255)       tGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACa                                                                                        >  2:2150867/1‑100 (MQ=255)        gAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAg                                                                                       >  1:243932/1‑100 (MQ=255)         agcagcAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGc                                                                                      <  2:166375/100‑1 (MQ=255)         agcagcAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGc                                                                                      <  2:2243862/100‑1 (MQ=255)         agcagcAAATGGAAGTCGTTAAAGCCACTGAAGACATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTCTTTCGACAGc                                                                                      >  1:1906148/1‑100 (MQ=255)           cagcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCgg                                                                                    >  2:1024925/1‑100 (MQ=255)           cagcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCgg                                                                                    >  1:2798910/1‑100 (MQ=255)           cagcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCgg                                                                                    <  2:744279/100‑1 (MQ=255)            agcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCgg                                                                                    <  1:1772007/99‑1 (MQ=255)            agcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc                                                                                   >  1:984991/1‑100 (MQ=255)             gcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCa                                                                                  >  1:971189/1‑100 (MQ=255)             gcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCa                                                                                  >  1:496453/1‑100 (MQ=255)             gcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCa                                                                                  <  2:1996744/100‑1 (MQ=255)             gcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCa                                                                                  >  1:1352973/1‑100 (MQ=255)             gcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCa                                                                                  <  2:2751382/100‑1 (MQ=255)               aaaTGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGa                                                                                >  2:1121186/1‑100 (MQ=255)                aaTGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGaa                                                                               >  2:2406002/1‑100 (MQ=255)                aaTGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGaa                                                                               >  2:2406178/1‑100 (MQ=255)                 aTGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGaa                                                                               >  2:2814676/1‑99 (MQ=255)                     aaGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCtt                                                                          >  2:1907545/1‑100 (MQ=255)                     aaGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCtt                                                                          >  2:1803946/1‑100 (MQ=255)                     aaGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCtt                                                                          >  1:553610/1‑100 (MQ=255)                       gTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTaa                                                                        <  1:2555524/100‑1 (MQ=255)                         cGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTcaaa                                                                      >  1:2603218/1‑100 (MQ=255)                         cGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTaaaa                                                                      <  2:1346716/100‑1 (MQ=255)                         cGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTaaaa                                                                      >  1:1017270/1‑100 (MQ=255)                          gTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAAc                                                                     <  2:1785214/100‑1 (MQ=255)                            tAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGc                                                                   >  2:1786813/1‑100 (MQ=255)                            tAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGc                                                                   >  2:4871/1‑100 (MQ=255)                             aaaGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGcc                                                                  >  2:1785813/1‑100 (MQ=255)                                gCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGc                                                               >  1:2091320/1‑100 (MQ=255)                                 ccACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCa                                                              <  2:2719295/100‑1 (MQ=255)                                  cACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaa                                                             <  2:1774263/100‑1 (MQ=255)                                  cACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaa                                                             >  2:1194725/1‑100 (MQ=255)                                  cACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaa                                                             >  1:2890927/1‑100 (MQ=255)                                   aCTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaaa                                                            >  2:623228/1‑100 (MQ=255)                                   aCTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaaa                                                            >  1:290680/1‑100 (MQ=255)                                     tGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAcg                                                          <  2:277356/100‑1 (MQ=255)                                     ngAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAcg                                                          >  1:2100222/2‑100 (MQ=255)                                        aGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGt                                                       >  1:188893/1‑100 (MQ=255)                                         gCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTc                                                      <  1:232069/100‑1 (MQ=255)                                           cATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAg                                                    >  1:1092571/1‑100 (MQ=255)                                           cATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAg                                                    >  2:2246157/1‑100 (MQ=255)                                              gCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGcc                                                 <  2:1477712/100‑1 (MQ=255)                                                  cGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAAc                                             <  2:2415119/100‑1 (MQ=255)                                                      cacaACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCa                                         >  1:954591/1‑100 (MQ=255)                                                      cacaACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCa                                         <  2:2344533/100‑1 (MQ=255)                                                      cacaACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCa                                         <  2:2305932/100‑1 (MQ=255)                                                       acaacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAg                                        >  2:112430/1‑100 (MQ=255)                                                       acaacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAg                                        <  1:708182/100‑1 (MQ=255)                                                        caacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAgg                                       >  2:1118659/1‑100 (MQ=255)                                                        caacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAg                                        >  1:1127379/1‑99 (MQ=255)                                                         aacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGc                                      <  2:2798946/100‑1 (MQ=255)                                                          acaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCt                                     >  2:494055/1‑100 (MQ=255)                                                          acaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCt                                     >  1:1371342/1‑100 (MQ=255)                                                          acaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCt                                     >  2:2076595/1‑100 (MQ=255)                                                           caGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTc                                    >  2:163657/1‑100 (MQ=255)                                                               cTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCt                                >  1:2559719/1‑100 (MQ=255)                                                                 gTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCAACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTgg                              >  1:1315212/1‑100 (MQ=255)                                                                  tAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGa                             <  2:1684729/100‑1 (MQ=255)                                                                  tAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGa                             >  1:1409849/1‑100 (MQ=255)                                                                   aacaacCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGAt                            <  1:2150840/100‑1 (MQ=255)                                                                     caacCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGAtgc                          >  1:404513/1‑100 (MQ=255)                                                                     caacCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGAtgc                          <  2:681822/100‑1 (MQ=255)                                                                            tACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGca                   >  2:2030835/1‑100 (MQ=255)                                                                             aCCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGcac                  <  1:1301662/100‑1 (MQ=255)                                                                               cGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAg                >  2:336705/1‑100 (MQ=255)                                                                                gTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGt               <  2:164230/100‑1 (MQ=255)                                                                                  tGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTg             <  2:2100250/100‑1 (MQ=255)                                                                                   ggCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTgg            <  1:2548153/100‑1 (MQ=255)                                                                                    gcgcATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGga           >  1:418159/1‑100 (MQ=255)                                                                                     cgcATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGgag          >  2:1858314/1‑100 (MQ=255)                                                                                     cgcATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGgag          <  1:2829810/100‑1 (MQ=255)                                                                                     cgcATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGgag          <  2:2091348/100‑1 (MQ=255)                                                                                         tCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAg      <  2:404518/100‑1 (MQ=255)                                                                                            tCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCg   >  1:2741848/1‑100 (MQ=255)                                                                                             ccAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGc  <  2:609318/100‑1 (MQ=255)                                                                                                 |                                                                                                TCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGC  >  NC_000913/4410951‑4411142

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