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

ACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCT  >  NC_000913/4622101‑4622399                                                                                                                                                        |                                                                                                                                                      aCCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCggggag                                                                                                                                                      <  1:1075641/151‑1 (MQ=255)   cGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCggggaggg                                                                                                                                                    <  1:1207773/151‑1 (MQ=255)        ggCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑GT‑TGTTGGGCgggggggggg                                                                                                                                                  >  1:972041/1‑147 (MQ=255)              ttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGTTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAg                                                                                                                                          <  2:702694/150‑1 (MQ=255)               tAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGGCGGGGGGGGGGATTATTAGc                                                                                                                                         >  2:1100372/1‑150 (MQ=255)                aGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCg                                                                                                                                        <  1:849890/150‑1 (MQ=255)                  ggCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGGGGGGGATTATCAGCGATg                                                                                                                                     >  2:304962/1‑151 (MQ=255)                   gcgcgcGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATgg                                                                                                                                    <  2:1115787/151‑1 (MQ=255)                   gcgcgcGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATgg                                                                                                                                    <  2:64539/151‑1 (MQ=255)                    cgcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGGCGGGGGGGGGGAATATCAGCGATgg                                                                                                                                    >  2:242091/1‑150 (MQ=255)                      cgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGAtt                                                                                                                                 <  2:990404/151‑1 (MQ=255)                        cgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc                                                                                                                               <  1:20333/151‑1 (MQ=255)                           ggggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTg                                                                                                                             <  1:644555/150‑1 (MQ=255)                            gggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑TGGTGGGGGGGGGGGGGGATTATCAGCGGTGGGATACTg                                                                                                                             >  2:335185/1‑149 (MQ=255)                               cAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGTCGGGGAGGGGAATTATCAACGATGGATTACTGcac                                                                                                                          >  2:265165/1‑149 (MQ=255)                                aGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGGGGGGGGGGGGGGTTATCAGCGATGGGTTAATGCacaa                                                                                                                        >  1:788060/1‑150 (MQ=255)                                 ggAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaaca                                                                                                                      <  1:284485/151‑1 (MQ=255)                                 ggAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑GT‑GGGTGGTCGGGGGGGGGGATTATTAGGGATGGGTTACTGCacaa                                                                                                                        >  1:440623/1‑149 (MQ=255)                                    aGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGGTGGGGGGGGGAATTATATACGGAGGAGTACTGcagcaaaccc                                                                                                                   >  2:78971/1‑143 (MQ=255)                                       tAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGGGGGGGATTATTAGCGATGGGTTTCTGCACAAAAccc                                                                                                                   >  2:865980/1‑148 (MQ=255)                                       tAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGGGGGGAATTATCAGCGAGGGATTACTGCACAACACCCTGc                                                                                                                >  1:1018718/1‑151 (MQ=255)                                           cgcgCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAAAACCCCGCgg                                                                                                              >  2:904821/1‑149 (MQ=255)                                                    gAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCctct                                                                                                    <  1:759833/150‑1 (MQ=255)                                                    gAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTg                                                                                                   <  1:247609/151‑1 (MQ=255)                                                    gAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGGCGGGGGGGGGGATTATCAGCGATGGATTAATGCACAACACCCTGCGGGAATTtctc                                                                                                     >  1:381229/1‑149 (MQ=255)                                                      aTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGGCGGGGGGGGGGAATATTAGCGGAGGATTTCTTCCCAACACCCTGCGGGAAATTCTCTgg                                                                                                  >  1:173542/1‑150 (MQ=255)                                                        ctgactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGCTT‑GGTT‑GTCGGGGGGGGGGATTATTAGCGAAGGGTTACTCCAAAAAACCCTGCGGAAATAcc                                                                                                       >  1:113719/1‑143 (MQ=255)                                                         tgactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTCCTCTGGATc                                                                                               >  2:218428/1‑150 (MQ=255)                                                           actgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGGTTGGGCGGGGGGGGGAA‑TATTAGCGCTAGGGTTCTCAACAAAACCCAGAGGAAATTtctc                                                                                                     >  2:491723/1‑142 (MQ=255)                                                                 gatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTc                                                                                       <  2:183565/150‑1 (MQ=255)                                                                     taCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAccc                                                                                   <  2:349716/150‑1 (MQ=255)                                                                      aCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCAt                                                                                 <  2:803270/151‑1 (MQ=255)                                                                              ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCCCCCATTACGa                                                                            <  2:476608/148‑1 (MQ=255)                                                                                cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGGCGGGGGGGGGAATTATCAGCGATGGATTTATGCACAACACCCTGCGGGAATTTCTCTGGATCAAAAAATCACCCCTTTCGAAc                                                                          >  2:999226/1‑148 (MQ=255)                                                                                cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑GT‑GGGTTGGCGGGGGGGGGGAATATCAGCGGTGGGTTACTGCCCAAAAACCCCCGGGAAATCCCCTGGGTTAAAAAATCAccc                                                                                   >  1:101572/1‑139 (MQ=255)                                                                                 aGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCt                                                                      <  2:447181/151‑1 (MQ=255)                                                                                  ggAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTg                                                                     <  1:990483/151‑1 (MQ=255)                                                                                       aagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGTCGGGGGGGGGAATTATCAGCGATGGGTTACTGCACCACACCCCGCGGGAATTCCTCTGGATCAAAAAATCACCCATTACGAACAGCTGGCa                                                                  >  2:6269/1‑149 (MQ=255)                                                                                        agaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTTCTGCACAACACCCTGCGGGAATTCCTCTGGATCAAAAACTCACCCCTTACGAACAGCTgg                                                                    >  2:395640/1‑146 (MQ=255)                                                                                          aTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTg                                                              <  1:251773/150‑1 (MQ=255)                                                                                          aTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGa                                                             <  2:445075/151‑1 (MQ=255)                                                                                              gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCa                                                         <  2:667219/151‑1 (MQ=255)                                                                                                  tACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑GT‑GGGAGGGCGGGGAGGGGGAATATCCGCGATGGGTTTCTGCTCAACCCCCTGTGGGAAATCCTCCGGATCAAAAAATCAGCCCTTCCGAACAACGGGCAAGCGAAGc                                                          >  2:1212969/1‑146 (MQ=255)                                                                                                   aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACt                                                       <  2:59910/148‑1 (MQ=255)                                                                                                   aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGGTGGTCGGGGGGGGGAATTATCAGCGAAGGGTTACTGCACAACACCCCGCGGGAATTTCTCTTGAGCAAAAAAACACCCATTACGAAAAGCTGGCAAGGGATGCACCgg                                                     >  1:859674/1‑150 (MQ=255)                                                                                                   aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑GT‑GGGTGGTTGGGGGGGGGGATTATCAGCGATGGATTACTGCACAAAACCCTGCGGCAATTTCTCTGGGCCAAAAACTCACCCATTAAGAAAAGCAGGCCAGGGATGCACCGGc                                                    >  2:811877/1‑151 (MQ=255)                                                                                                     gCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCCGGGAGGGGAAATATTAAGGATGGATTTCTGCCCAACAACCTCCGGGAAATCCTCTGGATCAAAAAATCACCCATTACGGACACATGGCAAGGGATGCACCGGc                                                    >  2:454664/1‑149 (MQ=255)                                                                                                       ccTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑TGTTGGCTGGGGGGGGGGATTATCAGCCATAGCTTTCTTGACAACACCCCGCGGAAATTTCTCTGGGTCAAAAAAACAACCCATACGAACAGCTGGCAAATGATGcccccgca                                                   >  2:14288/1‑141 (MQ=255)                                                                                                          gACCGTATTTAACGGTGAATATACTGGCGG‑TT‑TGGTGGTGGGGGGGGGGGAATATTAGAGATGGGTTAATTCACAACACCCCGCGGGAAATTCTCTTGGGCAAAAAATCACCCCTTACGAACAGCTGGCAAGTGGTGCACTgggcgg                                                 >  1:508893/1‑143 (MQ=255)                                                                                                          gACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGgag                                               <  1:512748/149‑1 (MQ=255)                                                                                                              gTATTTAACGGTGAATATACTGGCGG‑TT‑AGGTGGTCGGGGGGGGGGATTATCAGCGACGGGTTACTGCACCACACCCTCCGGAAATTTCTCTGGAACCAAAAATCACCCATTACGAACAACTGGGAAGGGATGAACACGCAGGCGaaa                                            >  2:586618/1‑148 (MQ=255)                                                                                                                   tAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGGGGGGAATTATTAGCGATGGATTACTGCACAACACCCTGCGGGAAATCCTCTGGGTCAAAAAATCACCCCTTACGAACAGCTGGCAAGTGATGCACTGGCAGGGGAAAATTTggg                                     >  2:988646/1‑150 (MQ=255)                                                                                                                   tAACGGTGAATATACTGGCGG‑TT‑GGGTGGTCGGGGGGGGGAATTATTAGCGATGGGTTACTGCCCAACACCCCGCGGGAATTCCTCTGGGTCAAAAAATCACCCATTACGAACCGCTGGCAAATGATGCACTCGCAGGGGGAAATGTggg                                     >  2:14724/1‑150 (MQ=255)                                                                                                                     aCGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGttc                                  <  1:909719/151‑1 (MQ=255)                                                                                                                        gTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGttcttc                               <  1:183328/151‑1 (MQ=255)                                                                                                                        gTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGttcttc                               <  2:378303/151‑1 (MQ=255)                                                                                                                           aaTATACTGGCGG‑TT‑GGTTGGGCGGGGGGGGGGATTATCAGCGATTGGTTTCTGCACAAAACCCCGCGGGAATTCCTCTGGGTTAAAAAATCACCCATTTCGAAAAACTGGCAAGGGGGGGACTGGAAGGGGGAAATGTGGGgtctttt                              >  1:942676/1‑147 (MQ=255)                                                                                                                              ataCTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcggg                         <  1:291750/151‑1 (MQ=255)                                                                                                                              ataCTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgg                          <  2:41505/150‑1 (MQ=255)                                                                                                                              ataCTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgg                          <  1:243279/150‑1 (MQ=255)                                                                                                                                   ggCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGggcggc                     <  2:1198355/150‑1 (MQ=255)                                                                                                                                   ggCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGgcggcg                    <  1:1218183/151‑1 (MQ=255)                                                                                                                                               tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCaa           >  2:88973/1‑150 (MQ=255)                                                                                                                                                ggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCaaa          <  1:812926/150‑1 (MQ=255)                                                                                                                                                 gtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAt        <  2:381241/151‑1 (MQ=255)                                                                                                                                                  tCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAtt       >  2:824542/1‑151 (MQ=255)                                                                                                                                                   cGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAttt      <  2:489709/151‑1 (MQ=255)                                                                                                                                                     gggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTa     >  1:395870/1‑150 (MQ=255)                                                                                                                                                       gaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACct  >  2:285486/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      ACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TT‑GGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCT  >  NC_000913/4622101‑4622399

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