BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,385,208	T→A	5.3%	V457E (GTG→GAG) ‡	yjeM →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,385,208	0	T	A	5.3%	386.8 / 7.8	132	V457E (GTG→GAG) ‡	yjeM	putative transporter
Reads supporting (aligned to +/- strand): ref base T (74/49); new base A (2/5); total (76/56)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 138 total aligned reads displayed.

CCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTATCG  >  NC_000913/4385059‑4385356

ccTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                      <  1:1516150/150‑1 (MQ=255)

 cTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgtg                                                                                                                                                     >  1:210377/1‑150 (MQ=255)

 cTCGCGCTGGCTTTCCCGTTCTATAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgtg                                                                                                                                                     >  1:1058364/1‑150 (MQ=255)

    gcgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                                  <  2:307620/150‑1 (MQ=255)

     cgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTg                                                                                                                                                 <  2:233883/150‑1 (MQ=255)

       cTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGaa                                                                                                                                               >  2:495818/1‑150 (MQ=255)

        tGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaat                                                                                                                                              <  1:578083/150‑2 (MQ=255)

        gggCTTTCCCGGTCTGTAAAGGACGTCAGGAGCTCGCCAGAGCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGGTGTCGTACTGGTGGTGACATTTGCGAACGTCGTCACCATCATTCAACCTGTGGTTGAAg                                                                                                                                              <  2:369042/149‑1 (MQ=255)

         ggCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaatg                                                                                                                                             <  1:451765/150‑3 (MQ=255)

         ggCTTACCCGTTCTTAAAAGCACGTCAGGATCTCGACAGACCGTTGGTGATTTTCAAAACGCATTAGTCGGCAATGATTGCGACAGTGGTTGTCGGACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGc                                                                                                                                             <  2:51334/150‑1 (MQ=255)

           cTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCg                                                                                                                                           <  1:47113/150‑1 (MQ=255)

               cccGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGgaga                                                                                                                                       >  2:510121/1‑150 (MQ=255)

                  gTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaat                                                                                                                                              <  2:1497558/140‑2 (MQ=255)

                    tCTTTAACGCACGGCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGGCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTggg                                                                                                                                  <  2:907228/150‑1 (MQ=255)

                    tCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTggg                                                                                                                                  >  1:623416/1‑150 (MQ=255)

                     cTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                                 >  2:1235235/1‑150 (MQ=255)

                     cTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                                 <  2:1127905/150‑1 (MQ=255)

                      tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGAc                                                                                                                                <  1:738664/150‑1 (MQ=255)

                       ttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACa                                                                                                                               >  1:95291/1‑150 (MQ=255)

                        tAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTCTGGTTGAAGCCGGAGACTGGGACac                                                                                                                              >  2:477696/1‑149 (MQ=255)

                        tAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaatgat                                                                                                                                           <  1:273455/137‑5 (MQ=255)

                        tAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaatgat                                                                                                                                           >  2:273455/1‑133 (MQ=255)

                         aaaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCCTCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGc                                                                                                                             >  2:695366/1‑150 (MQ=255)

                         aaaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGc                                                                                                                             <  1:93337/150‑1 (MQ=255)

                          aaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGca                                                                                                                            >  2:709206/1‑150 (MQ=255)

                          aaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGca                                                                                                                            >  1:2404/1‑150 (MQ=255)

                           aGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGcac                                                                                                                           >  1:887252/1‑150 (MQ=255)

                               cGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtg                                                                                                                       >  2:1643631/1‑150 (MQ=255)

                               cGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtg                                                                                                                       <  1:510121/150‑1 (MQ=255)

                                   aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGga                                                                                                                   <  1:653482/150‑1 (MQ=255)

                                    ggATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgat                                                                                                                  >  1:926882/1‑150 (MQ=255)

                                           gACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTggc                                                                                                           >  2:821027/1‑150 (MQ=255)

                                           gACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTggc                                                                                                           >  1:943132/1‑150 (MQ=255)

                                             cAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTggcgg                                                                                                         >  2:1488064/1‑150 (MQ=255)

                                             cAGACCGCTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTggcgg                                                                                                         <  1:1064425/150‑1 (MQ=255)

                                              cgACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTggcggc                                                                                                        <  2:1101559/149‑1 (MQ=255)

                                               gACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                       >  2:1530853/1‑150 (MQ=255)

                                               gACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                       >  1:934998/1‑150 (MQ=255)

                                                aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGccc                                                                                                      >  2:1337608/1‑150 (MQ=255)

                                                  cGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAATGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcccgg                                                                                                    >  2:968391/1‑148 (MQ=255)

                                                     ttGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGATGAAGACGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtct                                                                                                 <  1:821419/150‑1 (MQ=255)

                                                     ttGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGGGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtct                                                                                                 <  2:874840/150‑1 (MQ=255)

                                                        tgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttct                                                                                              <  1:148849/150‑1 (MQ=255)

                                                              tCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGc                                                                                        <  2:40139/150‑1 (MQ=255)

                                                                aaaaCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTg                                                                                      >  1:1377979/1‑150 (MQ=255)

                                                                  aaCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGtt                                                                                    >  2:881131/1‑150 (MQ=255)

                                                                   aCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCCTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTa                                                                                   >  2:300457/1‑150 (MQ=255)

                                                                    cGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAg                                                                                  >  1:264817/1‑150 (MQ=255)

                                                                      cATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCg                                                                                <  2:1182115/150‑1 (MQ=255)

                                                                       atgtggcGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGa                                                                               <  1:1375233/144‑1 (MQ=255)

                                                                        tttGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGAt                                                                              <  2:623416/150‑1 (MQ=255)

                                                                              ggCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCCTCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                        >  1:1494895/1‑150 (MQ=255)

                                                                              ggCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                        >  1:437746/1‑150 (MQ=255)

                                                                              ggCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                        <  1:507831/150‑1 (MQ=255)

                                                                                cAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGCCAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAttt                                                                      >  1:302616/1‑150 (MQ=255)

                                                                                cAATGATCGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAttt                                                                      >  1:1446786/1‑150 (MQ=255)

                                                                                    gATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCa                                                                  >  1:127785/1‑150 (MQ=255)

                                                                                       tGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGaa                                                               >  2:721413/1‑150 (MQ=255)

                                                                                        gCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAAc                                                              >  1:522110/1‑150 (MQ=255)

                                                                                              gTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGTAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgc                                                        >  1:1574081/1‑150 (MQ=255)

                                                                                               tGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                       <  1:1130042/150‑1 (MQ=255)

                                                                                                ggTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCTGAGCCTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgcag                                                      >  2:1663359/1‑150 (MQ=255)

                                                                                                ggTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgcag                                                      <  2:1324443/150‑1 (MQ=255)

                                                                                                 gTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGTATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTTTTgcagc                                                     <  1:147037/150‑1 (MQ=255)

                                                                                                   tGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGcagcag                                                   <  1:1337608/150‑1 (MQ=255)

                                                                                                   ggTCGGACTGGTGGGGACATAGGCGAACGTCTTCCCCATCATTCAACCTGGGGTTGAAGCCGGAGACTGGGACAGCACATTGTTGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCTATTTACCAGAACTATTGcagcag                                                   <  1:314514/149‑1 (MQ=255)

                                                                                                        tACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTgg                                              >  1:15773/1‑150 (MQ=255)

                                                                                                          cTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGCCAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCa                                            >  1:1124335/1‑150 (MQ=255)

                                                                                                          cTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCa                                            <  2:887252/150‑1 (MQ=255)

                                                                                                            ggtggtGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCaaa                                          <  1:502305/150‑1 (MQ=255)

                                                                                                               ggtgACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCaaaaaa                                       >  1:496241/1‑150 (MQ=255)

                                                                                                                gtgACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAAt                                      >  1:1223384/1‑150 (MQ=255)

                                                                                                                gtgACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAAt                                      >  1:1480177/1‑150 (MQ=255)

                                                                                                                  gACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATcc                                    >  1:1049901/1‑150 (MQ=255)

                                                                                                                   aCATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCg                                   >  1:383216/1‑150 (MQ=255)

                                                                                                                    cATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGc                                  <  1:881131/150‑1 (MQ=255)

                                                                                                                    cATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGACGCAGACTGGGACAGAACATTGTGGATGATTGGCGTACCTGTCTTCTTCACGCTGTTAGCGATGCCGAGTTACCAGCACTATTGCATCAGAGTTGCAAAAAATACGc                                  >  1:1124727/1‑150 (MQ=255)

                                                                                                                     aTTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCa                                 >  2:1604945/1‑150 (MQ=255)

                                                                                                                       ttGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTTGCGGCCCTGTCTTCTTCTCGCTGTTAGCTATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGt                               <  1:29183/150‑1 (MQ=255)

                                                                                                                         gCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTgg                             >  2:824734/1‑150 (MQ=255)

                                                                                                                         gCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTgg                             <  2:1478686/150‑1 (MQ=255)

                                                                                                                         gCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTgg                             >  1:1514046/1‑150 (MQ=255)

                                                                                                                          cGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTggg                            >  1:1148910/1‑150 (MQ=255)

                                                                                                                           gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGc                           <  2:496241/150‑1 (MQ=255)

                                                                                                                           gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGc                           >  2:1282187/1‑150 (MQ=255)

                                                                                                                           gAACGTCTTCACCATCATTCAAACTGTGGTTGAAGCCGGAGACTGGGACAGCACATTTTGGATTATTGGAGGCCCTGTCTTCTTCTCGATGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGc                           <  1:1318509/150‑1 (MQ=255)

                                                                                                                              cGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGt                        <  2:934998/150‑1 (MQ=255)

                                                                                                                              cGTCATCACCATCATACAACCTGTGGTTGAAGCCCGAGACTGGGACAGCACAATGTGGATGATTGGCGGCCCTTTCTTCTTCTCGCTGTTAGCGATGTCGATTTACAAGAACTATTGCCGAAGAGTGGCAATAAATCCCCAGTGGGAGGt                        >  1:784822/1‑150 (MQ=255)

                                                                                                                                tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTgg                      <  2:95291/150‑1 (MQ=255)

                                                                                                                                 cTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTAGGCGGTGGa                     <  2:926882/150‑1 (MQ=255)

                                                                                                                                    cACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATa                  >  1:943187/1‑150 (MQ=255)

                                                                                                                                      ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc                >  2:644418/1‑150 (MQ=255)

                                                                                                                                      ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc                >  2:650689/1‑150 (MQ=255)

                                                                                                                                      ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc                >  1:1514590/1‑150 (MQ=255)

                                                                                                                                              cAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGcc        >  1:1287920/1‑150 (MQ=255)

                                                                                                                                              cAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGcc        >  2:232755/1‑150 (MQ=255)

                                                                                                                                               aaCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGccc       >  1:705077/1‑150 (MQ=255)

                                                                                                                                                aCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCt      >  2:1502312/1‑150 (MQ=255)

                                                                                                                                                  cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTAt    >  1:1712050/1‑150 (MQ=255)

                                                                                                                                                   tgtgGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTATc   <  2:1446786/150‑1 (MQ=255)

CCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTATCG  >  NC_000913/4385059‑4385356

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