Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,250	0	A	G	42.1%	78.6 / 18.3	57	E161G (GAG→GGG)	yjjJ	putative protein kinase
Reads supporting (aligned to +/- strand):  ref base A (7/26);  new base G (24/0);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 3.60e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTA  >  NC_000913/4622103‑4622400                                                                                                                                                    |                                                                                                                                                       cGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGATGGGTGGGCgggggggg                                                                                                                                                     >  1:806811/1‑151 (MQ=255)  gCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggagggg                                                                                                                                                    <  2:263757/151‑1 (MQ=255)  gCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggagggg                                                                                                                                                    <  2:130206/151‑1 (MQ=255)      ggCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAAtt                                                                                                                                                <  1:538199/151‑1 (MQ=255)           tttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCa                                                                                                                                            <  1:593338/150‑1 (MQ=255)            ttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATctg                                                                                                                                           >  1:808651/1‑148 (MQ=255)            ttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAg                                                                                                                                           <  2:422629/150‑1 (MQ=255)              aGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCGTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCg                                                                                                                                         <  2:614901/150‑1 (MQ=255)                 gcgcgcGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg                                                                                                                                     <  2:806827/151‑1 (MQ=255)                   gcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGAt                                                                                                                                   <  2:871253/151‑1 (MQ=255)                        tGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACt                                                                                                                               >  2:605311/1‑150 (MQ=255)                         ggggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGGATTATCAGCGATGGATTACTGc                                                                                                                             >  2:961836/1‑151 (MQ=255)                             cAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGAATTTCAGCGATGGATTACTGcaca                                                                                                                          >  1:558805/1‑150 (MQ=255)                              aGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCacaa                                                                                                                         >  1:323250/1‑150 (MQ=255)                                 aaGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATTAGCGATGGGTTACTGGacaacaac                                                                                                                     >  1:471809/1‑149 (MQ=255)                                      aGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGGTTACTGCCCAACACCCCGc                                                                                                                 >  2:50417/1‑150 (MQ=255)                                       gCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTTGTCGGGGGGGGGGATTATTCGCGATGGGTTACTGGACAACAACCTGCgg                                                                                                               >  1:951462/1‑151 (MQ=255)                                        ccGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGa                                                                                                              <  2:962663/151‑1 (MQ=255)                                              aaCTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTATCAGCGATGGATTATTGCACAAAACCCCGCGGgaa                                                                                                            >  1:131586/1‑147 (MQ=255)                                              aaCTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGGAATTcc                                                                                                        >  2:340814/1‑151 (MQ=255)                                                  gAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGTCGGGGGGGGGGATTATCAGCGGTGGATTACTGCACAAAACCCCGCGGAAATTCctct                                                                                                     >  1:663001/1‑150 (MQ=255)                                                  gAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTg                                                                                                    <  2:464426/151‑1 (MQ=255)                                                    aTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAAGGAAGGGTTATTGCACAACACCCTGCGGGAAATCctcctg                                                                                                   >  1:637577/1‑147 (MQ=255)                                                     tCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATTAGCGATGGATTAATGCACAAAACCCTGCGGAAAATCctct                                                                                                     >  1:60045/1‑147 (MQ=255)                                                     tCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATTAGCGATGGATTACTGCACAACACCCTGCGGGAATTCCTCTgg                                                                                                   >  2:496912/1‑149 (MQ=255)                                                        gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGa                                                                                                  <  2:341902/147‑1 (MQ=255)                                                        gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaa                                                                                              <  2:945013/151‑1 (MQ=255)                                                             atgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACt                                                                                         <  2:161175/151‑1 (MQ=255)                                                               gatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCa                                                                                       <  1:579803/151‑1 (MQ=255)                                                                  ataCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTGGTCGGGGAGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAAATCAccc                                                                                    >  2:33870/1‑151 (MQ=255)                                                                      gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTa                                                                                <  1:860252/151‑1 (MQ=255)                                                                      gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGAGGGGGAATTTCAGCGGTGGATTTCTGCACAACACCCCGCCGAAATTCCTCTGGATCAAAAAATCACCCAtt                                                                                 >  2:239609/1‑150 (MQ=255)                                                                        ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTTCTCTGGGTCAAAAACTCACCCATTAc                                                                               >  2:322915/1‑150 (MQ=255)                                                                         tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATTAGCGAAGGATTACTGGACAACACCCCGCGGAAATTTCTCTGGATCAAAAAATCACCCCTTa                                                                                >  1:879601/1‑148 (MQ=255)                                                                                 gaagaaGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGGTGGATTACTGCACAACACCCTGCGGGAATTCCTCTGGGTCAAAAACTCACCCATTACGAACAGCTgg                                                                     >  2:295123/1‑151 (MQ=255)                                                                                            gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTGGGCGGGGGGGGGGATTATTAGCGATTGATTTATTGACAAAACCCTTCGGGAAATCCTCTTGATCAAAAAAACTCCCCTTACGAACAACTTGCCAGAGAt                                                             >  1:547492/1‑148 (MQ=255)                                                                                                tACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTCCTCTGGATCAAAAACTCACCCCTTACGAACAGCTGGCAAGTGATGCACTgg                                                      >  2:638066/1‑151 (MQ=255)                                                                                                  cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTgg                                                      >  2:79673/1‑149 (MQ=255)                                                                                                     ccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGga                                                 <  2:29422/151‑1 (MQ=255)                                                                                                     ccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGga                                                 <  2:693387/151‑1 (MQ=255)                                                                                                             tATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAAAACCCCGCGGGAATTCCTCTGGATCAAAAAATCACCCATTACGAACAACTGGCAAGTGATGCACTGGCAGGAGAAAtt                                           >  1:783626/1‑149 (MQ=255)                                                                                                              aTTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGAGGGGAATTTTCAGCGATGGGTTACTGCACAACACCCTGCGGGAAATCCTCTGGATCAAAAACTCACCCATTACGAAAAGCTGGCAAGTGATGCACTGGCAGGGGGAAAtgt                                         >  2:114314/1‑150 (MQ=255)                                                                                                               tttAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATTAGGGGGGGATTACTTGCCAAAACCCTTCGGAAAATCCTTTTGTGCAAAAAAACACCCCATACGAACAGCTGGTAAAGGGTGGACTGGCAGGGGAAAtt                                           >  1:884216/1‑147 (MQ=255)                                                                                                                ttAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGAAAATCCTCTTGATCAAAAAAACACCCCTTACCAACAACTGGCAAGTGATGCACTTGCAGGAGAAATTTTggg                                      >  1:518069/1‑151 (MQ=255)                                                                                                                  aaCGGTGAATATACTGGCGGTTTGTTGGTCGGGGGGGGGGATTATTAGCGATGGGTTACTGCACAACACCCTGCGGGAATTTCTCTGGATCAAAAAATCACCCCCTACGAACAGCTGGCAAATGATGCACCGGCCGGAGGAAATTTggg                                      >  1:912440/1‑149 (MQ=255)                                                                                                                  aaCGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGGGATGGATTACTGCCCAACACCCTGCGGAAATTTCTCTGGATCAAAAACTCACCCCTTACGAACAGCTGGCAAGTGATGCACCGGCAGGGGGAATTGTggggt                                    >  1:642669/1‑149 (MQ=255)                                                                                                                  aaCGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGt                                     <  2:879619/150‑1 (MQ=255)                                                                                                                  aaCGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGt                                     <  2:558816/150‑1 (MQ=255)                                                                                                                  aaCGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATTAGCGATGGATTACTGCACAACACCCCGCGGGAATTTCTCTTGGTCAAAAAACCACCCATTACCAACACCGGGCAAATGATGCACTTGCAGGGGAAATTGTggg                                      >  1:286058/1‑149 (MQ=255)                                                                                                                        gAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTtcttct                               <  1:696876/150‑1 (MQ=255)                                                                                                                            ataCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgg                           <  2:471819/150‑1 (MQ=255)                                                                                                                              aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgggcg                        <  1:727193/151‑1 (MQ=255)                                                                                                                                 ggCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgggcgg                       <  2:60043/149‑1 (MQ=255)                                                                                                                                 ggCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGggcggc                      <  2:710565/150‑1 (MQ=255)                                                                                                                                  gCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGggcggc                      >  2:780556/1‑149 (MQ=255)                                                                                                                                  gCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGgcggcg                     >  1:407066/1‑150 (MQ=255)                                                                                                                                        ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGagcagc               <  1:656302/150‑1 (MQ=255)                                                                                                                                         gttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCa             <  1:295119/151‑1 (MQ=255)                                                                                                                                            ggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCaaa           <  1:340807/150‑1 (MQ=255)                                                                                                                                               cGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAttt       <  1:638051/151‑1 (MQ=255)                                                                                                                                                    aggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACcta  <  1:961815/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       CGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTA  >  NC_000913/4622103‑4622400

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